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Iht  pi'obltmt  )'«Ut(Dil  to  undi^ryruund  rApU!  OMCAVutton  <ind  their 
reiivihcd  iro  dticribed  in  NAS  Report  1690,  "Repid  Uciva- 
tion:  Siyni ficence,  Needt,  Opportunitieti )96H,  Ihe  reioeroh  de* 
scribed  herein  is  pert  of  An  ongoing  progreni  in  disintegration  of 
hard  rock  for  the  rapid  excavation  of  tunnels ,  power  stations,  mines, 
and  underground  military  installations.  The  rate  of  productivity 
in  hard  rock  in  1968  was  increasing  only  about  3  to  3-1/1?  percent 
per  year,  which  Is  insufficient  to  keep  hard  rock  excavation  costs 
constant,  and  projected  rates  of  excavation  are  still  minimal  in 
1973  for  the  methods  now  in  use. 

The  NAS  Report  1690  gave  first  priority  to:  "Development  of  pro¬ 
cesses  and  equipment  for  hot  ing  tunnels  and  shafts  in  hard  abrasive 
rock  and  for  reducing  down  time  experienced  with  current  mining  and 
tunneling  machines  in  soft  or  medium  rock  masses".  It  was  further 
noted  that  of  the  several  novel  rock  disintegration  techniques  then 
being  studied  only  a  few  showed  promise  for  application  within  two 
decades.  There  is  a  good  Indication  that  the  method  described  in 
this  report  may  be  developed  for  practical  application  within  a 
period  of  five  years. 

Any  means  of  excavation  that  can  contribute  to  the  reduction  of 
costs  of  mining  minerals  or  to  increasing  their  production  rate 
directly  contributes  to  the  Improvement  of  the  standard  of  living, 
the  increase  of  mineral  reserves,  etc.  The  availability  of  underground 


t 


hl9h  trintportAtlon  uf  p«»op1o,  «.g,i  tn  the  Northeitt  Corridor 
(Weihington,  New  Vork.  (lu»ton),  with  the  incidoot  requirtmtnt  for  rtpid, 
Acreiitkhly  ecuiogicel  methtvU  of  OMCAvetion  of  hArti  rork,  ii  ah  e»x»r>* 
ciAi  fActur  in  th*  functioning  of  lociety  in  highly  popuiAted  Ar«AS.  Cx- 
CAVAtion  of  openingi  in  rock  undtr  popuUtod  Areoa  will  ripuire  rApid* 
quiet  methodi  of  rock  dikinUgrAtion. 

Achievement  of  the  requit'od  improveeientk  in  rapid  excAvation  of  hard 
rock  will  contribute  to  the  projected  improveiiient  of  effective  mask 
transportation  in  heavily  populated  areas,  will  ex.'end  our  mineral  re¬ 
sources  and  make  their  production  more  econoo'ical,  aod  will  improve  methods 
of  excavation  of  underground  civil  and  military  installations.  The  latter 
two  factors  directly  affect  the  whole  of  our  society. 

The  successful  development  of  a  new  method  of  fragmentation  will 
lead  to  the  design  of  the  basic  elements  of  a  new,  rapid  excavation 
system  for  hard  rock.  It  will  advance  the  knowledge  of  I'ock  disintegra¬ 
tion,  thermal  response  of  earth  materials,  energy  transfer  to  rock,  and 
the  relationship  of  those  factors  to  rock  and  mineral  properties.  This 
will  in  turn  lead  to  the  design  and  construction  of  a  prototype  system 
for  pilot  studies  and  the  later  development  of  full  scale  equiptiient  and 
excavation  systems. 

The  excavation  Industry  as  a  whole  Is  basically  Interested  in  the 
improvement  or  methods  of  mining  of  earth  materials  as  indicated  by  re¬ 
presentatives  of  government  and  industry  in  NAS  Report  1690,  at  the 
OECD  Tunneling  Conference  (Washington  1971)  and  the  Conference  on  Rapid 
Excavation  and  Tunneling  (Chicago  1972),  by  the  Interagency  Committee 
on  Rock  Mechanics,  and  the  committees  sponsored  by  the  National  Academy 
of  Science. 
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ThttmA)  (fir«  rotting)  hAvt  h9«n  (rniploycid  $1nc«  Antiquity 

In  th»  frigmtntdtlon  of  rock.  The*  i«ott  active  use  of  boat  In  recent 
and  current  operations  has  been  that  of  the  flame  Jet  for  drilling  hard 
rock  and  cutting  and  facing  of  spallable  monu«nent  stone.  Its  successful 
employment  led  to  the  investigation  of  the  feasibility  of  a  flame  Jot 
tunneling  machlnei  I.e.,  the  utllliatlon  of  flame  Jets  for  weitkenlng 
rock  to  make  It  more  susceptible  to  roller  cone  cutting.  Slot  cutting 
by  flame  Jets  and  wedging  of  ridges  between  kerfs  were  also  suggested 
as  a  possible  excavation  process.  Large  amounts  of  excess  heat  and 
combustion  generated  gas,  Inefficiency  of  heat  transfer  from  a  flame 
Jet  to  rock  and  other  factors  appeared  to  mitigate  against  its  use. 
as  well  as  the  resulting  limited  Increase  of  rock  removal  efficiency 
by  roller  cones. 

The  advantages  of  the  use  of  tunnel  boring  machines  in  soft  to 
medium  hard  rock,  which  they  can  excavate  economically,  and  the  disad¬ 
vantages  of  drill  and  blast  methods  (cyclic  character,  overbreak  Into 
walls,  and  vibration  and  noise  problems)  In  hard  rock  have  led  to  an 
Intenjive  search  for  novel  method(s)  which  will  Increase  tunneling 
capal'11  ities  primarily  In  hard  rock.  Of  the  several  methods  recently 
Investigated,  the  use  of  heat  has  been  the  subject  of  considerable 
research,  some  fourteen  techniques  having  been  proposed  and  investigated 
for  drilling  and  breaking  rock.  Flame  jets  have  been  in  use  for  a  long 
time,  but  are  limited  to  application  In  drilling  and  quarrying.  Various 
types  of  heat  drills  have  not  been  applied  to  date  because  of  non¬ 
adaptability  of  equipment  to  underground  conditions,  environmental  prob¬ 
lems  created  by  the  excavation  device,  high  energy  consumption,  or  be¬ 
cause  results  obtained  in  the  laboratory  could  not  be  applied  to  rock 
in  place. 
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It  shouVi  be  noted  thet»  while  recent  technical  literature  has 
etnphasl/ed  the  specific  energy  of  drilling,  breakage  and  other  types 
of  rock  removal,  the  absolute  value  of  this  factor  for  a  given  rock 
and  process  combination  Is  seldom  a  deciding  factor  by  Itself  unless 
a  critical  value  for  specific  energy  Is  exceeded.  The  most  critical 
operational  factor,  assuming  specific  energy  requirements  are  not  near 
or  above  maximum  permissible  values,  Is  the  rate  of  advance,  for  both 
drilling  and  excavation. 

While  some  novel  methods  of  rock  drilling  and  excavation  have 
been  found  to  be  technically  feasible,  equipment,  energy  or  other 
operational  costs  prohibit  their  use.  Hence,  this  investigation  Is 
designed  to  evaluate  both  the  technical  and  the  economic  feasibility 
of  the  application  of  the  method  of  excavation  by  use  of  thermal  rounds. 

Future  research  Is  designed  to  evaluate  further  the  factors  of  Im¬ 
portance  In  the  development  of  a  thermal  round  (of  the  same  concept 
as  a  high  explosive  blasting  round)  for  fragmentation  of  hard  rock 
In  such  a  manner  that  the  process  can  be  employed  In  operating  ex¬ 
cavation  systems.  This  includes  a  study  of  heat  sources,  hole  place¬ 
ment  and  geometry,  pertinent  mineral  and  rock  properties,  heat  trans¬ 
fer,  temperature  distribution,  thermal  stress  distribution,  fracture 
initiation  and  propagation,  required  strain  relief,  methods  of  slot 
cutting,  and  related  factors.  Based  upon  these  findings  appropriate 
existing  equipment  will  be  modified  or  new  equipment  designed  to  ac¬ 
complish  the  operation  of  the  fragmentation  process.  Three  years  of 
research  (as  of  June  1973)  have  been  completed  on  this  problem. 
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A  reasonable  assessment  of  a  method  of  thermal  fragmentation  and 
excavation  can  be  made  by  Investigation  of  the  following  factors: 

(1)  Current  status  of:  a)  theory,  b)  research,  and  c)  development. 

(2)  Energy  and  fracture  characteristics:  a)  specific  energy  require¬ 
ments,  b)  efficiency  of  energy  transfer,  c)  availability  of  energy,  and 
dj  control  of  energy  source. 

(3)  Application  characteristics:  a)  control  of  energy  within  the  rock, 

b)  preservation  or  strengthening  walls  and  roof,  c)  smoothness  of  opening 
excavated,  d)  efficiency  of  energy  utilization,  e)  environmental  factors, 
f)  adaptability  to  variable  rock  properties  and  conditions,  g)  feasibility 
of  equipment  design  and  construction,  h)  continuity  of  method  -  down  time, 
1)  research  and  development  required,  j)  time  to  develop  and  apply  new 
method,  and  k)  rate  of  advance. 

(4)  Cost  factors;  a)  research  and  development,  b)  capital  equipment, 

c)  operating:  labor,  materials,  etc.,  and  d)  overhead. 

The  above  factors  have  been  studied  in  relation  to  the  research  pro¬ 
ject  under  this  project  (Contract  No.  USDI-M0220068  [ARPA]),  Theoretical 
and  experimental  studies  are  in  advanced  stages  on  heat  transfer,  thermal 
stress  analysis,  and  thermal  properties  of  minerals  and  rocks.  Field 
tests  of  thermal  fragmentation  of  granite  have  reached  a  stage  which 
indicates  a  high  probability  of  successful  application.  Rate  of  ad¬ 
vance  and  cost  analysis  studies  are  in  progress. 
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Scientific  and  Engineering  Background 

While  there  have  been  specialized  studies  made  of  some  of  the 
thermal  properties  of  rocks,  mostly  in  relation  to  flame  jet  drilling, 
only  a  limited  number  of  investigations  have  been  made  of  the  specific 
scientific  and  engineering  factors  which  are  directly  related  to  the 
problems  posed  by  this  investigation.  There  is  a  large  fund  of  in¬ 
formation  from  ceramics  studies  of  heating  and  melting  of  synthetic 
silicates  and  selected  clays,  studies  of  the  calcination  of  carbonates 
in  their  sintering  for  manufacture  of  Portland  cement,  and  the  physical 
chemistry  of  some  types  of  rock  in  geologic  studies  of  igneous  action 
and  metamorphic  process*  ,  There  is  little  information  on  the  thermal, 
physical,  or  electrical  properties  of  molten  rock. 

Thermal  Rock  Drilling  and  Fracture 

Twelve  methods  of  drilling  or  breaking  rock  by  thermal  means  were 
reviewed  by  Maurer  in  1968  (1).  These  included  flame  jet  (jet  piercing), 
electric  disintegration  by  low  frequency  current,  cyclic  heating  and 
cooling  with  superheated  steam  and  liquid  nitrogen,  high  frequency 
currents,  microwaves,  induction,  fusion  and  vaporization  (LASL),  nuclear 
heat  (penetration),  electric  arc  (drill),  plasma  drill,  electron  beam, 
and  lasers.  Carstens  (2)  reviewed  several  of  these  methods  again  in 
1972,  and  described  an  additional  method  of  forming  an  internal  thermal 
inclusion  and  fracturing  rock  (3).  Considerable  research  has  been  ac¬ 
complished  in  the  past  four  years,  particularly  on  electron  beams  (4), 
lasers  (5  and  6),  with  flame  jets  (7,  8  and  9).  high  temperature  pene- 
trators  (lO),  plasma  jets  (ll),  resistance  wire  and  electric  arc 
heaters  (12). 
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Several  other  approaches  have  been  suggested  in  the  last  decade 
for  the  possible  use  of  heat  application  within  the  rock  to  cause  frag¬ 
mentation.  The  fracture  of  rock  has  been  the  subject  of  several  patents 
by  Sarapu  (13)  starting  in  1965.  In  a  study  of  laser  effects  on  rock 
Zar  (5)  proposed  that  a  heated  annulus  of  the  surface  of  a  hard  rock 
face  would  cause  the  rupture  of  a  section  of  the  rock  face.  However, 
this  assumes  that  the  rock  is  free  from  fractures,  joints,  etc. 

Thirumalai  (14)  was  the  first  to  report  (1970)  on  a  method  of  form¬ 
ation  of  an  internal  therr/ial  inclusion  by  dielectric  heating.  Labora¬ 
tory  test  blocks  of  granite  and  basalt,  somewhat  larger  than  1  foot 
cubes  were  fragmented  without  melting  by  localized  heating  below  600°C, 
the  heated  volume  being  less  than  2  percent  of  the  total  rock  volume. 

As  a  first  approximation,  a  thermal  inclusion  was  considered  to  exert 
a  pressure  on  the  inside  of  a  cylindrical  hole  in  an  infinitely  thick 
cylinder.  However,  quartzite  blocks  could  not  be  fragmented  by  dielectric 
heating.  The  total  electrical  energy  required  to  fracture  Dresser  basalt 
blocks  varied  from  0.067  to  0.100  kwhr,  and  the  volume  of  the  thermal  in¬ 
clusion  from  105  to  206  cc. 

Coiled  wire  resistance  heaters  of  Kanthal  wire  at  1000°C  were  found, 
during  the  investigation  described  herein,  to  fracture  hard  granite  when 
they  were  placed  in  pneumatically  drilled  holes,  but  the  heaters  failed 
after  short  term  usage.  Electric  arcs  between  carbon  electrodes  generate 
an  effective  thermal  inclusion  in  multiple  holes  in  solid  granite  and 
create  thermal  stress  fields  which  will  cause  fractures  over  distances 
of  2  feet  or  more  at  reasonable  energy  levels.  Tests  to  date  have  yielded 
promising  results  for  possible  application. 


For  granite  the  electrical  energy  was  0.10  kwhr  and  the  volume 
of  the  thermal  inclusion  was  206  to  350  cc.  Time  for  fracture  was 
less  than  2  minutes.  In  tests  2  foot  cubes  of  Missouri  red  granite 
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were  easily  fractured  by  means  of  a  rock  (partially  fused)  melt  created 
in  a  drill  hole  in  the  center  of  the  rock  with  carbon  electrodes.  The 
total  electrical  energy  was  approximately  1  kwhr  and  the  volume  of  the 
melt  about  380  cc,  the  total  heated  volume  of  rock  being  larger.  This 
represents  about  0.2  percent  of  the  total  volume  of  the  block.  Similar 
results  were  obtained  in  granite  in  situ. 

The  amounts  of  electrical  energy  per  unit  volume  of  rock  are  in 
the  ratio  of  1.6:1  for  Missouri  red  granite  as  compared  to  the  basalt 
and  granite  tested  by  Thirumalai.  The  time  for  fracture  and  the  total 
energy  applied  vary  approximately  as  the  volume  of  the  block,  indi¬ 
cating  that  the  cube  root  law  for  energy  usage  may  apply  for  similar 
rocks  and  geof.ietries. 

In  each  of  the  above  tests,  as  in  investigations  reported  by 
Nixon  (4)  with  electron  beams,  radial  cracks  as  well  as  cracks  normal 
to  the  axis  of  the  hole  or  inclusion  were  formed.  Nixon  (4)  reports 
that  temperature  distributions  were  calculated  by  empirical  formulas, 
with  the  assumption  of  a  1700®C  melting  temperature  at  the  boundary  of 
a  finite  element  prov.'am  utilizing  values  of  Young's  modulus,  Poisson's 
ratio  and  thermal  coefficient  of  expansion  as  functions  of  temperature. 
Both  tangential  and  radial  stresses  near  the  cavity  were  found  to  be 
compressive,  while  in  the  cold  zone  the  stresses  were  tensile. 

Nixon  also  came  to  the  conclusion  that  for  electron  beam  penetra¬ 
tion  and  heating  the  compressive  stresses  and  gradients  are  typically 
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larger  by  an  order  of  magnitude  than  the  tensile  stresses  and  gradients. 
The  expansion  of  the  cavity  along  its  axis  and  the  radial  growth  of 
the  heated  region  cause  tensile  stresses  in  the  axis-radius  plane  and 
in  the  tangential  direction.  At  depths  equal  to  the  melt  cavity  and 
below,  the  tensile  stresses  in  the  r-axis  plane  are  larger  than  the 
tangential  stresses  creating  favorable  conditions  for  a  crack  normcTl 
to  the  axis.  The  calculations  were  confirmed  only  by  the  behavior  of 
blocks  of  quartzite,  however,  and  not  with  other  types  of  rock. 

Thirumalai  proposed  that  dielectric  internal  heating  to  form  a 
thermal  inclusion  be  applied  to  ore  and  rock  in  stopes  where  two  free 
faces  exist  for  displacement.  However,  the  feasibility  of  breakage 
by  this  means  was  not  tested.  He  (15)  also  suggested  two  types  of 
tunnel  rounds  using  thermal  inclusions,  one  with  a  cut  at  the  bottom 
of  the  face,  and  a  V-cut  type  round  with  no  strain  relief  other  than 
the  tunnel  face.  The  feasibility  of  such  tunnel  rounds  also  was  not 
investigated. 

Most  of  the  novel  methods  which  have  been  proposed  for  the  rapid 
excavation  of  hard  rock  are  being  considered  in  terms  of  several  years 
for  possible  development  and  application.  Some  have  inherent  disad¬ 
vantages  which  may  not  be  easy  to  overcome.  For  example,  the  x-rays 
created  by  electron  beams  are  difficult  to  shield.  Current  laser  equip¬ 
ment  is  not  adaptable  for  use  underground.  Plasma  Jets  are  inefficient 
for  breaking  rock. 

On  the  other  hand,  many  of  the  principles  of  the  thermal  fragmenta¬ 
tion  method  proposed  herein  have  been  Investigated  in  sufficient  detail 
to  indicate  a  strong  probability  of  successful  application  within  a 
relatively  short  period  of  time. 
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The  principles  of  breakage  employed  are  similar  to  those  which 
have  been  successfully  applied  in  blasting.  A  large  part  of  the  equip¬ 
ment  to  be  used  in  further  research  and  development  is  available  for 
purchase,  i.e.,  air  compressors,  drills,  power  sources,  carbon  elec¬ 
trodes,  etc.  Some  adaptation  of  equipment  will  be  required,  and  only 
a  few  of  the  components  will  require  new  design  and  development.  All 
of  the  components  envisioned  are  adaptable  to  use  in  underground  en¬ 
vironments.  Environmental  problems  will  be  minimal,  i.e.,  less  venti¬ 
lation  will  be  required  than  with  explosives.  Some  of  the  problem 
areas  are  to  be  found  in  the  development  of  rugged  carbon  electrodes, 
working  in  very  wet  tunnels,  investigation  of  other  high  temperature 
heaters,  etc. 

Several  methods  utilizing  rock  melting  combined  with  fracture  pro¬ 
cesses  have  been  proposed  by  Los  Alamos  Scientific  Laboratory  (10,  16 
and  17)  for  the  excavation  of  tunnels.  These  include  (a)  full  face 
melting  (Subterrene  2),  (b)  partial  face  melting  and  solid  core  excava¬ 
tion  (Subterrene  1),  and  (c)  a  multiple  hot  prong  melting  and  fractur¬ 
ing  head.  Experimentation  has  been  performed  on  penetration  of  rocks 
and  soils  with  2  in.  models  which  form  a  glass  lined  hole  in  porous 
rock  and  soil,  and  eject  melted  material  from  the  hole  in  non-porous 
rocks. 

Maurer  (1)  notes  that  in  early  tests  (17)  a  5  kw  heater  drilled 
a  hole  5.1  cm  diam  in  basalt  at  a  rate  of  1.2  cm/min,  which  he  further 
reports  corresponds  to  a  specific  energy  of  12,300  Joules/cm^  How¬ 
ever,  the  basalt  used  at  Los  Alamos  is  porous,  some  of  it  having  a 
density  of  1.43  whose  heat  of  fusion  would  be  2720  Joules/cc  of  hole 
created.  This  ratio  of  heat  required  of  over  4  to  1  should  be  expected 
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because  some  of  the  rock  around  the  hole  is  melted.  If  the  work  re¬ 
quired  to  force  the  drill  into  the  rock  is  taken  into  account,  the 
ratio  would  be  greater. 

Such  a  device  has  two  disadvantages.  The  minimum  energy  required 
to  melt  holes  in  granite,  assuming  all  of  the  energy  goes  into  the 
melted  rock,  is  in  the  neighborhood  of  5000  joules/cc,  while  that  for 
a  pneumatic  drill  is  250  to  400  joules/cc.  The  drilling  (melting)  rate 
for  electric  heaters  is  (1)  1  to  3  cm/sec  whHe  that  for  pneumatic 
drills  is  50  to  75  cm/sec. 

The  feasibility  of  molten  rock  penetration  into  the  rock  mass  around 
a  tunnel  has  been  shown  for  porous  earth  materials  and  a  granite.  The 
stability  of  unreinforced,  brittle  ceramic  glass  tunnel  liners  has  not 
been  investigated,  nor  the  effects  of  heat  upon  surrounding  rock  and  soil. 

The  multiple  hot  prong  drilling  heater  suggested  by  LASL  (16)  appears 
to  have  two  disadvantages.  The  first  is  the  amount  of  energy  required 
to  make  holes  in  the  rock  face  by  melting.  The  second  concerns  the 
mode  of  fracture  from  internal  thermal  inclusions.  For  effective  frac¬ 
turing  lateral  displacement  relief  is  required,  which  cannot  be  pro¬ 
vided  by  the  drilling  head.  In  principle,  a  "no  cut"  round  is  proposed, 
which  has  been  found  to  be  feasible  for  explosives  in  some  types  of  rock. 
However,  experimentation  has  shown  that  more  displacement  relief  is 
required  for  crack  propagation  caused  by  thermal  stress  than  for  stress 
waves  produced  by  explosives. 

Because  of  the  high  heat  of  fusion  of  rock,  the  energy  requirements 
and  costs  for  melting  of  any  significant  portion  of  the  rock  are  criti¬ 
cal  factors.  For  example,  for  the  7  meter  diameter  tunnel  excavation 
proposed  by  LASL  (10)  50  Mw  of  power  are  required  for  melting  (of  one- 
fourth  of  the  face)  alone,  but  heat  absorption  and  losses  may  increase 
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this  by  as  much  as  three  limes.  Thus,  at  3  cents  per  kwhr  the  power 
cost  would  be  from  $750  to  $2250  per  foot  of  tunnel.  Ventilation 
costs,  also,  would  be  very  high. 

The  required  heat  flux  of  500  w/cm^  to  accomplish  necessary  rates 
of  advance  are  well  beyond  the  state-of-the-art,  e.n  ,30  w/cm^  for 
lithium  heat  pipes  and  200  w/cm^  for  electrical  heaters. 

Recent  investigations  (18,  19)  have  also  shown  that  complex  re¬ 
actions  take  place  in  radiant  heaters  with  carbon  heating  elements, 
molybdenum  shells,  and  boron  nitride  insulation.  The  average  life¬ 
time  of  such  a  penetrator  is  about  six  hours,  while  those  without  boron 
nitride  have  been  tested  to  approximately  45  hours. 
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Chapter  H 

ROCK  CHARACTERISTI C S 

The  response  of  rocks  to  the  application  of  heat  depends  upon 
their  composition,  structure,  density,  thermal  properties,  and  mechani¬ 
cal  properties.  The  percentage  of  SlOy  In  a  rock  affects  the  thermal 
behavior  markedly  as  well  as  the  strength,  brittleness,  and  other  fac¬ 
tors  . 

While  the  r'esearch  program  described  herein  Is  concerned  primarily 
with  hard  Igneous  and  motamorphic  rocks,  the  applicability  of  thermal 
fragmentation  to  sedimentary  rocks  Is  also  of  Interest. 

Igneous  Rocks 

Johannsen  (20)  proposed  a  classification  for  Igneous  rocks  which 
will  be  particularly  useful  in  the  present  study.  It  Is  based  on 
double  tetrahedron  (F1g.  l)  with  quartz  at  the  apex  at  the  top  and  K, 

Na  and  Ca  feldspars  at  the  common  base  corners  and  feldspatholds  at 
the  bottom  apex.  Thus,  rocks  containing  free  quartz  (or  without)  may 
be  classified  by  their  cojnpositlon  and  their  position  in  the  double 
tetrahedron  according  to  their  composition  and  class  name. 

In  a  consideration  of  the  physical  chemistry  of  rock  meUs,  the 
phase  rule  of  mineralogy  may  be  useful  a^  a  guide.  This  rule  states 
that  (21)  the  number  of  different  minerals  In  a  rock  should  not  ex¬ 
ceed  but  may  be  equal  to  the  number  of  components  In  the  magma  from 
which  the  rock  Is  formed.  It  Is  noted  that  over  90  percent  of  the 
earth's  crust  consists  of  silicate  minerals,  Including  quartz.  Hence, 
the  behavior  of  these  minerals,  which  has  been  studied  extensively  by 
ceramists  and  geologists.  Is  of  basic  interest  in  this  Investigation. 
Eltel  (22)  has  compiled  a  useful  comprehensive  treatment  of  the  thermal 
behavior  of  silicate  minerals  (to  1965).  However,  most  of  the  silicate 
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lyitttmi  ilu<lle(l  do  not  cornttpond  to  th«  compoiUluiu  of  Igntout  nocki. 
I^rlnclpil  vonttituontt  my  io  In  tho  «om«  propontlont,  but  imll  omountt 
of  occmory  wintriU  or  compoundi  wiiy  groiily  «ffoct  thn  nf'oportioi 
of  rocHi)  «1ntfrii  moUii  ond  glititK  produt;»d, 

Tho  princloil  typoi  of  pUitonk  ignooui  rocki  v«ry  In  ovorogo  corn* 
potUlon  from  26  ptrctnt  to  no  guortii  ond  moy  contiln  foldtpart  or  othtr 
111 ir.ato  mlnoroh  In  varying  amount!  (Tabla  1).  Tha  vulcanK  con\po!l~ 
tlonal  oqulvalant!  of  granltti  Includa  thu  rhyolltei,  dacUo!  and  por> 
phyrlei.  That  lii  typoi  of  plutonic  rocks  may  hava  approKlinato  volcanic 
•quivalanti  In  chaa^lca)  coi»po$1t1on.  Volcanic^  may  ba  hard  and  dense i 
or  soft  and  porous,  a.g.,  tuffs  and  pyruclasts. 

The  composition,  density  and  related  physical  properties  determine 
melting  temperature  and  general  response  to  application  of  heat.  Molt¬ 
ing  tetnperatures  of  noneutectic  compositions  are  higher  than  those  at 
or  near  the  eutectic  condition,  and  other  thermal  properties  are  closely 
related  to  rock  type,  texture  and  origin. 

Sedimentary  Rocks 

Sedimentary  rocks  are  easily  fractured  by  heat  in  small  test  blocks 
and  it  may  be  desirable  to  apply  the  proposed  method  to  them  in  the 
future.  They  vary  in  composition  more  than  Igneous  types.  An  average 
chemical  composition  (Table  2)  shows  that  they  contain  many  of  the  same 
combinations  of  elements  as  igneous  rocks.  However,  the  minerals  and 
rock  structure  are  the  result  of  large  scale  chemical  and  mechanical 
fractionation  p»'ocesses.  Their  thermal  and  mechanical  properties  ^  > 
different  from  igneous  rocks,  and  tunnels  can  be  excavated  in  most  of 
them  with  boring  machines.  The  most  common  sedimentary  rocks  are  sand¬ 
stone,  limestone,  and  shale. 


table  1 

ApproxiHate  Mineral  C<*position  of  Prtncipal  Types 
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TABLE  2 

Avtraige  Chemical  Composition  cf  Rocks  (21) 

Constit-  Igneous  Shale  Sandstone  Limestone  Sediment^  Sediment^ 
uent 


S10, 

T10* 

A1,0, 

FoaOs 

FeO 

MgO 

CaO 

NaaO 

KjO 

HaO 

P2O5 

CO2 

SOa 

BaO 

C 

MnO 


59.14 

1.05 

15.34 

3.08 

3.80 

3.49 

5.08 

3.84 

3.13 

1.15 

0.30 

0.10 

006 


58,10 

0.65 

15.40 

4.02 

2.45 

2.44 

3.11 

1.30 

3.24 

5.00 

0.17 

2.63 

0.64 

0.05 

0.80 


78.33 

0,25 

4.77 

1.07 

0.30 

1.16 

5.50 

0.45 

1.31 

1.63 

0.08 

5.03 

0.07 

0.05 


5.19 

0.06 

0.81 

0.54 

7.89 

42.57 

0.05 

0.33 

0.77 

0.04 

41.54 

0.05 


57.95 

0.57 

13.39 

3.47 

2.08 

2.65 

5.89 

1.13 

2.86 

3.23 

0.13 

5.38 

0.54 

0.66 


44.5 

0.6 

10.9 

4.0 

0.9 

2.6 

19.7 

1.1 

1.9 

0.1 

13.4 


0.3 


Total....  99.56  100.00  100.00  99.84  99.93  100.0 

^Shale  82,  sandstone  12,  limestone  6;  after  Leith  and  Mead. 

^Poldervaart,  Geol.  Soc.  Amer,  Spec,  Paper  62,  132,  1955. 
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Sandstone  may  be  composed  largely  of  particles  or  grains  of  quartz, 
of  quartz  and  feldspar,  or  they  may  be  graywackes  (basalt,  slate,  fel- 
site,  rhyolite,  schist,  quartz  and  feldspaths),  or  tuffaceous  sandstones. 
Shales  are  clastic  rocks  containing  silt  and  clay.  The  carbonate  rocks 
(general  term  limestone)  include  dolomites  and  limestone,  whose  primary 
constituents  are  calcium  carbonate,  calcium  magnesium  carbonate,  and 
other  carbonates.  They  may  also  contain  fragments  of  many  other  types 
of  minerals.  Perhaps  the  most  important  reaction  of  these  rocks  to  appli¬ 
cation  of  heat  is  their  loss  of  water  and  COj. 

Metamorphic.  Rocks 

The  mineral  content  of  metamorphic  rocks  is  somewhat  similar  to 
that  of  igneous  rocks.  They  include  schists,  gneisses,  quartzite, 
slate,  marble,  hornfels  and  many  others.  Many  of  them  are  hard,  strong, 
and  difficult  to  excavate.  The  conceptual  tunnel  system  for  the  NE  Cor¬ 
ridor  plans  for  considerable  excavation  through  metamorphic  rocks.  Their 
response  to  the  application  of  heat  is  generally  similar  to  that  of  ig¬ 
neous  rocks,  with  notable  exceptions.  One  is  a  quartzite  which  may  be 
90+  percent  quartz.  Sioux  quartzite,  for  example,  has  been  found  to 
be  difficult  to  melt  in  tests  at  UMR,  but  fractures  readily  when  heat 
is  applied  because  of  the  high  coefficient  of  thermal  expansion  of 
quartz  and  its  crystallographic  phase  change  at  STS^C.  Other  meta¬ 
morphic  rocks  have  their  own  individual  characteristic  response  to 
high  temperatures  depending  upon  composition,  grain  size,  structure. 


etc. 
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Thermal  Properties  of  Molten  Rocks  and  Glasses* 

Solidus  and  Hquidus  temperatures  of  representative  rock  types. 
Sufficient  information  about  chemical  compositions  of  rocks  and  phase 
equilibria  in  rock-forming  systems  is  available  to  provide,  on  comple¬ 
tion  of  a  thorough  synthesis,  a  reliable  estimate  of  the  temperatures 
of  initial  liquid  formation  and  of  complete  melting.  Ternary  systems 
involving  SiOj,  A«,203  ,  Na20,  K2O,  CaO,  MgO,  "FeO",  and  Fe203  have  been 
rather  well  investigated  and  quaternary  relationships  of  importance  are 
understood.  Reasonable  extrapolation  to  higher  order  systems  should  be 
possible,  though  some  direct  experimental  work  will  be  necessary  to 
fill  the  gaps. 

Nonequilibrium  assemblages.  Because  the  process  of  rapid,  localized 
melting  of  rocks  involves  a  considerable  thermal  gradient,  the  phase 
assemblages  will  necessarily  vary  from  equilibrium  or  near-equilibrium, 
through  all  stages  of  thermal  equilibria,  outward  to  unaltered  rock. 
Nonequilibrium  assemblages  formed  during  the  rapid  heating  and  cooling 
processes  are  expected,  however,  and  will  be  sensitive  to  secondary  rock 
properties  such  as  particle  size,  porosity,  presence  of  accessory 
nucleating  agents,  and  the  nature  of  the  transitory  environment  of 
melting,  particularly  water  and  carbon  dioxide  pressures.  Though  a 
good  general  picture  of  equilibrium  phase  assemblages  can  be  obtained 
from  available  phase  equilibria  data,  direct  investigation  of  real  rock 
will  be  necessary  to  determine  the  importance  of  nonequilibrium  assem¬ 
blages* 


*  Section  written  by  Dr.  Charles  A.  Sorrell,  Dept,  of  Ceramic  Engineering, 
University  of  Missouri-Rolla 
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Glass  transition  temperatures.  Experimental  determinations  of 
glass  transition  temperatures  have  been  largely  restricted  to  commer¬ 
cial  glass  compositions,  none  of  which  approximate  major  rock  types. 
Because  the  rate  of  melting  and  heat  transfer  will  be  affected  by 
the  manner  in  which  a  glass  inclusion  is  formed,  it  may  be  necessary 
to  obtain  glass  transition  temperatures  of  rock  compositions.  Of 
equal  importance  are  temperature-viscosity  relationships  in  liquid- 
crystal  assemblages  likely  to  be  encountered  because  of  incomplete 
melting  or  partial  crystallization  during  cooling  of  the  inclusion. 

Devitrification.  Strengths  of  glasses  are  strongly  dependent  not 
only  on  crystals  present  because  of  partial  melting  but  also  on  very 
small  crystals  formed  by  devitrification  below  glass  transition  tempera 
tures.  In  many  cases,  as  in  commercial  glass -ceramics  (Pyroceram), 
considerable  devitrification,  induced  by  nucleation  of  primary  phases 
on  precipitates  of  a  minor  component  (Ti02),  strengthens  the  material 
appreciably  and  improves  thermal  shock  resistance.  The  effects  of  de¬ 
vitrification  on  strength  and  other  properties  of  glasses  should  be 
studied  in  real  rocks  and  some  consideration  given  to  addition  of 
nucleating  agents  during  the  melting  operation  in  order  to  promote 
beneficial  devitrification.  This  is  important  in  relation  to  the  LASL 
studies. 

Volume  changes  associated  with  melting.  Though  some  data 

are  available,  the  bulk  volume  changes  which  occur  on  melting  are  com¬ 
plicated  by  variable  porosity  and  by  loss  of  volatile  components.  Bulk 
densities  of  rocks  and  densities  of  glasses  formed  from  them  should 
be  measured.  Johannsen  (20)  has  indicated  that  glasses  occupy  from 
5  to  20  percent  greater  volume  than  their  crystalline  rocks. 
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Alteration  of  rock  composition  by  additives.  Nonsilicate  rocks, 
such  as  limestone,  and  some  silicate  rocks,  notably  those  rich  in 
Si02,  will  not  form  glasses  or  will  do  so  only  at  very  high  tempera¬ 
tures.  In  general,  addition  of  other  oxide  components  will  lower 
eutectic  and  liquidi'o  temperatures  considerably.  Deliberate  intro¬ 
duction  of  Na20  and  Si02  to  the  area  of  melting  may  be  expected  to  re¬ 
sult  in  glass  formation  in  a  limestone  in  the  500-700°C  range. 

Transitory  Environments 

Oxidation-reduction  reactions.  Some  elements,  notably  iron,  are 
readily  reduced  or  oxidi2ed,  with  associated  structural  changes,  dis¬ 
sociation  and  recrystallization,  wnich  affect  the  physical  strength 
of  the  rock.  Though  the  iron  oxides  are  well  understood,  the  oxidation- 
reduction  relationships  in  iron-bearing  silicates  have  not  been  studied 
thoroughly.  Provision  of  a  transitory  environment  during  melting,  more¬ 
over,  complicates  the  situation  even  more,  with  the  nature  of  the  environ¬ 
ment  In  large  part  determined  by  the  composition  of  the  heating  element. 

All  experiments  with  materials  likely  to  undergo  oxidation  or  reduction 
should  be  done  in  at  least  two  ways,  one  with  a  locally  reducing  en¬ 
vironment,  as  in  graphite  crucible,  and  the  other  in  an  oxidizing  en¬ 
vironment,  as  in  a  platinum  crucible. 

Effects  of  volatile  components.  Melting  and  glass  formation  in 
some  materials  will  be  preceded  or  accompanied  by  formation  of  volatile 
components  which,  in  addition  to  changing  the  amount  of  condensed  material, 
will  provide  a  transitory  environment.  Unfortunately,  only  meager  in¬ 
formation  regarding  phase  equilibria  in  complex  systems  containing 
water  and  CO2  is  available  and  even  less  information  on  partial  pres¬ 
sures  or  vaporization  rates  of  volatile  components  of  the  silicates. 


22 


Attempts  to  predict  the  magnitudes  of  the  effects  of  volatile  com¬ 
ponents  on  melting  and  glass  formation  must  be  secondary  to  experi¬ 
mental  measurements. 

Heating  element-melt  compa/cioility.  In  general,  silicate  melts 
do  not  react  with  metals  or  carbon  unless  the  metals  are  pemitted 
to  oxidize.  Use  of  metallic  heating  elements  in  contact  with  melts 
should  cause  no  problem  unless  the  stability  of  the  element  is  de¬ 
pendent  on  formation  of  an  oxide  layer.  The  same  generalizations  are 
appropriate  to  the  refractory  carbides.  Experimental  evaluations  of 
different  heating  element  materials  in  contact  with  different  melt 
compositions  should  be  made. 

Experimental  Studies  of  Melts  and  Glasses 

A  reasonably  thorough  understanding  of  the  chemical  processes 
which  occur  during  rapid,  localized  melting  of  major  rock  types,  with 
or  without  additional  components,  should  be  obtainable  through  a  de¬ 
tailed  synthesis  of  available  physical  chemical  data.  A  limited, 
carefully  selected  set  of  experiments  is  planned  in  the  future, 
utilizing  natural  rock  and  soil  samples,  to  determine  the  following: 

1.  temperatures  of  initial  and  complete  liquid  formation  in  re¬ 
ducing  and  oxidizing  environments, 

2.  phase  assemblages  developed  as  functions  of  temperature  and 
distance  from  the  heating  element  -  rock  interface, 

3.  volume  changes  associated  with  the  melting  process, 

4.  compatibility  of  heating  element  materials  with  silicate  melts, 

5.  glass  transition  temperatures  of  common  rock  compositions, 
viscosity,  dissociation,  electrical  properties, 

6.  physical  changes  in  unmelted  or  partially  melted  rocks  re¬ 
sulting  from  decrepitation,  volume  changes,  thermal  fracturing. 
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Closely  related  to  these  are  studies  of  the  electrical  character¬ 
istics  of  natural  rock  melts  In  arc  created  environments. 

Heated  and  melted  rocks.  The  heat  capacities,  melting  points,  heats 
of  fusion,  and  thermal  coefficients  of  expansion  of  hard  rock  forming 
minerals  are  of  primary  Importance  In  the  present  study  and  some  have 
been  tahfilated  (23).  The  thermal  coefficient  of  expansion  of  the  sili¬ 
cate  minerals  varies  with  temperature,  crystallographic  direction  and 
composition  (24).  The  coefficient  for  quartz  (plus  its  crystallo¬ 
graphic  change  at  573®C)  Is  larger  than  that  for  most  silicates 
(Table  3)  which  Is  one  of  the  primary  causes  of  thermal  spallation  of 
rocks  containing  a  relatively  high  percentage  of  free  silica  (25).  How¬ 
ever,  experiments  performed  at  UMR  under  the  current  program  show  that 
rocks  composed  of  minerals  with  low  coefficients  of  expansion,  such  as 
basalts  and  dolomites,  will  fracture  readily  under  the  Influence  of 
internal  thermal  stress  fields. 

Data  on  the  melting  temperatures  of  silicate  systems  have  been 
acquired  In  studies  of  ceramic  glasses,  and  have  been  compiled  because 
of  geological  Interest  (26).  Although  the  melting  point  of  some  of 
the  constituents  In  the  A1I2O3— K2O— SIO2  system,  as  with  other  systems, 
are  quite  high,  eutectic  mixtures  may  have  lower  melting  points  (Fig.  2). 
Certain  combinations  of  glass  systems  are  possible  (27).  It  Is  reason¬ 
able  to  assume  that  systems  which  result  from  glass  formation  from 
some  kinds  of  rock  will  be  more  corrosive  than  others  on  heaters  to 
be  employed  for  heating  and  melting,  e.g,,  if  iron  sulphide  Is  present. 

Carbonaceous  rocks  which  are  heated  under  pressure  are  decarbonated 
(Table  4)  and  are  a  source  of  volatiles  which  may  cause  problems. 
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TABLE  3 


Coefficients  Of  Volume  Expansion,  _  l/dV|, 

Of  Selected  Silicates  At  Atmospheric  Pressure  (23) 


Compositions,  where  stated,  are  in  mole  per  cent.  The  coefficients  were 
derived  by  analytically  differentiating  polynomial  functions  fitted  to  the 
volume-temperature  data  by  the  method  of  least  squares.  Most  sets  of 
measurements  conmence  at  room  temperature  giving  poor  control  of  the  de¬ 
rived  functions  at  20°C,  compared  to  the  control  at  higher  temperatures, 
and  leading  to  an  uncertain  estimate  of  at  20°C.  The  values  enclosed 
in  parentheses  were  derived  by  extrapolating  from  temperatures  where  ade¬ 
quate  control  on  the  functions  could  be  established.  Less  credence  should 
be  given  the  values  at  20°C  than  those  at  400°C  or  800°C. 


Compound 

20°C 

av  X  10*°C'* 
400®C 

800°C 

Ref 

Akermanite 

(26) 

30 

33 

11 

Andalusite 

16 

29 

43 

19 

Coesite 

8 

11 

14 

18 

Cordierite  (synthesized  at 

(14) 

18 

22 

11 

1420OC) 

Feldspars 

Adularia  0re**3Ab9 ‘lAnj *4 

14 

20 

24 

3 

Microcline  Ores'sAbie’s 

(7) 

17 

23 

15 

Plagioclase  Ab^sAni 

18 

27 

33 

6 

Plagioclase  AbyrAnjs 

12 

19 

24 

6 

Plagioclase  Ab$tAn44 

13 

17 

20 

6 

Plagioclase  AbsAn«s 

12 

12 

20 

6 

Garnets 

Almandine 

15 

25 

30 

17 

Andradite 

21 

25 

29 

17 

Grossularite 

18 

23 

28 

17 

Pyrope 

19 

26 

30 

17 

Spessartite 

15 

28 

34 

17 

Gehlenite 

(23) 

25 

26 

11 

Hornblende 

23 

28 

33 

8 

Kyanite 

11 

28 

30 

19 

Merwinite 

(29) 

38 

42 

11 

Mullite 

9.5 

15 

17 

1 

Nephel ine 

Ne7»KS2i 

31) 

53 

72 

16 

NejiKSm 

(66) 

51 

49 

16 

25 


TABLE  3  -  continued 


Compound 

0 

20^0 

V  X  10®°C"* 
400OC 

800°C 

Ref. 

01  ivine 

Fai  00 

(27; 

30 

31 

12 

FAs  0  FO2  0 

(26! 

32 

34 

13 

FA41 FO59 

27| 

32 

35 

13 

Fdi sFOas 

25 

32 

39 

13 

Faio*iF089*9 

(23  j 

31 

39 

10 

FOioo 

24 

38 

44 

18 

Monticell ite 

(32] 

36 

39 

11 

Fe-Monticell ite 

(25j 

32 

38 

12 

Pseudo -wollastonite 

(30] 

32 

36 

11 

Pyroxenes 

Aug ite 

18 

25 

32 

9 

Clinoenstatite 

25 

29 

33 

11 

Oiopside 

24 

28 

32 

7 

Jadeite 

(20) 

29 

(38) 

20 

Quartz 

34 

69 

-  3 

2,4,14 

Slllimanite 

10 

18 

26 

19 

Topaz 

14 

20 

25 

5 

Zircon 

9.3 

13.7 

17.7 

1 

26 


TABLE  4 


Reactions  Involving  Decarbonation  (23) 


Temperatures  (°C) 

Pressure 

(kb) 

.25 

.5 

1.0 

2.0 

Ref. 

Calcite  -  lime  +  CO2 

CaCOa  *  CaO  +  CO2 

Magnesite  -  periclase  +  CO2 

MgCOa  «*  MgO  +  CO2 

650 

710 

775 

850 

4 

Rhodochrosite  -  manganosite  +  CO2 

MnCOs  “  MnO  +  CO2 

615 

660 

710 

775 

1 

Smithsonite  -  zincite  +  CO2 

ZnCOs  ■  ZnO  +  CO2 

•  * 

285 

340 

410 

3 

Dolomite  -  calcite  +  periclase  +  CO2 
CaMg(C03)2  “  CaCOa  +  MgO  +  CO2 

710 

765 

825 

900 

4 

Calcite  +  quartz  -  wollastonite  +  CO2 

CaCiOa  +  Si02  “  CaSiOs  +  CO2 

•  • 

630 

680 

745 

5 

Wollastonite  +  calcite  -  spurrite  +  CO2 
2CaS10$  +  SCaCOa  ■  Ca5Si208(C0i)  +  2CO2 

960 

*• 

•  • 

•  • 

8 

Calcite  +  diopside  +  forsterite 
-  monticellite  +  CO2 

ZCaCOa  +  CaMgSiaOg  +  Mg2S10i, 

=  3CaMgSi04  +  2  CO2 

780 

865 

965 

•  * 

6 

Calcite  +  diopside  -  akermanite  +  CO2 

CaCOa  +  CaMgSi206  “  Ca2MgSi207  +  C02 

845 

890 

965 

•  » 

6 

Calcite  +  wollastonite  -  Tilleyite  +  CO2* 
3CaCOj  +  2CaSi0a  "  Ca5Si207(C03)2  CO2 

890 

•  » 

•  ■ 

•  • 

2 

Tilleyite  -  spurrite  +  COa* 

CasSl207(C0j)2  ■  CasSia Oe (CO3 )  CO2 

935 

•  • 

t  • 

•  • 

2 

*  Tilleyite  could  be  synthesized  only  in  the  presence  of  small  amounts 
of  AUOj  and  fluorine,  although  the  latter  may  have  acted  merely  as  a 
flux. 


X«Kj0.Al,0,.4Si0, 

(leucite) 

W«Ka0.A1,0,,.6S10, 

(potash  feldspar) 

E'»Kj0.Ala0j.2S10, 

(hexagonal  an^ 
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tunotl  liortoQ  iMCh1n«ii,  iom«  mt,y  b»  of  (o  tbU  A)4o» 

wh<»n  biio*  t.«rboo  eUt;t»'oo»<i  Ar'e  tNnployeil  in  .1  tlrlU  hole  in  e  iHuk, 
tU)  end  CO;  miy  be  fomert  ei  the  elet:t»’(Hloi»  ere  lovolveil  In  cnewu\tl 
reletlons. 

Ai  Imllceted  ebove.  other  Importent  behevlorel  fectorik  of  »lllt:e 
«i  sllicete  nielts  ere  concerned  with  the  voUtllee  which  e»\>  prvhluced, 
end  corroilvlty  of  the  melt,  Witer  content  heii  e  merhed  effect  on 
melting  pointt,  end  when  mlnereli  heve  combined  wetar  their  (hydi'o* 
themel)  behavior  It  different  from  that  of  the  dry  illlcetes, 

Molten  tllicatet  ere  believed  to 

fall  In  the  category  of  being  Newtonian  fluldt,  whote  vUcotUy  It 
Independent  of  thear  rate.  Viscosity  Is  0  simple  activated  process 
which  may  be  represented  by  the  function  »>  •  iioexp(f  yHT)  (where  no 
Is  the  limiting  viscosity  at  high  temperatures  and  Is  the  activa¬ 
tion  energy),  This  law  Is  obeyed  by  silica,  and  more  complex  silicate 
liquids  may  be  represented  by  a  sum  of  terms. 

The  viscosities  of  SIOj,  obsidian,  granite  end  andesite,  ere  higher 
than  thosQ  of  less  siliceous  rochs  (Table  5),  while  those  of  basalts 
and  lava  are  quite  low  at  1400*C,  This  has  been  found  to  be  true  of 
the  rocks  melted  with  carbon  arcs  In  research  at  IWR.  It  was  difficult 
to  form  a  melt  from  nearly  pure  quartxite,  and  the  molt  from  granite 
containing  about  20  percent  fret  quartz  was  very  viscous,  Under  the 
same  tost  conditions  the  molt  forewd  from  a  basalt  flowed  easily, 

The  amount  of  SIOj  present  in  silicate  systent  (Table  [5])  has 

a  marked  effect  on  its  viscosity,  the  viscosity  Increasing  with  the 
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i  percentage  of  SiO?.  Thus,  more  siliceous  rocks  would  require  higher 
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ttmptraturis  and  more  heat  to  cause  them  to  flow,  or  artificially  in- 
troductd  fluxes  would  be  required  to  reduce  the  viscosity. 
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Chapter  III 
ELECTRIC  ARCS 
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Arcs  -  Rock  DrllUng 

A  review  of  the  important  characteristics  of  electric  arcs  was 
made  by  McMaster,  et  al  (28)  In  relation  to  an  Investigation  of  the 
feasibility  of  their  use  in  drilling  oil  wells.  Their  performance 
in  air  and  water  was  of  direct  interest  in  these  studies,  and  much 
of  this  information  is  of  basic  interest  in  the  rock  fragmentation  re¬ 
search  with  which  this  project  is  concerned. 

The  low  voltage  required  for  arcs  may  be  advantageous  for  rock 
breakage,  i.e.,  it  is  generally  easier  and  safer  to  use  in  underground 
environments.  McMaster  (28)  reported  that  the  voltage  drop  (dc)  at 
the  cathode  is  of  the  order  of  the  least  ionization  potential  of  the 
gas  or  vapor  in  which  it  burns.  The  arcs  employed  in  drilling  experi¬ 
ments  were  subject  to  severe  transients  and  rapid  fluctuations  in 
length,  position,  current  and  voltage  so  that  the  average  character¬ 
istics  might  not  correspond  to  a  steady  state  operating  condition. 

The  temperature  of  the  arc  column  in  air  was  reported  (28)  to  vary 
from  4330°C  to  at  atmospheric  pressure,  and  as  high  as  10,760°C 

in  the  center  of  the  arc  column  for  high  current  arcs.  It  is  noted 
that  the  temperature  of  the  surrounding  gas  just  a  few  hundredths  of 
an  inch  from  the  boundary  of  the  arc  may  be  nearly  ambient.  Also,  the 
temperature  of  the  electrodes  cannot  exceed  their  boiling  temperature, 
and  hence,  very  large  temperature  gradients  exist  at  the  ends  of  the 
arc  and  at  its  boundaries.  Conduction,  convection,  and  ionization 
potential  affect  the  voltage  gradient  and  arc-column  diameter.  If  the 
ionization  potential  is  lowered  by  the  presence  of  easily  ionized  sub¬ 
stances,  the  voltage  gradient  and  temperature  are  reduced. 
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Heat  is  lost  from  the  arc-column  by  conduction,  convection,  and 
diffusion  through  the  surrounding  gas.  It  is  estimated  that  about  15 
percent  of  the  heat  is  given  off  by  radiation.  Most  chemical  compounds, 
particularly  gases,  are  completely  dissociated  at  arc  temperatures, 
and  contribute  significantly  to  transfer  of  heat.  That  is,  dissociated 
atoms,  ions,  and  electrons  diffuse  from  the  arc  column,  and  when  they 
recombine  in  the  cold  surroundings  they  give  up  their  energy  of  dis¬ 
sociation. 

In  a  survey  of  literature  on  electric  arcs  made  by  Cooper,  et  al 
(29)  in  1950,  they  report  that  Filkelnberg  showed  that  burn-off  of 
carbon  electrodes  is  due  almost  entirely  to  evaporation,  i.e.,  it  occurs 
only  when  the  boiling  point  of  the  carbon  is  reached  in  the  cathode 
spot. 

One  equation  relates  the  voltage  gradient  of  the  arc  to  the  current 
by 

E  =  B/i"  (3.1) 

where 

E  =  voltage  gradient 
i  “  current 
n  and  B  =  constants 

For  nitrogen  and  air,  n  »  0.60,  and  for  steam  n  *  0.59.  The  equivalent 
value  for  the  material  in  the  arc-column  in  molten  or  vaporized  rock 
is  not  known,  but  for  basalt,  granite,  and  quartzite  the  arc  becomes 
stable  after  a  melt  is  formed.  The  above  mathematical  law  is  an  oppo¬ 
site  effect  than  that  described  by  Ohm's  law.  This  negative  character¬ 
istic  requires  a  series  resistance  to  maintain  a  stable  arc. 
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The  current  densities  in  high  current  arcs  reported  to  1949  were 
as  high  as  20,000  atnps/in*.  An  amount  of  heat,  approximately  equal  to 
the  voltage  drop  at  the  electrode  times  the  arc  current,  is  produced 
at  each  electrode,  and  an  additional  amount  from  the  recombination  of 
dissociated  molecules. 

•  Arc  stability,  or  its  ability  to  reignite  after  short  circuit, 
or  after  each  reversal  of  alternating  current,  is  determined  by  the 
arc  characteristics  and  the  recovery- voltage  capabilities  of  the  power 
source.  It  has  been  found  that  superposition  of  high  frequency  volt¬ 
ages  may  add  to  the  stability  of  arcs.  An  arc  is  usually  started  by 
bringing  the  electrodes  in  contact,  or  by  "striking  the  arc",  preferred 
voltages  for  dc  being  60  to  150  volts  and  for  ac,  80  to  200  volts. 

The  kinetics  of  an  arc  formed  in  a  molten  rock  environment  and  the 
associated  heat  transfer  processes  are  very  complex.  Hence,  experi¬ 
ments  will  be  designed  primarily  to  measure  gross  phenomena.  Unless 
the  physical  and  chemical  processes  which  take  place  in  the  rock  near 
the  arc  are  of  such  a  character  that  much  of  the  generated  heat  cannot 
be  recovered,  the  high  temperature  of  the  arc  provides  such  a  large 
temperature  difference  that  heat  flow  is  rapid.  This  factor  plus  the 
heat  transfer  by  recombination  of  atoms,  etc.,  at  the  surface  of  the 
cavity  makes  for  efficient  flow  of  heat  energy  into  the  rock.  Carbon 
arcs  need  not  be  cooled,  a  proce*:s  which  often  represents  a  large 
heat  loss.  Pertinent  information  from  an  abundant  literature  on  arcs 
will  be  obtained  and  utilized  in  support  of  further  research  effort. 

The  only  problems  encountered  to  date  in  the  use  of  carbon  arcs 
are  (a)  the  electrodes  fracture  easily,  and  (b)  it  is  difficult  to 
maintain  stable  arcs  without  high  frequency,  high  voltage  current 
similar  to  that  employed  in  TIG  welding. 
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It  has  been  found  (30)  that  arcs  may  be  initiated  (1)  by  transi¬ 
tion  from  a  low  current  discharge  such  as  a  glow,  (2)  by  means  of  a 
transient  non-steady  spark  discharge,  or  (3)  by  separating  two  current- 
carrying  electrodes  initially  in  contact,  usually  called  striking  an 
arc.  Molecules  of  gases  and  other  substances  are  electrically  neutral, 
and  will  conduct  electricity  only  if  charged  particles  are  introduced 
or  created  in  them.  In  most  common  gases  the  carriers  are  electrons 
and  ions,  while  in  the  complex  environment  created  by  an  arc  inside 
of  a  rock  additional  types  of  carriers  may  be  involved. 

An  arc  channel  in  mercury  vapor  has  been  shown  to  be  at  5900°K 
and  contains  approximately  equal  densities  of  positive  and  negative 
charge.  Several  other  conditions  are  shown  (30)  to  hold  in  arc  be¬ 
havior  including  types  of  equilibria,  number  of  collisions,  etc. 

Under  extreme  conditions  the  temperature  of  the  central  arc  column 
may  exceed  50,000'’K. 

Energy  is  transferred  out  of  the  arc  column  to  a  surrounding  wall 
in  several  ways: 

1)  By  impact  of  molecules,  ions,  or  electrons,  at  high  pressures 
by  thermal  conduction. 

2)  By  ions  and  electrons  diffusing  to  the  wall,  recombining  and 
giving  up  their  energy  of  ionization. 

3)  By  excited  atoms  giving  up  energy  of  excitation. 

4)  By  dissociated  atoms  re-associating  and  giving  up  energy  of 
dissociation. 

There  are  also  several  types  of  radiation  occurring,  Including 
those  of  atomic  and  molecular  spectra,  and  radiation  involved  in 
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volume  recombination.  The  distribution  of  energy  transfer  is  known 
for  mercury  vapor  and  other  gases  as  a  function  of  pressure,  but  no 
results  of  such  studies  for  vaporized  minerals  could  be  found  in  the 
literature. 

In  an  arc  between  carbon  electrodes  in  a  drill  hole  in  granite  the 
arc  must  be  initiated  in  an  atmosphere  of  air.  As  the  temperature  of 
the  rock  increases  it  is  melted  and  vaporized  and  the  arc  environment 
becomes  a  combination  of  melt  and  vapors  from  the  rock. 

It  is  pointed  out  by  Somerville  that  electrode-column  junctions 
have  sharp  discontinuities,  electrically,  thermally,  and  other  ways. 
Electrically,  a  transition  is  made  from  a  conductor  (electrode)  in 
which  current  is  carrie^^  by  electrons  to  one  (gas)  in  which  both  elec¬ 
trons  and  ions  carry  current.  At  the  anode  electrons  may  enter  the 
metal,  but  ions  do  not  move  from  the  metal  to  the  gas,  while  at  tht 
cathode  gaseous  ions  move  to  the  metal,  but  electrons  must  overcome 
a  potential  barrier.  Thermally,  the  arc  column  is  very  hot  and  the 
electrode  relatively  cool.  These  and  associated  problems  are  of  funda^ 
mental  interest  in  heater  design. 

Arc  Experimentation  -  Recent 

In  1965  the  Bureau  of  Mines  (31  )  published  a  comprehensive  review 
of  the  literature  and  some  experimental  results  of  tests  on  high  cur¬ 
rent  metallic  arcs.  Much  of  this  information  is  pertinent  to  the  cur¬ 
rent  project  including  the  Glossary  of  terms. 
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GLOSSARY 

Anode  Fall:  A  very  thin  space-charge  region  in  front  of  an  anode  surface, 
characterized  by  a  steep  potential  gradient  through  the  region. 

Arc  Plasma:  The  space  between  arc  terminals  in  which  gaseous  conduction  of 
electricity  takes  place.  Approximately  equal  numbers  of  electrons  and  ions 
virtually  neutralize  the  space  charge,  but  ionization  and  excitation  are 
generally  intense.  The  potential  gradient  is  low. 

Cathode  Fall:  A  very  thin  space-charge  region  in  front  of  a  cathode  surface, 
characterized  by  a  steep  potential  gradient  through  the  region. 

Cathode  Spot:  A  bright  candescent  spot  on  the  surface  on  a  cathode,  apparently 
an  active  area  through  which  current  flows  from  the  plasma. 

Cathoae-Spot  Mode:  A  particular  mode  of  arc  operation  as  observed  in  laboratories 
of  the  Linde  Company.  The  cathode  spot  is  unusually  bright  and  the  arc  is  ab¬ 
normally  constricted  at  the  cathode  terminus.  The  same  behavior  has  been  called 
the  "contracted  mode." 

Cold-Cathode  Arc:  An  indefinite  term  descriptive  of  any  arc  with  a  cathode 
that  is  not  candescent. 

Electrode  Burn-Off  Rate:  The  rate  at  which  an  electrode  is  consumed  by  an  arc 
in  units  of  mass  per  time  per  arc  power. 

Electrode  Consumption  Rate:  The  rate  at  which  an  electrode  is  consumed  by  an 
arc  in  units  of  mass  per  time  per  arc  current. 

Electrode  Melting  Rate:  The  rate  at  which  an  electrode  is  consumed  by  an  arc 
in  units  of  mass  per  time. 

Globular  Transfer:  Also  known  as  "drop  transfer"  or  more  specifically  as  "large- 
drop  transfer."  The  term  describes  the  transfer  of  metal  as  relatively  large 
drops  or  globules  during  consumable-electrode  arc  welding. 

High-Current  Metallic  Arc:  An  electric  arc  between  metal  terminals  and  at  a 
high  enough  current  to  cause  appreciable  melting  of  at  least  one  terminal, 

High-Intensity  Arc:  Similar  to  a  high-current  metallic  arc,  except  not  neces¬ 
sarily  between  metal  terminals.  The  term  originally  was  used  to  describe 
arcs  of  several  hundred  amperes  between  carbon  electrodes  in  air. 

N-Mode:  A  corruption  of  "normal  mode",  a  term  that  contrasts  with  "cathode- 
spot  mode";  a  particular  mode  of  arc  operation  as  observed  in  laboratories 
of  the  Linde  Company.  The  normal  mode  has  a  normally  diffuse  cathode  spot 
and  no  un  .jal  constriction  at  the  cathode. 

Normal  Arc:  A  snecifically  intended  to  differentiate  between  the  arcs 
that  are  commo.  *  served  and  the  low-pressure  "skittering  arcs." 

Positional  Stabii1t>:  A  Vjsls  employed  at  the  Westinghouse  Research  Laboratories 
to  describe  the  behavio;  of  vacuum  arcs.  The  cathode  spot(s)  of  positionally 
stable  arcs  operate  at  the  tips  of  cathodic  electrodes.  The  cathode  spot(s) 
of  posit ioral ly  unstable  arcs  operate  up  and  down  the  sides  of  cathodic 
electrodes. 
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GLOSSARY  (continued) 


Skittering  Arc:  An  anomalous  low-pressure  arc  form  characterized  by  multiple 
cathode  spots  of  small  size  and  in  rapid  motion.  Skittering  arcs  are  further 
described  as  "stable"  if  the  cathode  spots  operate  at  the  tips  of  cathodic 
eler.trodes,  and  as  "displaced"  if  the  cathode  spots  operate  up  and  down  the 
sides  0^  ''?t‘'odic  electrodes. 

Spray  Transfer:  Also  known  as  "small -drop  transfer."  The  term  describes  the 
transfer  of  metal  as  a  spray  of  droplets  during  consumable-electrode  arc 
welding.  Spray  transfer  is  further  described  as  "axial"  when  the  stream  of 
droplets  is  along  the  arc  axis,  and  as  "rotating"  if  it  rotates  at  an  oblique 
angle  to  the  arc  axis  thus  generating  a  conical  surface. 
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Experimental  equipment  utilized  by  the  Bureau  of  Mines  consisted 
of  a  small  scale  arc  furnace  with  appropriate  auxiliary  instrumentation 
and  controls.  Arc  potential  was  found  to  vary  with  material  of  the 

cathode  and  anode,  gas  pressure,  type  of  gas,  arc  gap,  current,  and  some 

other  factors,  for  a  configuration  utilizing  a  button  anode  and  a  con¬ 
sumable  metallic  electrode. 

The  following  trends  were  noted:  1)  current  densities  may  be  in 
the  order  of  10*  amp/cm*  for  normal  arcs,  and  10*  to  10®  amp/cm^  for 
skittering  arcs,  2)  increase  of  arc  current  results  in  an  increase  of 
cathode-spot  area,  which  decreases  before  it  increases  in  size, 

3)  cathode-spots  for  normal  arcs  diffuse  at  lower  pressures,  4)  cathode- 

spot  densities  are  different  for  different  gases,  5)  the  arc  gap  af¬ 

fects  the  cathode-spot  area. 

The  electrode  consumption  (melting)  rate  of  zirconium  in  argon 
gas  at  250  torr,  520  to  540  amps,  varied  from  4.19  to  4.42  grams/sec 
or  about  0.48  gm/min  per  amp.  The  cathode-spot  temperature  for  thori- 
ated  tungsten  varied  from  2322°C  to  2867®C,  the  latter  for  500  amps 
at  32  volts.  No  data  are  given  for  carbon  electrodes. 

In  early  research  the  cathode-spot  temperature  for  carbon  arcs , 
measured  by  optical  pyrometry,  was  near  3600®K.  Most  of  the  equations 
developed  for  relationships  between  potential,  are  gap,  and  current 
are  for  low  currents  with  electrodes  other  than  carbon.  However, 
they  offer  only  approximations.  The  arc  potential  includes  the  cathode 
fall,  the  anode  fall,  and  the  plasma  potential  component.  There  is 
little  direct  experimental  data  on  these  parameters  for  any  types  of 
electrodes. 
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For  studies  of  current  densities  made  with  water-cooled  copper 
electrodes  and  total  arc  currents  of  100  amperes,  the  current  density 
ranged  from  about  2  by  10^  to  2.5  by  10*  amp/ cm*  at  the  center  of  the 
anode  area,  with  an  approximate  Gaussian  distribution  about  the  center. 

The  temperatures  of  arc  plasmas  are  determined  by  measurements  of 
the  electron  density  from  which  the  temperature  may  be  calculated.  Ioni¬ 
zation  temperatures  are  measured  by  observation  of  line  spectra,  the 
Intensity  of  which  is  related  to  the  temperature.  All  such  measurements 
usually  are  accepted  with  qualification  because  of  several  uncertain 
measurement  and  calculation  factors. 

Maximum  arc  temperatures  are  a  function  of  arc  current  (Fig.  3) 
and  range  from  about  4000®K  to  2  by  10®®K.  It  Is  reported  that  some 
toroidal -discharge  plasma  temperatures  continue  to  Increase  at  higher 
current  reaching  about  2  by  10®  to  4  by  10®®K  at  150,000  amperes,  al¬ 
though  arc  plasmas  may  not  follow  the  same  path.  Also,  if  proper  nuc¬ 
lei  are  present  nuclear  Interactions  may  begin  to  occur  at  80,000  am¬ 
peres  In  the  toroidal  discharges,  yielding  neutrons. 

Axial  temperatures  of  the  arc  are  lower  near  the  anode,  differences 
ranging  from  a  few  hundred  degrees  to  several  thousand  degrees.  Molec¬ 
ular  gases  appear  to  be  related  to  higher  temperature  differentials 
than  atomic  gases.  Although  pressure  has  an  effect  on  arc  performance, 
other  research  In  the  thermal  fragmentation  project  Is  being  conducted 
largely  at  or  near  atmospheric  pressure. 

The  arc  plasma  Is  largely  electrons  and  ions  of  non-uniform  distri¬ 
bution,  and  may  be  classified  as  a  normal  or  skittering  arc  according 
to  mode. 
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Figure  3.  Maximum  Temperatures  in  Arcs  and  Toroidal  Discharges 
as  a  Function  of  Current  (31). 
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El«ctrodtt  consumption  rittt  for  motsl  Ar«  Affoctod  hy  th@  (,'om- 
posltlon  of  the  eloctrodti,  the  ere  currents  the  pressure  end  compo¬ 
sition  of  ere  environmental  gaii  the  electrical  conductance  of  the 
electrodes,  the  arc  gap,  the  anode  size,  and  other  dimensional  fac¬ 
tors.  The  first  three  appear  to  be  the  most  dominant. 

The  characteristics  of  the  cathode,  plasma,  and  anode  regions 
of  the  arc  are  discussed  In  some  detail  by  Wood  and  Beal  (31)  and 
this  Information  will  be  of  Importance  In  projected  research  studies 
cf  the  behavior  of  arcs  In  environments  composed  of  the  vapors  created 
by  an  arc  In  a  closed  hole  In  rock.  It  is  noted  that  even  though  the 
power  input  Into  an  arc  can  be  estimated,  the  expenditure  of  power 
is  determined  by  the  properties  of  the  system,  which  Include  thermal 
conductivity,  emissivlty,  vapor  pressure,  electron  work  function, 
enthalpy,  heats  of  fusion  and  vaporization,  the  melting  and  boiling 
points  of  the  electrode,  and  the  transmittance  and  thermal  conductance 
of  the  gas  surrounding  the  electrode. 
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IV 

Tht  cott  of  •xp»)'1m«ot«ny  ttuidytno  tho  vorloui  PArtmoUt'i  of 
•ng1nt«f1n9  lyitwni  JuitlfUi  <i<wpHfUd  thtorH»tk«)  itodioi  wtHfh 
pravldo  guldinco  for  tht  opttmt)  dtiign  of  lueh  lytttmi.  It  1i  ofttn 
of  vtlut  to  bo  ibU  to  confftm  thoorttlr.oMy  cortoln  ui^unutl  phonomtnt 
which  hovt  boon  obittvod  oMptrlioontilly.  Tht  mochonlci  uf  tho  thonnil 
intchinki)  rock  frigioontotlon  lystom  con«koro(l  horo  oii*  t^ypittol  k 
th«t  full  ikCtlo  tfitt  oro  difficult  tnd  okpthilvo. 

Tht  tctuti  frigmtntitlon  lyittm  k  thi'oo*dlnton«iQn«l ,  Figure  4 
ihowi  i  iwni" Infinite  region  with  tqutlly  sptced  ImIos.  The  holei 
«r«  drived  to  i  conttent  depth  over  the  ero«  where  rock  k  to  be 
removed,  Heat  $ource&  are  placed  at  the  bottoma  of  the  drilled 
holes.  The  reiulting  thennal  Inclualont  (heat<  '  lonea)  cause  two 
types  of  fracturea,  one  of  which  k  along  planes  containing  the  axes 
for  the  rows  of  holes,  A  more  Important  fracture  occurs  parallel 
to  tho  working  face  (perpendicular  to  the  hole  axee)  at  approximately 
the  depth  of  the  Inclusloni,  This  fracture  pattern  makes  possible  the 
removal  of  rectangular  blocks  of  rock,  the  block  dimensions  being 
determined  by  the  hole  spacing  and  the  depth  of  the  thermal  Inclusions. 
To  ficllltate  fracture  end  removal  of  the  rock,  free  surfaces  In  the 
form  of  slots  (Figure  S  )  or  as  would  be  provided  by  a  spiral  round 
are  required.  The  slots  can  possibly  be  cut  by  several  different 
means.  Drilling  tangent  holes  Is  one  method  which  has  been  explored. 
Other  possible  methods  Include  water  Jets,  electron  boams,  plasma 
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Jets,  lasers  and  other  novel  techniques.  Employing  heaters  for 
pyramid  cuts  or  In  a  spiral  round  similar  to  the  way  explosives  are 
used  Is  possible.  However,  the  theoretical  studies  assume  no  dis¬ 
placement  relief. 

A  complete  theoretical  analysis  of  the  mechanics  of  the  process 
required  a  transient  temperature,  thermal  stress  and  fracture  study 
of  the  three-dimensional  geometry  shown  In  Figure  4.  The  problem 
Is  nonlinear  in  the  sense  that  the  thermal  and  elastic  properties  of 
nK)st  rocks  are  highly  dependent  on  temperature. 

To  study  the  fracture  or  fragmentation  problem  two  alternatives 
to  the  method  selected  have  been  considered.  The  first  and  best 
alternative  from  the  viewpoint  of  accuracy  would  have  been  to  study 
the  actual  three-dimensional  geometry  with  anisotropic  properties  and 
temperature  dependent  boundary  conditions.  From  the  three-dimensional 
geometry  the  finite  element  method  can  be  used  to  study  most  of  the 
parameters  which  affect  fracture.  Three-dimensional  codes  capable  of 
performing  the  calculations  would  have  had  to  be  developed  at  much 
greater  cost.  The  second  alternative  would  have  been  to  treat  a 
single  hole  and  heater  in  a  semi  or  half  space.  If  the  interactions 
between  the  Individual  stress  fields  can  be  neglected,  the  stress  at 
any  point  in  the  rock  can  be  found  by  adding  the  contributions  to 
the  total  stress  at  that  point  of  each  of  the  surrounding  holes. 

The  approach  used  allows  the  displacement  boundary  conditions  to 
supply  the  interaction  between  the  stress  fields.  Isotropic  material 
properties  are  used  and  the  effect  of  temperature  on  the  thermal 
properties  has  been  considered  only  In  the  one-dimensional  studies. 
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For  the  two-dimensional  studies  a  particular  rock  (Dresser  basalt) 
was  selected  whose  thermal  properties  are  relatively  Independent  of 
temperature.  Consider  the  hole  pattern  shown  In  Figure  6.  This 
pattern  can  be  assumed  to  extend  Indefinitely  In  all  directions.  The 
dashed  lines  Indicate  the  Intersections  with  the  working  face  of  the 
planes  of  symmetry.  Over  each  of  these  planes  of  symmetry,  material 
displacements  do  not  take  place  perpendicular  to  the  plane.  Thus, 
considering  the  location  of  these  planes  the  three  models  shown  In 
Figure  7  can  be  used  to  study  the  significant  process  parameters. 
Figures  8-10  show  one  one-dimensional  and  three  two-dimensional  models 
which  have  been  used  In  this  study.  The  Influence  of  hole  diameter, 
depth  and  spacing  has  been  explored  using  the  two-dimensional  models. 

In  these  studies  the  elastic  properties  have  been  treated  as  functions 
of  temperature.  The  temperature  distributions  were  obtained  assuming 
constant  thermal  properties.  The  one-dimensional  models  have  been 
used  to  study  appropriate  types  of  boundary  conditions  as  well  as  the 
effect  of  temperature  dependent  thermal  properties. 

The  slot  model  has  been  used  to  study  the  fracture  that  occurs 
parallel  to  the  working  face  of  the  tunnel  or  excavation.  In  the  slot 
model  such  a  fracture  passes  generally  along  the  bottom  of  the  slots. 
This  fracture  depends  on  the  nature  of  the  thermal  Inclusions.  The 
model  Is  a  good  Indicator  of  the  effect  of  hole  spacing  as  well  as  the 
effect  of  variations  In  the  heated  length. 

The  hole  model  has  been  used  to  study  the  fracture  that  occurs 
between  holes  on  a  plane  containing  the  hole  axes.  The  model  provides 
Information  about  the  effect  of  different  hole  diameters  and  spacings. 
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Figure  7.  Two  Dimensional  Models 


Axisyninetric  Cylinder  Models 
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The  cylinder  model  has  been  treated  both  as  one  and  two- 
dimensional.  The  one-dimensional  cylinder  model  was  selected  to 
obtain  more  basic  Information  such  as  the  effect  of  nonlinear  properties 
and  the  proper  type  boundary  conditions  for  different  types  of  heaters. 
The  two-dimensional  cylinder  model  In  this  study  was  used  to  verify 
some  of  the  conclusions  obtained  from  the  slot  and  hole  models. 


.if. 
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Chapter  V 

ROCK  CHARACTERISTICS  AND  APPROXIMATIONS 

Rocks  are,  in  general,  complex  materials  with  inherent  inhomo¬ 
geneity  and  anisotropy  as  well  as  temperature-dependent  material 
properties.  For  reck  related  studies  one  must  also  consider,  among 
other  factors,  effects  of  porosity,  moisture  content,  microcracks, 
bedding  and  joints,  stress  relief  and  stress  absorption.  An  exact 
thermal  stress  analysis  with  all  of  these  factors  will  be  ext  mely 
difficult  even  for  the  simplest  one-dimensional  geometry.  Simpli¬ 
fying  assumptions  are,  therefore,  necessary  in  order  to  obtain  a 
workable  solution.  A  description  of  general  characteristics  of 
different  kinds  of  rock  is,  thus,  in  order. 

A.  General  Classification  (32) 

The  most  general  rock  classification  system  is  based  on  the  mode 
cf  origin.  Rocks  are  divided  into  three  main  groups:  igneous  rocks, 
sedimentary  rocks,  and  metamorphic  rocks. 

Igneous  rocks  originate  from  solidification  of  hot,  molten 
material  below  the  earth's  crust.  These  are  very  hard,  massive  rocks 
with  granular  nonporous  structure  and  exhibit  very  high  compressive 
strength.  Their  tensile  strength  is  many  times  smaller  than  their 
v.ompressive  strength,  yet  their  tensile  strength  is  higher  than  the 
tensile  strength  of  most  other  rocks. 

The  most  predominant  characteristic  of  sedimentary  rocks  is 
stratification  or  bedding.  As  a  result,  they  possess  strong  trans¬ 
versely  isotropic  properties,  and  have  much  less  compressive  strength 
than  the  igneous  rocks. 
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The  term  "metamorphism"  denotes  an  aging  process  which  brings 
about  a  complete  change  in  the  original  rock  characteristics  as  a 
result  of  very  long,  continuous  action  of  pressure,  temperature,  and 
moisture  and  chemical  reactions  of  atmospheric  gases.  Most  meta- 
morphic  rocks  possess  the  highly  crystalline  texture  of  igneous  rocks 
with  a  parallel  structure  which  closely  resembles  the  stratification 
of  sedimentary  rocks.  Fracture  generally  occurs  due  to  cleavage 
along  these  planes.  Metamorphic  rocks  that  do  not  possess  the 
cleavage  planes  usually  have  hard,  massive  structure  with  compressive 
strength  on  the  same  order  as  granite.  Most  hard  metamorphic  rocks 
exhibit  some  characteristics  resembling  closely  those  of  igneous 
rocks.  Quartzite  is  a  well  known  example  of  this  type. 

Since  this  investigation  is  concerned  with  thermal  fragmentation 
of  hard  rocks,  mathematical  models  will  be  based  on  properties 
characteristic  of  igneous  rocks. 

B.  Nonlinearity  Considerations 

As  mentioned  in  the  beginning  of  this  chapter,  approximations  are 
necessary  in  order  to  obtain  workable  solutions  for  the  two-dimensional 
analytical  models.  These  approximations  will  be  based  on  the  following 
three  major  factors: 

i.  influence  on  fracture  stresses, 

11.  complications  involved  in  theoretical  studies,  and 
iii.  accuracy  of  available  property  data. 

Fracture  stresses  are  influenced  differently  by  different  non¬ 
linearity  related  factors.  Material  nonhomogeneity,  for  example, 
will  cause  Internal  loading  whereas  a  work-hardening  stress-strain 
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condition  will  load  to  stnass  Absorption.  A  «1ni1Uf  i1  tuition  U 
tru«  regardlnp  thr  dogroo  of  cotnplvxl f.y  of  tbtoretlcil  inilysli. 

Thu$»  for  oxample.  u&ing  finito  olwnont  mathodi,  mittrlil  nonhonto* 
geneltv  and  anisotropy  can  b«  handlod  r«thAr  oislly,  wborois  tb® 
study  of  crack  affects  or  the  Influence  of  rock  porosity  will  require 
a  prohibitively  largo  ainount  of  work.  This  Is  true  regardless  of  the 
fact  that,  for  a  particular  kind  of  nonlinearity,  certain  methods 
are  better  than  others. 

Availability  and  accuracy  of  t'ock  properties  are  perhaps  two 
of  the  most  Important  factors  in  the  consideration  of  theoretical 
analysis.  Review  of  laboratory  techniques  for  measuring  rock  properti 
Indicates  that  there  are  no  general  specified  standards  except  for 
a  very  few  mechanical  property  measurements  133-34).  Besides,  rocks 
are  known  to  exhibit  slightly  different  properties  In  situ  (35).  Also, 
the  characteristics  of  a  given  rock  type  from  different  geological 
locations  vary  considerably.  Most  measured  properties,  as  a  result, 
show  about  ±B%  deviation  from  tlie  average  value.  This  deviation,  how¬ 
ever,  is  considered  to  be  rather  small  for  engineering  purposes,  and 
Iv  usually  neglected. 

The  above  considerations  are  used  in  approximating  the  following 
nonlinear  rock  characteristics; 

1.  Stress-strain  Behavior 

For  most  hard,  crystalline  rocks,  the  stress-strain  curve  Is 
approximately  linear  and  ends  abruptly  indicating  a  brittle  failure 
(36).  Some  of  the  rocks  do  exhibit  a  nonlinear  behavior.  However, 
this  nonlinearity  is  attributed  to  the  presence  of  microcracks  which 
close  under  confinement  (37-40),  The  stress-strain  behavior  of  in 
situ  rock  can,  therefore,  be  assumed  to  be  linearly  elastic. 


And  (tniN  (4))  fondut'Ud  An  i,Nt»nA<VA  ^tudy  on  th4 
•ffnct  of  t«MA|m»'Atu»'n  en  *'«Utiuniih<oA  t'f  a  vArInty  of 

rocKi.  AccoiMlIna  to  tholo  noAuUii  moil  t  ockA  aHom  a  itOAdly  dnuooAie 
In  ylold  |)o(ni  with  Inct'OAiA  In  tpi^oyAlwon.  ThU  of  foot,  howov»r*  1i 
not  ImportAnt  in  tho  itudy  of  nock  fAUuno  OAsgUIno  fnow  th*rwA) 
Incluilon,  ThU  1i  duo  t('  a  numbon  of  ntAionA.  Mrttt  the  HtreiA- 
Atriln  uunvei  for  rockt  AiHt  uAiiAMy  obtAlned  unden  unlAMUl  cowpreAilon, 
whereAi  fnACture  InItlAtes  from  tenille  itrekiot  thet  ore  mony  timei 
imtlUr  thAn  the  yleUt  strength  In  compretolon.  Alio,  during  the 
prtllmlnAry  field  teiti  conducted  In  thli  InveitIgAtlon,  only  a  imeU 
percent  of  the  frectured  rock  voluine  wai  obnervud  to  hAve  experienced 
Any  ippredeble  chAnge  in  tempereture.  Thus,  the  effect  of  the  lowered 
yield  itrength  Is  locAllied  In  nAture,  In  the  close  vicinity  of  the 
Inclusion,  However,  very  high  tetnperetures  In  this  region  diAstlcAlly 
reduce  the  elestlc  resUtence  of  the  meterlAl, 

2-  H9I!mn£l-^.Ah.<1  .UftlrW 

Herd,  crystalline  rocks  consist  moinly  of  guArti,  feldspAr,  Augite. 
inAgnetIte,  And  micA.  For  a  given  rock  type,  eAch  constituent.  In 
yenerAl,  differs  from  others  In  both  texture  And  grain  slje.  At  the 
same  time,  average  grain  sixe  of  a  constituent  varies  with  rock  types. 
Most  hard,  crystalline  nonporous  rocks  posse'^s  flne-tO'medium  grained 
texture  end  the  normal  range  of  grain  sixe  Is  usually  from  0,03  mm  to 
1.0  mm.  Individual  crystals  are  chemically  heterogeneous  and,  in 
general,  possess  different  physical  properties.  However,  the  grain 
distribution  Is  very  random.  As  a  result,  on  a  macroscopic  scale, 
hard  rocks  usually  have  homogeneous  properties.  Any  slight 


h»tiir0i)0ni(ty  wit),  of  (TOuriiOt  tw  <iuiU  iniignificont  for  the  in  «itu 
rock  iny«KtiQAtion>i  (J6). 

Anitotropir  hehevior  ii  pf^nouneed  only  fur  ledimentiry,  Ktretlfied 
I'ockt  ind  metemornhie  »'ucki  which  h«ve  we))  defined  cteevege  p)inei. 
Thete  ruckt  ere  generetiy  toft  to  medium  herd,  end  ere  not  inciuded 
ill  thit  ene)ytifi,  For  herd  I'ochi^  propertiet  ere  ettumed  to  be 
i«otropic<  At  e  metier  of  feet,  ieboretory  meeiurtMnenit  ere  genere))y 
obteined  in  rendom  directions  end  the  properties  ere  then  evereged, 
without  recording  the  veriotions  in  individue)  directions. 

Creeks  ere  creeled  in  the  form  of  voids  es  e  resuit  of  remove)  of 
getes  end  weter  vepor  which  ere  entrepped  et  high  pressure  end  temp* 
ereture  during  I'ock  formation.  These  voids  ere  extremeiy  sma))  com¬ 
pered  to  feuits  end  Joints  end  ere  usueiiy  in  the  fonn  of  mici'ocrocks 
in  nonporous,  herd,  dense  rocks.  Nevertheiess,  they  have  very  definite 
effects  on  rock  properties,  due  meiniy  to  meterie)  discontinuity  and 
the  pressure  of  eccumuieted  moisture  in  the  cracks.  For  example,  the 
initia)  noniinearity  In  the  eiestic  behavior  of  rock  has  been  shown 
to  be  due  to  the  existing  microcracks  (36-40). 

Compressive  sti'ength  of  rocks  is  observed  to  decrease  with 
increasing  moisture  content.  Aithough  not  compietely  understood,  this 
weakening  effect  is  be)ieved  to  resu)t  from  the  interna)  )oeding 
caused  by  pore  pressure.  Apart  from  this  weakening  effect,  cracks 
act  as  stress  raisers  and  have  a  very  significant  effect  on  fracture 
propagation. 

As  the  equations  of  thormoelastici ty  are  based  on  the  theory  of 
continuum,  crack  effects  are  probably  the*  most  difficult  to  study. 
However,  for  comparatively  nonporous,  hard  rocks,  cracks  occupy  less 
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than  \%  of  toll)  volutne  as  compartd  to  10-20X  for  highly  porous 
rocks  likt  sands ton*  and  llmastona  (42).  Mtncf«  for  hard  rocks  which 
posstss  a  vary  low  apparant  porosity,  tha  mkrocrack  af facts  can  be 
naglactad. 

^ •  Ttffiptratura  Effacts 

Cyan  though  application  of  haat  as  a  technique  to  fracture  rocks 
has  baan  usad  slnca  antiquity,  llttla  was  known  about  the  mechanisms 
Involved.  Howavar,  In  recant  years,  a  large  amount  of  effort  has  baan 
put  forth  on  the  part  of  numerous  researchers  to  Investigate  the 
mechanisms  responsible  for  thermal  fragmentation  of  rock  (4.  14.  43-47). 
Heating  devices  being  considered  Include  such  modern  techniques  as  high 
frequency  electric  heating.  Induction  and  microwave  heating,  electric 
ar'cs,  Infra-red  rays,  electron  beam,  lasers,  and  plasma  Jets.  Although 
parameters  such  as  the  heat  transfer  efficiency,  energy  Input  rate, 
and  fracture  time  vary  considerably  with  the  mode  of  heating,  the 
fracture  mechanism  Involved  remains  unchanged. 

As  opposed  to  the  mechanical  fragmentation  mode,  the  thermal 
fragmentation  mechanism  depends  on  the  following; 

1.  Variation  of  physical  properties  with  temperature, 

11.  Thermal  weakening  effects  as  a  result  of  spallation.  Inter¬ 
granular  crack  growth,  chemical  changes,  and  expansion  of  entrapped 
gases  and  moisture, 
a .  Thermoelastic  Properties 

Properties  of  importance  In  the  thermal  fragmentation  analysis 
of  rocks  Include  thermal  conductivity,  diffusivity,  coefficient  of 
thermal  expansion,  modulus  of  elasticity,  and  Poisson's  ratio.  For 
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most  rocks,  these  properties  are  highly  temperature  dependent.  This 
Is  to  be  expected  since  the  properties  of  different  rock  constituting 
minerals  exhibit  varying  degrees  of  temperature  dependency.  Formulating 
the  closed  form  thermoelastlclty  equations  and  their  solution  for  mater¬ 
ials  with  temperature  dependent  properties  Is  an  extremely  difficult 
task.  However,  using  numerical  finite  element  methods,  these  effects 
can  be  studied  without  a  great  amount  of  difficulty. 

Thermal  conductivity  of  most  Igneous  rocks  decreases  with 
Increasing  temperature.  Figure  11  shows  the  variation  of  thermal 
conductivity  for  some  typical  Igneous  rocks  as  obtained  by  Birch  and 
Clark  (48). 

Diffusivlty,  In  general.  Is  obtained  mathematically  from  the  values 
of  conductivity,  specific  heat,  and  density.  Most  hard,  non-porous  Ig¬ 
neous  rocks  show  little  variation  In  density  with  temperature. 

Mean  value  of  specific  heat  over  the  temperature  range  from  25°C 
to  625°C  was  obtained  by  Geller,  et  al,  (25).  Only  a  small  variation 
was  observed  among  the  specific  heats  of  fourteen  rock  types  tested 
and  the  values  ranged  from  approximately  0.23  to  0.25  cal/gm-°C. 

However,  no  investigations  were  made  to  study  the  variations  in  specific 
heat  with  temperature. 

Temperature  dependency  of  specific  heats  for  six  hard  rock  types 
was  later  investigated  by  Lindroth  and  Krawza  (49).  Tests  were  con¬ 
ducted  at  temperatures  up  to  1000°C.  According  to  this  study,  tem¬ 
perature  dependence  of  specific  heat  Is  strongly  influenced  by  the 
percentage  of  quartz  content.  Results  of  this  study  are  shown  In 
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Sy  Syenite,  Ontario  RG 

WG  Westerly  Granite  Br 
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Figure  11.  Thermal  Conductivity  of  Hoi ocrystal line  Rocks  (48) 
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Figure  12  .  The  break  1n  these  curves  occurs  as  a  result  of  a 
phase  change  of  quartz  at  approximately  573®C.  For  rocks  containing 
less  than  5%  quartz,  this  break  occurs  at  higher  temperatures. 

Similar  discontinuities  can  be  observed  in  curves  shown  in 
Figure  13,  which  shows  the  plot  of  coefficient  of  thermal  expansion 
for  hard  crystalline  bodies  as  a  function  of  temperature.  These  are 
part  of  the  results  obtained  by  the  Canadian  Mines  Branch  of  the 
Department  of  Mines  and  Mineral  Surveys  which  conducted  experiments 
on  37  rock  types  to  determine  the  variation  of  linear  thermal  expansion 
with  temperature  (25).  Studies  on  fusion  were  also  performed. 

Fusion  temperatures  of  45  rock  types  were  obtained;  fusion  temperatures 
for  hard  crystalline  rocks  were  observed  to  lie  in  the  range  from  1150°C 
to  1300°C. 

Studies  conducted  on  the  effect  of  temperature  on  elastic  properties 
of  rocks  indicate  that  for  hard,  nonspallable  rocks,  both  Young's 
modulus  and  Poisson's  ratio  show  a  marked  decrease  in  the  vicinity  of 
the  fusion  temperature.  For  spallable  rocks,  however,  the  upper 
temperature  limit  is  restricted  by  the  alpha  to  beta  phase  transition 
of  quartz,  and  above  this  temperature  range,  only  small  variations  in 
elastic  properties  are  observed  (50-52). 
b.  Fracture  Mechanisms 

As  opposed  to  the  well  defined  properties  involved  in  both  the 
temperature  and  stress  analyses,  the  fracture  phenomenon  depends  on  a 
number  of  mechanisms  for  which  the  variations  with  temperature  have 
not  been  described  in  tabular  or  graphic  form.  Little  is  known  about 
the  interactions  of  various  mechanisms  involved.  Nevertheless,  it  is 
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Figure  13.  Average  Linear  Expansion  of  Igneous  Rocks  (25) 
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well  established  that  temperature  has  definite  weakening  effects  on 
rock  strength,  and  the  mechanisms  responsible  include,  primarily, 
spallation,  intergranular  crack  growth,  chemical  changes,  and  gas  and 
water  pocket  expansion. 

The  relative  importance  of  each  of  these  weakening  factors 
depends  mainly  on  the  nature  of  rock;  porosity,  average  grain  size, 
microcracks,  and  mineral  constituents.  Except  for  spallation, 
relatively  little  information  Is  available  on  temperature  effects  on 
other  mechanisms  involved  in  fracture,  mainly  due  to  their  inseparation. 

The  term  "spallation”  refers  to  a  progressive  failure  of  rock  in 
the  form  of  chips  caused  by  thermal  stress.  The  basic  mechanism  involved 
in  formation  of  a  spall  depends  on  a  sudden  increase  in  the  thermal 
expansion  of  quartz  due  to  its  phase  transformation  at  approximately 
573^C.  Although  it  is  well  established  that  a  certain  percentage 
of  quartz  and  its  phase  transition  are  essential  to  the  inducement  of 
spalling,  review  of  various  research  efforts  on  thermal  spalling 
reveals  that  little  agreement  exists  as  to  the  exact  cause  of  the 
phenomenon  (  25,  33,  47,  53-55). 

According  to  Gray  (56),  spalling  can  be  controlled  by  proper 
choice  of  heating  rate.  For  high  heating  rates  which  induce  surface 
melting,  spalling  effects  become  trivial.  Studies  performed  by 
Moavenzadeh,  et  al,  (46)  indicate  that  for  igneous  rocks  complicated 
crack  growth  mechanisms  have  little  or  no  effect  on  fracture  compared 
to  those  due  to  the  thermal  stress  field.  Thus,  the  fracture  of  rock 
can  be  predicted  with  reasonable  accuracy  from  the  knowledge  of  the 
stress  field  alone. 
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C-  Fracture  Theories 

As  discussed  above,  in  thermally  Induced  rock  fragmentation 
studies,  the  thermal  stress  field  can  be  regarded  as  the  single  major 
factor  responsible  for  rock  failure.  Numerous  theories  have  been 
proposed  over  the  years  to  relate  the  stress  field  to  the  brittle 
failure  of  materials.  Of  the  various  theories,  the  Griffith  theory, 
in  its  different  modified  forms,  has  been  the  most  widely  accepted 
in  the  field  of  rock  mechanics. 

1.  Griffith  Theory 

Griffith's  theory  evolved  through  his  study  of  the  problem  of 
the  tensile  strength  of  glass  being  much  lower  than  theoretically 
expected  (57).  He  suggested  that  the  low  tensile  strength  was  due  to 
failure  caused  by  stress  concentrations  at  the  tips  of  minute  internal 
and  surface  ‘laws  whicn  have  come  to  be  known  as  Griffith  cracks.  In 
formulating  this  theory,  Griffith  assumed  the  material  to  be  isotropic, 
homogeneous  continuum  containing  randomly  oriented  sharp  ended  elliptical 
cracks.  The  mathematical  condition  for  the  initiation  of  a  macro¬ 
crack  was  derived  based  on  the  energy  approach  which  states  that  the 
work  done  by  externally  applied  forces  is  equal  to  the  sum  of  internal 
strain  energy  and  the  surface  energy  associated  with  the  rupturing 
of  atomic  bonds  when  cracks  are  formed. 

Specific  energy  is  a  difficult  quantity  to  evaluate  through  direct 
measurements.  Hence,  the  fracture  criterion  is  expressed  mathematically 
in  terms  of  principal  stress  values.  If  compressive  stresses  are  taken 
as  negative,  the  conditions  of  fracture  initiation  are  given  by  the 
following  equations. 
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If 


3a  1  +03  >  0 


(5.U) 


fracture  Initiation  occurs  when 


Oi  =  Ot  ; 


If 


3oi  +03  <  0, 


the  condition  of  fracture  Initiation  Is  given  by 

(aj  -  03)2  +  8o^(ai  +  03)  =  0  . 


{5.1b) 


(5.2a) 


(5.2b) 


In  these  equations 


Oj  =  major  principal  stress, 

Oj  =  minor  principal  stress,  and 

o.  =  uniaxial  tensile  strength  of  the 
^  material. 

For  failure  governed  by  Equation  5.Ja  ,  the  microscopic  crack 
begins  to  propagate  In  Its  own  plane  In  a  direction  perpendicular  to 
that  of  the  major  principal  stress.  For  failure  conditions  given  by 
Equation  5.2a  ,  however,  the  crack  begins  to  extend  In  a  plane  at  an 
angle  e  from  the  minor  principal  stress  axis,  given  by 

(5.3) 

It  should  be  noted  that  the  fracture  criterion  given  by 
Equations  5.1  a  and  5.2a  was  developed  for  a  biaxial  stress  state. 
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Sack  (58)  extended  Griffith's  theory  to  three  dimensions  by  con¬ 
sidering  a  penrvy-shaped  crack  under  a  tri axial  stress  state,  and  concluded 
that  the  fracture  initiates  as  a  result  of  growth  of  cracks  that  are 
parallel  to  the  intermediate  principal  stress  direction.  Thus,  the 
intermediate  principal  direction  is  the  most  critical  crack  orientation. 
However,  the  magnitude  of  this  principal  stress  component  has  no  ap¬ 
preciable  influence  on  the  crack  growth  and  the  Griffith  biaxial 
fracture  criterion  given  by  Equations  5* la  5.3  can  be  applied  to 
triaxial  stress  conditions  as  well.  This  result  has  experimentally 
been  verified  by  Brace  (59)  for  igneous  rocks  such  as  granite. 

2.  McClintock-Walsh  Modification 

Experimental  verification  of  Griffith's  criterion,  however, 
brought  to  attention  some  serious  flaws  in  Equations  5.&  and  5.3. 

Under  uniaxial  compression  (oi=0,  03  <  0),  the  formula  loses  its 
meaning  as  it  implies  that  the  compressive  strength  of  the  material 
is  equal  to  eight  times  its  tensile  strength.  For  most  rocks, 
however,  this  is  not  true  as  the  compressive  strength  is  much  higher 
than  eight  times  that  in  uniaxial  tension.  Also,  under  uniaxial 
compression,  according  to  Equation  S3  ,  the  most  critical  cracks 
would  be  those  at  45°  to  the  stress  direction.  However,  experimental 
studies  by  Brace  (59)  and  Bieniawski  (60)  show  that  under  a  compressive 
stress  state,  cracks  propagate  out  of  their  plane  in  the  direction  of 
the  major  principal  compressive  stress.  These  discrepancies  required 
some  modification  of  the  Griffith  theory  for  rock  failure  when  one  or 
two  of  the  principal  stress  components  are  compressive. 
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Experimental  studies  by  Brace  (59)  and  Bieniawski  (60) 
served  to  give  an  insight  into  the  more  complex  mechanism  involved  in 
rock  failure.  Based  on  their  experiments,  they  concluded  that  when 
compressive  stresses  are  present,  pre-existing  Griffith  cracks  close 
before  the  tensile  stress  at  the  crack  tip  reaches  the  critical  value 
for  fracture  initiation.  To  open  these  closed  cracks,  higher  stresses 
are  required  in  order  to  overcome  the  shear  resistance  resulting  from 
the  interlocking  of  irregular  crack  faces.  The  concept  of  crack 
closure  was  also  able  to  explain  the  initial  nonlinearity  in  the 
stress-strain  curve  of  rocks. 

McClintock  and  Walsh  (61)  were  the  first  to  modify  the  Griffith 
criterion  with  consideration  given  to  the  effects  of  crack  closure. 

The  fracture  criterion  in  mathematical  form  given  by  McClintock 
and  Walsh  is  as  follows: 

If 

'^n  “  7  (oi-t»3)  cos  <  0  ,  (5.4a) 


fracture  Initiation  occurs  when 


03  =  01 


/I  +  -  U  ^ 


(5.4b) 


where 


o  is  the  normal  stress  acting  aero::  tiie  crack  faces 
"  causing  the  cracks  to  close, 

0^  is  the  uniaxial  compressive  strength  of  the  material, 

\i  is  the  internal  coefficient  of  friction  between  the  crack 
faces ,  and 

0.,  0-  are  the  major  and  minor  principal  stress  components,  as 
before. 


For  this  co  lition  of  crack  initiation*  crack  propaoatlon  diroc- 


tion  is  given  by 


tan  '  1/n  (5,5) 

where  is  measured  tfout  the  minor  principal  stress  ajds. 

Comparison  of  Equations  B.2  and  5.4  shows  that  the  McClintock- 
Walsh  criterion  is  based  on  the  actual  compressive  strength  of  the 
material  and  also  takes  into  account  the  effect  of  crack  closure. 

Also,  note  that  th**  original  Griffith  criterion  given  by  Equations 
5.1-  5.3  is  valid  whenever  >  0,  so  that  the  normal  stress  acting 
across  the  crack  does  not  cause  crack  closure. 

It  should  be  ’"“alized  that  both  Griffith  and  McClintock-Walsh 
modified  Griffith  c  ...eria  refer  to  fracture  initiation  only,  which  is 
not  the  same  as  fracture  (62).  They  do  not  consider  factors  such  as 
the  energy  of  plastic  defonaation;  orientation,  density,  and  inter¬ 
actions  of  existing  microcracks,  difference  between  the  stress  levels 
causing  fracture  initiation  and  those  causing  the  ultimate  strength 
failure*,  crack  propagation  velocity  and  other  dynamic  effects. 
Nevertheless,  experimental  results  obtained  by  Hoek  and  Bieniawski 
(63)  for  a  wide  variety  of  rocks  show  a  remarkable  agreeniont  with  those 
predicted  theoretically  by  the  Griffith  and  McClintock-Walsh  modified 
Griffith  criteria.  They,  however,  replaced  the  internal  crack  friction 
coefficient,  v,  by  fracture  surface  coefficient  of  friction,  which 
was  obtained  experimentally  for  each  rock  type  tested.  For  igneous 
rocks,  the  value  of  has  been  found  to  lie  between  1.0  and  1.5  (39,  59} 
D.  Properties  Used  in  Analytical  Stuf^ies 

Even  though  field  tests  have  been  conducted  on  Missouri  red 
granite  for  the  theoretical  analysis.  Dresser  basalt  was  found  to 


fo 

bti  moiti  tu<t«b)«  It;  contAtni  or  ho  ()UArti,  AUo,  tho 

theimi)  AHd  »tA«tk  proportloi  for  hii»AU  ai  a  fuHCtton  of  tOAiporAturo 
woro  rpAtJlly  AVAUibk  (14*  4?,  j)k\  64). 

>-  JiS3eAWMciL,An^ 

Tho  AVillNbU  f\nit9  •lomont  twiiporAtui'o  cv>d«  U  not  rapAblo  of 
hondllng  toinptr«tur«  (iopondtnt  thorwa)  propert1p».  Howovor,  for 
OrtAior  b«»AU.  th«  them\«l  conductivity  «nd  dlffuilvlty  only 
smill  vAf'lAtlons  with  tixnporaturo*  is  shown  In  PlguiHts  14  ind  is. 

Thft  uppnr  limit  of  700^C  for  th«  tomporaturp  rang*  was  chostn  In  view 
of  the  fact  that  the  modulus  of  dlai*tklty  for  tetnperatures  above 
700*^0  becomes  very  small,  as  seen  frotn  Table  II  .  The  effect  of  this 
low  value  of  Young's  nKtdulus  Is  that,  regardless  of  the  steepness  of 
temperature  gradients.  In  the  region  whet'e  temperatures  are  higher  than 
700®C,  the  stresses  will  bo  very  small.  In  other  words,  le  upper 

limit,  700®C,  of  the  temperature  range  for  curves  shown  In  Figures 
14  and  15  Is  taken  as  the  temperature  at  which  the  rock  becomes 
plastic.  Hence,  the  non! Ineari tits  In  the  thermal  pi-^jpertles  will 
have  a  very  small  effect,  If  atiy,  on  the  stresses  which  contribute 
to  fracture. 

In  order  to  simplify  the  analysis,  the  assumption  of  temperature 
independent  properties  has  been  employed  In  most  of  the  rock  related 
Investigations.  Under  this  assumption,  properties  corresponding 
to  the  average  temperature  of  the  material  are  usually  selected.  This 
procedure  hes  been  shown  to  yield  analytical  results  well  within  the 
range  of  experimental  errors  (45,  65-68). 

For  this  analysis,  average  values  of  thermal  conductivity  and 
diffusivlty  were  obtained  from  Figures  14  and  15.  The  input  properties 
for  the  finite  element  code,  however,  are  the  conductivity,  specific 
heat  and  specific  gravity.  Even  though  the  specific  heat  values  as  a 


Figurt  14,  Virlitlon  of  Conductivity  of  Dresser  basalt  with  Temperature  (47) 
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Figure  15.  Variation  of  Diffusivlty  of  Dresser  Basalt  with  Temperature  (47) 
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function  of  ttmperaturo  w«ro  rtodlly  iv«11it)1§  for  bisiU  (49) >  th« 
following  |)roctdur«  wts  found  to  b«  nwro  offlclont.  The  ipecific 
grevlty  of  the  materiel  was  assumed  to  he  constant,  equal  to  that  at 
room  temperature.  The  average  value  of  the  specific  heat  was  then 
obtained  from  the  average  values  of  the  thermal  conductivity  and  the 
diffusivlty  using  the  well-known  relationship  (69) 

* 


where 

•  average  value  of  themal  diffusivlty 
k  •  average  value  of  thermal  conductivity 

•  the  specific  gravity  of  material,  and 

c  •  the  average  value  of  the  specific  heat. 

Values  of  the  properties  used  for  the  temperature  analysis  are 
tabulated  below. 

TABLE  6 

Property  Data  Used  In  Temperature  Analysis  of 


Dresser  Basalt  (47.  64,  70) 

Thermal  conductivity,  k,  cal/cm-sec  -  ®C  0.0042 

Specific  heat,  c,  cal/gm  -  °C  0.293 

Density,  p,  gm/cm^  2.97 

Melt  Temperature,  T^^  J250°C 

Surface  convection  coefficient,  h,  cal/cm^  -  sec  -  ®C  0.00021 


* 

For  solids,  no  distinction  Is  generally  made  between  specific  heat  at 
constant  pressure,  c.,  and  that  at  constant  volume,  c^.  Nevertheless, 
it  should  be  noted  that  in  this  equation  c  refers  to  c^. 
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2*  Stress  Analysis 

Unlike  the  temperature  analysis »  the  stress  analysis  was 
performed  using  temperature  dependent  material  properties.  This 
was  necessary  due  to  the  fact  that  the  variations  of  the  stress 
related  properties  with  temperature  greatly  Influence,  as  explained 
later  on,  the  resulting  stress  field  and  consequently,  the  fracture 
predictions. 

The  stress  related  properties  Include  the  coefficient  of  thermal 
expansion,  a,  modulus  of  elasticity  (Young's  modulus),  C,  and  the 
Po1ss„<n's  ratio,  v.  Table  7  gives  the  values  of  these  properties 
used  In  the  stress  analysis. 

TABLE  7 

Property  Data  Used  In  Stress  Analysis  of 
Dresser  Basalt  (4,  14,  52) 


Temperature,  T 

Young's  Modulus 

Poisson's 

Ratio 

Coeff.  of  Thermal 
Expansion,  a. 

°C 

E,  10®  psi 

V 

10'®/°C 

100 

14.5 

0.24 

2.6 

162 

14.1 

0.23 

6.0 

287 

13.5 

0.21 

8.2 

412 

12.6 

0.19 

10.2 

537 

10.7 

0.145 

n.2 

610 

8.3 

0.105 

11.7 

630 

6.4 

0.09 

1,1.8 

650 

4.6 

0.07 

11.9 

670 

2.8 

0.05 

12.0 

690 

0.9 

0.02 

12.0 
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As  was  mentioned  earlier,  the  rock  is  assumed  to  become  plastic 
at  700°C  as  its  elastic  resistance  for  higher  temperatures  becomes  very 
small . 

In  order  to  explain  the  strong  influence  of  the  stress  related 
properties  on  the  probability  of  fracture,  the  concept  of  fragmentation 
potential  was  introduced  by  Thirumalai  (14).  This  concept  is  based 
on  the  observation  that  the  thermal  load  vector  is  proportional  to  the 
factor 


and  thus,  in  the  absence  of  mechanical  loading,  the  resulting  stress 
field  depends  directly  on  F  (7l). 

In  the  nondimensional  form,  Thirumalai  defines  the  fragmentation 
potential  as 


F*  =  ^  (5.8) 

max 

where  F„,„  denotes  the  maximum  value  of  F. 
max 

Figure  16  shows  the  variation  of  F  with  temperature  for 
Dresser  basalt.  If  the  properties  were  assumed  to  be  constant,  the 
variation  in  F*  would  have  been  linear  as  shown  by  the  broken  line. 

The  difference  between  the  ordinates  for  the  solid  and  the  broken 
line  is  an  indication  of  the  influence  of  the  stress  field  of  the 
nonlinear,  temperature  dependent  properties  as  given  in  Table  7  • 
Another  method  of  predicting  whether  or  not  the  temperature  dependence 
of  the  thermoelastic  properties  will  have  a  significant  effect  on 
the  stress  field  is  through  the  use  of  what  will  be  referred  to  as 
the  nonlinearity  coefficient.  This  coefficient  will  be  def'  .jd  as 
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Figure  16.  Fragmentation  Potential  of  Dresser  Basalt  at  Elevated  Temperatures  (14) 


Figure  17.  Nonlinearity  Coefficient  of  Dresser  Basalt  at  Elevated  Temperatures 


76 


N 


* 


N 


max 


(5.9) 


where 


and  denotes  the  maximum  value  or  N. 

max 

The  usefulness  of  this  definition  lies  in  the  fact  that,  for  a 
homogeneous  solid  under  pure  thermal  loading  it  is  the  variation  of 
N*  which  influences  the  stress  field  rather  than  the  individual  variations 
of  each  of  the  thermoelastic  properties.  Thus,  for  materials  exhibiting 
thermoelastic  properties,  it  is  possible  to  obtain  quite  accurate  stress 
solutions  provided  the  nonlinearity  coefficient,  N  ,  remains  fairly 
constant. 

Variation  of  the  nonlinearity  coefficient  for  basalt  is  shown  in 
Figure  17.  The  highly  nonlinear  nature  of  this  curve  indicates  that 
for  a  reasonably  accurate  solution,  the  thermoelastic  properties  should 
be  allowed  to  vary  with  temperature  in  the  stress  analysis. 

3.  Fracture  Predictions 

Fracture  analysis  is  performed  based  on  the  Griffith  and  McClintock- 
Walsh  modified  Griffith  criteria.  Once  the  principal  stresses  have  been 
obtained  from  the  stress  analysis,  the  following  properties  are  required 
in  order  to  predict  the  initiation  of  fracture:  uniaxial  tensile  strength, 
0^,  uniaxial  compressive  strength,  o^,  and  the  fracture  surface  coef- 
ficient  of  friction,  of  the  material. 

Confining  pressure  is  known  to  increase  the  strength  of  the  rock. 

This  effect,  however,  is  neglected  here  since  this  investigation  is 
concerned  with  rocks  a  few  feet  below  the  surface  of  the  earth  (41 ). 
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Temperature,  on  the  other  hand,  has  a  weakening  effect  on  the 
rock  strength  (41).  This,  too.  Is  neglected  in  view  of  the  localized 
nature  of  heating  and  the  observation  that  the  major  portion  of  thermally 
fractured  rock  volume  experiences  only  a  small  average  increase  in 
temperature. 

Table  8  gives  the  properties  used  in  the  fracture  analysis. 

TABLE  8 

Propertiesof  Basalt  used  for  Fracture  Predictions  (39,  59,  63,  72) 

Uniaxial  tensile  strength,  2,195  psi 

Uniaxial  compressive  strength,  42,399  psi 

Fracture  surface  coefficient  of 
friction, 


0.9 


Chapter  VI 

MATHEMATICAL  FORMULATION  AND  METHOD  OF  ANALYSIS 
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The  three-dimensional  thermoelasticity  problem  and  Its  reduction 
to  three  two-dimensional  mathematical  models  were  considered  In  Chap¬ 
ter  VI.  In  Chapter  V,  hard  rock  characteristics  were  discussed  and 
simplifying  assumptions  Introduced  for  theoretical  studies.  In  this 
chapter,  analytical  aspects  such  as  mathematical  formulation  of  equa¬ 
tions  and  methods  of  their  solution  are  considered. 

The  thermal  fragmentation  analysis  considered  here  Involves  rela¬ 
tively  low  stress  levels  and  very  low  strain  rates  as  compared  to  con¬ 
ventional  explosive  and  nuclear  blasting  techniques.  Under  these  con¬ 
ditions,  the  effects  of  thermoelastic  coupling  are  negligible.  Thermal 
inertia  effects  are  also  neglected  as  the  time  rate  of  tetnperature 
change  Is  very  slow.  Thus,  the  general  problem  can  be  formulated  with¬ 
in  the  framework  of  the  linear,  uncoupled  theory  of  thermoelasticity. 

For  this  formulation,  the  problem  becomes  explicit,  that  solutions  can 
be  obtained  in  two  distinct  steps;  solution  to  a  well-defined  heat  con¬ 
duction  problem,  and  solution  of  the  stress  problem  with  known  tempera¬ 
ture  distribution. ( 68,  71). 

A.  Heat  Conduction  Problem 

The  geometries  for  the  slot  and  the  hole  models  considered  In  this 

investigation  are  shown  In  Figure  18a  and  b.  Parameters  used  to  describe 

these  geometries  are  as  follows: 

c  =  center-to-center  distance  between  the  hole  (or  slot)  axes, 
d  =  hole  diameter  (or  slot  width), 

A  =  total  hole  (or  slot)  depth, 
a  =  melt  deptli, 

L  =  fracture  propagation  length. 


X 

'  L/2 
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In  the  subsequent  discussion,  the  remaining  hole  or  slot  depth, 

A-a,  is  referred  to  as  convection  depth,  a^.  Parameter  B  is  chosen 
arbitrarily  but  large  enough  so  as  to  satisfy  both  the  traction  free 
and  ambient  temperature  conditions  on  boundary  y  =  -B  of  the  slot  model. 

In  the  absence  of  internal  heat  generation,  the  governing  equation 
of  two-dimensional  heat  conduction  for  a  homogeneous,  isotropic  solid 
with  constant  properties  is  given  by 


where 

k  =  thermal  conductivity 
p  =  density 

c  =  specific  heat,  and 

T  =  T(x,y,t)  denotes  the  temperature  at  point  (x,y)  at  time  t 

For  notational  convenience,  the  ambient  temperature  is  assumed 
to  be  zero  so  that  the  initial  condition  becomes 

T(x,y,0)  =  0  (6.2) 

The  boundary  conditions  for  the  slot  model  are: 

-  0  ,  y  <  0  (6.3) 

(  x  =  c/2 

x  *  d/2,  y  >  a  (6.4) 

y  =  A  (6.5) 

I  0  -  X  ^d/2.  y  =  0  (6.6) 

(  X  =  d/2,  0  ^  y  ^  a 

y  =  -B.  (6.7) 


—  =  0 
Dx  ^  ’ 


II 

3x 


+  hT  =  0, 


|I^hT  =  0. 


T  =  0. 
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The  boundary  conditions  for  the  hole  model  are: 


(  X  =  0.  y  <  d/2 
I  X  =  c/2 


(6.8) 


|I=0,  |x>d/2,y.O  (5  5) 

I  y  =  c/2 

T  =  T„,  x'.y=.dV4,  (6.:0) 

In  the  above  equations,  denotes  the  melt  temperature  of  the 
material  and  h  is  the  coefficient  of  convection  heat  transfer. 

Solution  Techniques 

A  great  deal  of  work  has  been  done  on  solving  problems  of  unsteady 
state  or  transient  heat  conduction  (69,  73-76).  Analytical  tools  include 
such  classical  techniques  as  separation  of  variables,  integram  transforms. 
Green's  functions  and  those  based  on  variational  principles.  The  integral 
methods  introduced  by  Goodman  (77)  and  the  Galerkin's  variational 
method  (78)  are  widely  used  to  obtain  approximate  closed  form  solutions. 
Each  method  has  some  advantages  over  the  others;  however,  each  has  its 
own  limitations  and  the  choice  of  the  method  to  solve  a  particular  prob¬ 
lem  depends  largely  on  the  nature  of  the  problem  itself. 

With  the  advent  of  high-speed,  large  scale  digital  computers, 
the  approximate  numerical  finite  difference  and  finite  element  methods 
have  b^icome  more  popular  due  to  their  ease  of  application  and  capability 
of  handling  today's  highly  nonlinear,  complex  problems. 

Even  though  the  temperature  problem  considered  in  this  investiga¬ 
tion  is  linear  in  that  constant  average  thermal  properties  are  used. 


the  boundary  conditions  are  quite  complex.  Also,  for  the  stress  analy¬ 
sis  part,  properties  are  allowed  to  vary  with  temperature  which  makes 
obtaining  a  closed  form  solution  extremely  difficult.  Thus,  closed 
form  temperature  solutions  will  be  of  little  advantage  and  since  a  two- 
dimensional  code  was  readily  available,  the  numerical  finite  element 
solution  was  chosen. 

Application  of  the  finite  element  method  for  the  solution  of  prob¬ 
lems  governed  by  a  general  quasi -harmonic  differential  equation  is 
discussed  by  Zienkiewicz( 79) .  The  well-known  Laplace  and  Poisson  equa¬ 
tions  are  particular  cases  of  the  general  quasi -harmonic  equation  and 
govern  such  frequently  encountered  problems  as  those  of  heat  conduction, 
seepage  flow,  distribution  of  electromagnetic  potential,  and  torsion 
and  bending  of  prismatic  members.  The  finite  element  conduction  code 
used  in  this  analysis  was  developed  based  on  the  formulation  given  by 
Wilson  and  Nickell  (89).  Details  of  this  code  are  proprietary  in  nature 
and  hence  are  not  given  here. 

B.  Stress  Problem 

The  problem  of  linear,  uncoupled  thermoelasticity  involves  fifteen 
equations  with  fifteen  unknowns;  six  stress  components,  six  strain 
components,  and  three  displacement  components.  The  field  equations  in 
index  notation  are  as  follows: 

Equilibrium  equations: 

",j_j  *  ■  0  {6.11) 

Stress-strain  relations: 

“^ij  ^ij^‘'kk  ^  j(3X+2u)aT  (6.12) 

Strain-displacement  relations: 


-  th(t  stress  tensor 
(  y  «  th*  strain  tensor 
0^  *  the  tlisplaceiMfnt  vector 
*  the  body  force  vector 
<S^j  “  the  Kronecker  delta 
It  ■  the  coefficient  of  thertnal  expansion »  and 
X,p  “  Lame's  constants. 


It  should  be  noted  that  the  stress  tensor, 

'U 


*  ^'jl 


J1 


1  f*  J. 


Also,  since  six  strain  components  are  expressed  In  tenns  of  only 
three  displacement  components,  the  following  contpatlbll Ity  relations 
must  be  satisfied; 


where 


Y  Y  • 

'pi'm’qi* 


rs,  mn 


10.14) 


^IJk  vlternaiing  tensor. 

Plane  Theory  of  Thermoelasticity 

For  simply  comiectoo  regions.  In  the  absence  of  body  forces,  the 
general  formulation  given  above  simplifies  for  the  plane  theory  of 
thermoelasticity  to  the  solution  of  the  biharmonic  equation 

V*V*F  +  cEV*(aT)  ■  0  (6.15) 


where 

V*  '  ~Y  •  l^he  Laplaclan 

3x*  3y* 

^  ^  1,  plane  stress 

y~  ,  plane  strain 

V  ■  Poisson's  ratio 

E  *  Young's  modulus,  and 

a  »  coefficient  of  thermal  expansion 


Once  the  Airy  stress  function  I-  is  detennined,  the  stresses  cen 


be  obtrtined  frotn  the  fol lowing  ruleiions: 


(6.16) 

•  .i*f 

(6,17) 

■‘x-ly 

(6.ia) 

It  is  important  to  note  that  the  solution  to  the  biharmonic  equa> 
tion  is  to  bo  obtained  subject  to  the  boundary  corditions.  Since 
both  the  slot  and  the  hole  models  considered  in  tnis  investigation 
involve  mixed  boundary  conditions,  that  is,  specif  led  tractions  as 
well  as  specified  displacements,  solution  of  the  bihannonic  equation 
in  itself  becomes  quite  complicated.  Further  complications  are  intro¬ 
duced  due  to  thP  t-'mpi,rature  v-.riatlons  of  ‘■ho  thenm'''last1r  proper¬ 
ties,  E,  V,  and  n.  No  attempt  was  therefore  made  to  obtain  closed 
form  solutions.  Rather,  approximate'  numerical  solutions  wore  sought 
u'llng  the  finite  element  method  which  due  to  its  generality  of  appli¬ 
cation  provides  a  powerful  tool  for  the  solution  of  complex  structural 
problems  with  arbitrary  geometries  and  nonlinear  material  properties. 

Numerous  finite  element  codes  have  been  developed  in  recent  years 
for  the  solution  of  elasto-plastic  and  thermal  stress  problems.  How¬ 
ever,  no  suitable  code  was  available  for  the  transient  thermoelastic 
stress  analysis.  The  available  codes  for  thermoelastic  analysis  are 
quasi -static  in  nature  and  can  handle  stationary  temperature  distribu¬ 
tions  only.  This  means  that  for  a  given  problem  geometry,  the  tempera¬ 
ture  distribution  at  each  time  step  must  be  run  separately  through  the 
stress  program  to  detennine  the  stress  distribution  as  a  function  of 


time.  This  procedure  becomes  highly  Inefficient  es  the  program  re- 
cofnputes  the  entire  stiffness  matrix  end  load  vector  for  each  of  the 
sets  of  temperature  data. 

A  now  finite  element  code  was»  therefore,  developed  for  the  tran¬ 
sient  thermal  stress  analysis.  This  code  is  Identified  by  the  acronym 
"TRATSA"  (TRAnslent  Thermal  Stress  Analysis)  and  can  be  used  for  plane 
or  axisyiinietric  bodies  with  temperature-dependent  material  properties. 

The  code  Is  based  on  theoretical  formulations  given  by  Zienklewicz  (79) 
and  Jones  and  Crose  (81).  These  formulations  assume  linear  displace¬ 
ment  between  nodes,  resulting  in  constant  stress  elements.  Both  tri¬ 
angular  and  quadrilateral  elements  can  be  used.  The  Input  Instructions 
and  the  program  listing  are  given  In  Reference  02. 

C.  fracture  Analysis 

The  fracture  predictions  made  In  this  Investigation  are  exploratory 
1n  nature  and  are  based  on  the  str'ess  field  which  Is  obtained  by  ojsum- 
Ing  the  material  to  remain  a  continuum.  In  practice,  however.  Initia¬ 
tion  of  a  crack  will  change  the  temperature  field  which.  In  turn,  effects 
the  stresses.  Cracks  also  tend  to  relieve  the  stresses,  but  at  the  same 
time,  they  also  act  as  stress  raisers.  A  rigorous  treatment  ot  the 
problem  would  require  progressive  Introduction  of  additional  convection 
and  traction  free  boundary  conditions  with  the  creation  of  new  fracture 
surfaces. 

If  the  propagation  of  a  crack  Is  to  be  followed,  the  numerical 
treatment  of  the  thermoelasticity  problem  becomes  coupled  In  the  sense 
that  once  the  crack  Initiation  Is  predicted  from  the  known  stress 
field,  one  must  go  back  and  recompute  the  temperature  distribution  for 
the  next  time  step  allowing  for  the  newly  developed  fracture  surfaces. 
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For  an  extensive  investigation  such  an  approach  becomes  Impracticable 
due  to  the  amount  of  work  Involved  In  the  reformulation  of  the  prob¬ 
lem  at  each  time  step  and  the  resulting  p>*oh1b1tively  large  amount 
of  machine  time.  Apart  from  this,  for  problems  not  Involving  dynamic 
fracture,  there  Is  little  Justification  In  using  a  progressive  fracture 
approach  as  the  crack  propagation  of  nondynamic  nature,  In  actual  prac¬ 
tice,  Is  largely  governed  by  the  orientation  of  the  pre-existing  micro¬ 
cracks  and  the  Interactions  between  them;  factors  of  which  no  theoreti¬ 
cal  knowledge  Is  available. 

The  fracture  predictions  In  this  analysis  are  based  on  the  Griffith 
and  the  McCIIntock-Walsh  modified  Griffith  criteria.  These  fracture 
theories  are  themselves  approximate  as  they  are  derived  from  an  energy 
formulation  and  neglect  the  effects  of  stress  concentrations  and  the 
Interactions  between  the  cracks.  These  theories,  therefore,  can  at 
best  be  considerel  statistical  In  nature.  This  Is  evidenced  by  the 
results  of  Laurlello  (68)  and  Bienlawski  (60)  who  obtained  "statisti¬ 
cally"  good  correlations  between  the  experimentally  observed  fracture 
fields  and  those  predicted  theoretically  using  the  Griffith  and  the 
modified  Griffith  fracture  criteria. 

Fracture  Propagation  Studies 

During  the  preliminary  analysis,  1t  was  found  that,  predicting 
the  Initiation  of  a  crack  and  following  its  propagation  even  under  the 
assumption  that  the  presence  of  a  crack  does  not  effect  the  temperature 
and  stress  fields,  requires  a  rather  large  number  of  stress  runs  with 
very  small  time  steps.  Stress  runs  with  large  time  steps,  on  the 
other  hand,  resuH  in  large  fracture  zones  giving  no  indication  what¬ 
soever  of  the  actual  crack  propagation  path.  This  observation 
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necessitated  some  sort  of  approximate  procedure  that  could  serve 
as  a  compromise  between  the  two  extremes  and  would  yield  practically 
comprehensible  propagation  data. 

A  procedure  was  devised  based  on  a  concept  referred  to  hereafter 
as  the  concept  of  fracture  intensity  level.  With  this  concept,  the 
crack  can  be  predicted  to  lie  within  a  small  narrow  band  rather  than 
a  whole  large  fractured  zone  as  obtained  from  the  stresses  using 
large  time  steps.  The  fracture  intensity  level  is  defined  in  terms 
of  the  stress  magnitude  in  excess  of  that  necessary  for  fracture. 

Thus,  a  fractured  element.  A,  having  a  higher  excess  stress  than  some 
other  fractured  element,  B,  will  be  identified  by  a  fracture  intensity 
level  higher  than  that  of  element  B. 

This  procedure,  of  course,  does  not  eliminate  the  trial  and  error 
method  necessary  in  finding  the  solution  associated  with  the  complete 
fracture  of  the  solid;  that  is,  solution  when  fracture  reaches  a  free 
surface.  Once  this  solution  is  obtained,  excess  stresses  associated 
with  fractured  elements  are  computed.  The  total  range  from  minimum 
to  maximum  excess  stress  magnitude  is  divided  into  a  number  of  inter¬ 
vals  and  labeled  in  terms  of  the  fracture  intensity  levels.  The  com¬ 
puter  code  is  written  to  plot  the  fractured  elements  together  with  their 
respective  stress  intensity  levels.  From  the  plotted  output,  it  is 
possible  to  trace  the  approximate  crack  propagation  path;  with  the 
crack  initiating  at  the  point  of  the  highest  stress  intensity  level 
and  progressing  in  the  direction  of  decreasing  fracture  levels. 
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D.  Summary 

The  thermal  fragmentation  analysis  is  perfonned  in  three  steps. 

The  first  step  involves  the  temperature  solutions  obtained  by  using 
a  finite  element  conduction  code.  Average  thermal  properties  are 
used  in  this  part  of  the  analysis.  The  temperature  results  are  then 
used  as  input  into  the  finite  element  stress  code  TRATSA.  Thermo¬ 
elastic  properties  are  allowed  to  vary  with  temperature  for  the  stress 
analysis.  Principal  stresses  are  computed  and  read  through  the  frac¬ 
ture  code.  Fracture  predictions  are  based  on  the  Griffith  and  the 
McClintock-Walsh  modified  Griffith  criteria.  The  solid  is  assumed 
to  remain  a  continuum  and  the  effects  of  crack  initiation  on  tempera¬ 
ture  and  stress  fields  are  neglected.  The  concept  of  fracture  inten¬ 
sity  level,  based  on  the  stress  magnitude  in  excess  of  that  required 
for  fracture  was  introduced  to  obtain  an  approximation  to  the  frac¬ 
ture  propagation  path. 
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Chapter  VII 
SLOT  MODEL  ANALYSES 

As  mentioned  in  Chapter  IV,  the  slot  model  is  obtained  by  passing 
a  cutting  plane  along  a  series  of  holes  and  observing  the  geometry 
that  is  projected  on  the  cutting  plane.  This  model  was  designed  for 
the  parametric  study  of  the  fracture  which  occurs  on  a  plane  perpendi¬ 
cular  to  the  hole  axes  passing  approximately  through  the  thermal  inclu¬ 
sions.  This  fracture  is  parallel  to  the  working  face  and  makes  possible 
the  removal  of  a  layer  equal  to  the  depth  of  the  thermal  inclusions. 

In  the  subsequent  discussion  this  subsurface  fracture,  due  to  its  orienta 
tion  with  respect  to  the  working  face,  is  referred  to  as  the  parallel 
fracture.  The  significance  of  the  process  parameters  which  are  the 
hole  diameter,  hole  depth,  hole  spacing  and  the  depth  of  heating  on 
the  actual  removal  of  the  layer  can  be  approximated  by  studying  their 
influence  on  the  parallel  fracture.  The  slot  model  studies,  therefore, 
are  of  primary  importance. 

In  order  to  study  the  significant  process  parameters,  twenty  differ¬ 
ent  slot  models  were  analyzed.  Effects  of  hole  diameter  and  hole  spac¬ 
ing  were  studied  by  using  three  different  values  of  each;  the  combina¬ 
tions  giving  nine  different  fracture  lengths.  Effects  of  convection 
depth  and  melt  depth  were  investigated  using  various  combinations  of 
these  values  for  a  given  hole  diameter  and  hole  spacing.  Parametric 
description  of  these  models  is  given  in  Table  9  . 

A.  Grid  Size  Effects 

Accuracy  of  the  finite  element  solution  depends  largely  on  the 
size  of  the  elements;  the  smaller  the  element  size,  the  more  accurate 
the  solution.  However,  computer  time  increases  with  an  increase  in 
the  number  of  elements.  Hence,  it  is  necessary  to  adjust  the  finite 
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TABLE  9 


Parametric  Description  of  Slot  Models 


Model 

Hole  Diam 
d  ,in. 

Spacing 
c ,  i  n . 

Fracture  Length 

L  =  c-d 

Melt  Depth 
a, in. 

Conv.  1 

a  =  i 
c 

lA 

1.0 

4.0 

3.0 

0.5 

1.5 

lA-1 

2.0 

0.0 

1A-2 

2.0 

2.0 

lA-3 

4.0 

0.0 

IB 

1.0 

8.0 

7.0 

0.5 

3.0 

1C 

1.0 

12.0 

11.0 

0.5 

5.0 

2A 

2.0 

4.0 

2.0 

0.5 

1.0 

2A-1 

1.2 

0.8 

2A-2 

0.5 

0.1 

2B 

2.0 

8.0 

6.0 

0.5 

3.0 

2B-1 

0.1 

2C 

2.0 

12.0 

10.0 

0.5 

5.0 

2C-1 

3.0 

2C-2 

0.1 

3A 

1.5 

4.0 

2.5 

0.5 

1.5 

3B 

8.0 

6.5 

3.0 

3C 

12.0 

10.5 

5.0 

llA 

2.0 

5.0 

3.0 

0.5 

1.5 

13B 

8.5 

6.5 

3.0 

13C 

12.5 

10.5 

5.0 
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element  grid  so  that  acceptable  solutions  can  be  obtained  economically. 
Fairly  large  elements  can  be  used  in  the  areas  where  low  stresses  are 
expected,  whereas  regions  of  high  stress  concentrations  should  be 
divided  into  very  fine  grids.  Thus,  the  problem  geometry  and  the  boun¬ 
dary  conditions  should  be  given  primary  consideration  in  drawing  the 
finite  element  grid.  For  thermal  stress  problems,  the  magnitude  of 
the  temperature  gradient  must  also  be  taken  into  consideration  as  it 
greatly  influences  the  stress  field. 

For  the  slot  model  analysis,  effects  of  element  size  and  time  in¬ 
crement  were  investigated.  A  value  of  0.2  seconds  for  the  time  incre¬ 
ment  for  a  typical  grid  shown  in  Figure  19  was  found  to  yield  quite 
accurate  results.  Size  of  the  smallest  elements  used  in  this  grid 
is  0.025  in.  square. 

B.  Melt  Depth  Studies 

The  stress  State  in  the  elastic  rock  surrounding  the  heater  and 
molten  rock  inclusion  depends  on  the  temperature  field  and  mechanical 
constraints  in  the  form  of  specified  displacement  boundary  conditions. 
Due  to  the  poor  thermal  conductivity  of  rocks,  the  resulting  thermal 
gradients  are  very  steep.  This  requires  very  small  sized  elements 
for  a  reasonably  accurate  analysis.  It  is  necessary,  therefore,  to 
reduce  the  semi -infinite  problem  geometry  to  a  finite  region  for  numeri¬ 
cal  treatment. 

The  experimentally  observed  fracture'-  occur  in  the  vicinity  of 
the  slot  base.  This  location  for  the  secondary  fractures  is  highly 
desirable  from  an  economical  viewpoint.  Hole  depths  used  in  the  field 
tests  varied  from  12.0  in.  to  about  30.0  in.  and  the  thermal  inclusions 
were  observed  to  be  in  the  range  from  4.0  in.  to  12.0  in.  in  length. 
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These  values  were  much  too  high  to  be  used  for  the  numerical  finite 
element  analysis.  However,  for  a  given  geometry  under  identical  mechan¬ 
ical  constraints,  the  stress  field  in  an  area  depends  only  on  the  tem¬ 
perature  field  in  that  area.  Thus,  it  is  possible  to  obtain  a  reason¬ 
ably  accurate  stress  field  in  the  vicinity  of  the  slot  base  by  using 
small  melt  depths  provided  the  temperature  field  in  this  area  remains 
fairly  constant  for  different  melt  depths. 

Effects  of  different  melt  depths  and  convection  depths  on  the  tem¬ 
perature  field  in  the  vicinity  ot  the  slot  base  were  investigated  using 
models  lA-1 ,  lA-2,  and  lA-3.  Parametric  description  of  these  models 
is  given  in  Table  9. 

Results  of  these  studies  are  plotted  in  Figures  20-25  which  show 
the  temperature  plots  for  two  different  time  values.  Temperature  con¬ 
tours  are  plotted  in  increments  of  100°C  except  for  the  50®C  contour. 

The  melt  temperature  useJ  in  the  analysis  was  1250"C.  Comparison  of 
these  results  show  that  the  temperature  field  in  the  vicinity  of  the 
slot  base  does  remain  fairly  constant  for  different  combinations  of 
melt  depth  and  convection  depth  values.  However,  it  should  be  realized 
that  for  the  above  statement  to  be  true,  there  exists  a  lower  limit  on 
the  value  of  the  melt  depth,  a.  A  value  of  0.5  in.  was  found  to  be 
sufficiently  high  in  order  to  obtain  an  invariant  temperature  field  in 
the  immediate  surrounding  of  the  slot  base.  The  subsequent  analysis 
was,  thetefore,  completed  using  a  melt  depth  valiie  of  0.5  in. 

C .  Convection  Depth  Studies 

As  seen  from  the  results  shown  in  Figures  20-25  ,  the  temperature 
distribution  in  the  vicinity  of  the  slot  base  is  independent  of  the 
convection  depth.  However,  for  small  melt  depths,  the  stress  field 


Figure  22.  Temperature  Distribution,  Model  lA-3,  t  ■  60  sec 


Figure  24.  Temperature  Distribution,  Model  1A-2,  t  ■  120  sec. 
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In  this  rpqlon  i;t  highly  dependent  on  the  convection  depth.  This 
effect  results  from  two  factors.  Although  the  tetnperature  distribu¬ 
tion  as  seen  from  Figures  20-25  Is  approximately  one-dimensional  In 
the  melt-depth  section,  thermal  gradients  near  the  point  of  transi¬ 
tion  from  melt  condition  to  convection  condition  are  highly  localized. 
These  highly  localized  gradients  give  rise  to  very  high  stress  con¬ 
centrations.  Also,  due  to  the  localized  heating,  the  thermal  expan¬ 
sion  of  this  area  Is  restrained  by  the  relatively  large  cold  zone. 

This  thermally  Induced  constraint  induces  very  high  compressive  stresses 
In  the  heated  zone  and  tensile  stresses  in  tHa  cold  zone.  For  a  given 
temperature  distribution,  the  severeness  of  the  thermal  constraint 
Increases  with  Increase  in  the  volume  of  the  cold  zone,  up  to  a  cer¬ 
tain  critical  value.  Once  the  cold  zone  volume  reaches  this  critical 
value,  any  further  Increase  has  no  effect  on  the  stress  field.  For 
the  slot  models,  since  the  volume  of  the  cold  zone  depends  directly 
on  the  convection  depth,  there  exists  a  critical  value  of  the  convec¬ 
tion  depth  beyond  which  any  increase  In  the  convection  depth  has  no 
significant  effect  on  the  fracture  inducing  stresses. 

The  convection  depth  studies  were  conducted  using  models  2A,  2A-1, 
2A-2;  2B,  2B-1 ;  and  2C,  2C-1 ,  2C-2.  The  hole  diameter  for  these  models 
was  kept  constant.  Three  values  of  hole  spacing  were  considered,  giv¬ 
ing  three  different  fracture  lengths.  For  models  of  equal  fracture 
length,  convection  depths  were  varied,  the  highest  value  being  approxi¬ 
mately  equal  to  half  the  fracture  length.  Parametric  description  of 
these  models  is  given  in  Table  9. 

The  results  of  these  studies  are  shown  in  Figures  26-33.  In  view 
of  the  assumption  that  the  material  above  700°C  has  no  elastic  resis¬ 
tance,  the  problem  geometry  is  modified  using  the  700°C  isotherm. 
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Figure  27.  Fracture  Zone,  Model  2A-1 
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Figure  31.  Fracture  Zones,  Model  2C 
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Figure  33.  Fracture  Zones,  Model  2C-2 
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The  fracture  completion  time  Is  denoted  by  t^f,  and  thu  probable  frac¬ 
ture  zones  at  this  time  are  shown  by  the  cross -hatched  area.  The  doubly 
cross-hatched  area  shows  the  fracture  zone  governed  by  the  McClIntock- 
Walsh  modified  Griffith  criterion,  Equation  2.4,  The  remainder  of  the 
fracture  zone  Is  governed  by  the  original  Griffith  criterion,  Equations 
2.  land  2.2.  In  the  subsequent  discussion,  the  two  types  of  fracture 
zones  are  referred  to  as  the  McClIntock-Walsh  zone  and  the  Griffith 
zone,  respectively. 

The  McCl Intock-Walsh  zone  Is  associated  with  very  high  compressive 
stress  components,  l-or  very  small  or  no  convection  depths,  the  fracture 
Is  characterized  by  a  McCl Intock-Walsh  zone  1n  the  high  temperature 
region  followed  by  a  small  Griffith  tension  zone  In  the  low  temperature 
area  and  a  large  secondary  McClIntock-Walsh  zone  extending  to  the  line 
of  synriietry.  The  Increase  In  the  cold  zone  volume  and  the  consequent 
thermal  constraint  as  a  result  of  the  Increase  1n  the  convection  depth 
has  relatively  little  influence  on  the  McCl Intock-Walsh  zone  In  the 
high  temperature  area  and  the  small  Griffith  zone  Immediately  follow¬ 
ing  it.  This  Is  due  to  the  fact  that  the  stress  state  In  these  zones 
result  primarily  from  the  temperature  gradients  and  that  changes  In 
the  convection  depth  have  very  little  effect  on  the  temperature  field 
In  this  region. 

The  McCl Intock-Walsh  zone  in  the  cold  region  Is,  however,  greatly 
influenced  by  convection  depth  changes.  As  the  convection  depth  Is  In¬ 
creased,  the  thermal  constraint  becomes  more  severe.  As  a  result, 
the  high  compressive  stress  components  in  the  cold  region  start  to  de¬ 
crease  in  magnitude  and  the  tensile  stress  components  increase  sharply. 
The  net  effect  of  these  changes  in  the  stress  components  is  that  the 
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ittcondAry  HtiCHntock-W«nh  ion«  begins  to  ihrlnk,  If  th<>  uonvttctlon 
dtpth  U  contimiuuily  a  It  rt»Ach»d  when  the  tecondAry 

McCl Intuck'Walth  tone  O' iHettly  d(k«|)|)eA»*s  «nd  trAnifomk  into  the 
Griffith  tenilon  tone  which  now  oxtendt  to  the  line  of  lyonwtry.  In 
thit  ditcuitlont  the  convection  depth  Attocloled  with  the  completion 
of  the  ibove  trentfomwtlon  It  »'eferrod  to  et  the  critical  convection 
depth. 

The  critical  convection  depth  wat  found  to  be  approtlniatnly  equal 
to  half  the  fracture  length,  From  the  retultt  for*  convection  depths 
approximately  equal  to  the  critical  value,  It  It  teen  that  the  total 
fracturo  zone  consist:*  of  a  small  McCl Intock-Walsh  zone  In  the  high 
temperature  region  followed  by  a  Griffith  tension  zone  In  the  low  tem¬ 
perature  region  and  extending  all  the  way  to  the  line  of  symmetry. 

For  convection  depths  smaller  than  the  critical  value,  the  fractu>*e 
time  shows  only  small  variations  with  changes  In  convection  depth, 
However,  as  the  convection  depth  approaches  the  critical  value,  the 
McClintock-Walsh  ione  rapidly  transfoms  into  the  Griffith  tension 
zone,  and  because  of  the  tensile  strength  of  the  rock  being  much  lower 
than  Its  compressive  strength,  the  fracture  time  is  drastically  reduced. 
Once  the  fracture  mode  transfommtion  is  completed,  the  fracture  time 
shows  a  negligibly  small  decrease  with  any  further  ir:rease  in  the  con¬ 
vection  depth. 

The  fracture  zones  in  Figures  26-33  show  another  very  important 
effect  of  the  convection  depth  changes.  This  effect  involves  the  loca¬ 
tion  of  the  parallel  fracture.  The  optimum  location  of  the  parallel 
fracture  requires  the  crack  initiation  to  occur  very  close  to  the 
hole  base  and  the  crack  propagation  in  a  plane  approximately  parallel 
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In  oNor  to  obtain  th«  opltmunt  location »  tht  incluilon  thoulti 

b«  eonctntr«t(»cl  «t  tht  vory  b«ii  of  th»  hoU  while  the  convection  depth 
iHouUli  ot  leeit.  be  equil  to  helf  the  frecture  lenqth,  4nd  thet  eny 
further  Increete  In  the  convection  depth  hes  only  a  trlvUl  Influence 
on  the  locAtlon  of  the  perellel  crecki  end  thp  frecture  time. 

0 .  Hale.  JUiinf^r  JiMiPAQiOlL 

Hole  diemeter  end  the  hole  specinu  ei'e.  perheps.  the  two  most 
significent  parameters  In  the  study  of  rock  fragmentation  using  subsur¬ 
face  thennal  Inclusions.  The  hole  diameter  affects  the  fracture  length 
as  well  as  the  site  of  the  thermal  Inclusion,  The  hole  spacing »  how¬ 
ever,  Is  the  more  Important  factor  as  It  directly  controls  the  volunwj 
of  the  rock  removed, 

In  order  to  Investigate  the  effects  of  the  hole  diameter  and  the 
hole  spacing,  nine  fracture  lengths  were  considered,  Hole  spacings 
of  4.0  In.,  8.0  In.,  and  12.0  In.  were  used  and  for  each  hole  spacing, 
hole  diameter  values  of  1,0  In.,  1.5  In.,  and  2.0  in.  were  considered. 

In  Table  9,  the  models  used  for  the  hole  diameter  and  the  hole 
spacing  studies  are  Identified  as  1A,  IB,  1C;  2A,  2B,  2C;  and  3A,  3B, 

3C.  The  melt  depth  and  the  convection  depth  values  were  selected 
based  on  the  results  of  the  melt  depth  and  the  convection  depth  studies 
described  in  the  beginning  of  this  chapter.  The  melt  depth  value 
was  kept  constant  at  0.6  In.  for  all  models.  For  a  given  model,  the 
total  hole  depth  considered  was  such  that  the  convection  depth  was 
approximately  equal  to  the  critical  value,  that  Is,  approximately 
half  the  fracture  length. 
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The  tefl(per«tu«'e  dlttrlbutioh  At  the  tinm  of  frectut'e  contpletlon 
weti  found  to  be  fer  from  «tettonei*y  end  hiphty  locellied  in  the  vicini¬ 
ty  of  the  itot  beie,  Plott  of  the  typicAl  tewpereture  field*  «t  the 
time  of  frecture  completion  ere  thoMn  in  fiQuret  34  -  36 ,  As  teen  from 
these  plots,  the  major  portion  of  the  fractured  volume  does  not  ex¬ 
perience  any  tee^peroture  change.  Thus,  the  stresses  in  this  region 
occur  entirely  due  to  the  load  vector  resulting  fro»n  the  themal  con¬ 
straint.  The  natural  thermal  expansion  of  the  high  temperature  region 
is  resisted  by  the  cold  volume.  As  a  consequence  of  this  constraint, 
the  hot  region  undergoes  a  contpressive  loading  while  the  cold  zone  is 
under  a  tensile  load  condition.  In  view  of  the  equilibrium  of  forces, 
the  ccjtnpressive  load  vector  integrated  over  the  entire  high  tempera¬ 
ture  region  must  equal  the  tensilo  load  vector  integrated  over  the  cold 
zone  volume.  Since  the  rock  fracture  is  attributed  to  the  tensile 
failure  because  of  the  tensile  strength  of  rock  being  much  lower  than 
its  compressive  strength,  only  a  small  hot  region  is  required  to  in¬ 
duce  fracture  causing  stresses  in  a  comparatively  large  cold  zone. 

As  the  hole  spacing  is  increased,  the  cold  zone  volume  also  in¬ 
creases.  A  higher  tensile  load  vector  is  therefore  necessary  in  order 
to  cause  the  fracture  Inducing  stresses  in  the  cold  zone.  This,  in 
turn,  requires  a  larger  natural  expansion  of  the  hot  region  and,  there¬ 
fore,  a  larger  volume  of  the  hot  zone  and  a  longer  time  for  fracture 
completion. 

Thus,  as  the  hole  spacing  is  increased,  the  fracture  completion 
time  also  increases.  However,  as  seen  from  Figures  34 -36 ,  the  increase 
in  the  fracture  completion  time  is  much  higher  than  the  corresponding 


Q 

ro 


25  •€ 
50  ®C 
I00*C 
300  ®C 
500  ®C 
700  ®C 


t  *0.97  sec. 


Figure  34,  Temperature  Distribution, 
Model  lA 
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Figure  36.  Temperature  Distribution,  Model  3C 
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increase  In  the  hot  zone  volume.  This  is  explained  by  the  fact  that 
while  the  thermal  expansion  and  the  resulting  force  vector  are  direct¬ 
ly  dependent  on  the  volume  of  the  hot  zone,  the  volume  of  the  hot  zone 
itself  is  a  complex  function  of  the  thermal  conductivity  and  the  ther¬ 
mal  diffusivity,  and  for  rocks  the  valuej  of  both  these  properties  are 
extremely  low. 

2.  Fracture  Length  Effects 

From  the  above  discussion,  it  is  clear  that  while  for  a  given 
hole  diameter,  increasing  the  hole  spaci  icreases  the  fracture  com¬ 

pletion  time,  it  must  also  be  true  that  for  a  given  hole  spacing,  the 
fracture  completion  time  will  decrease  if  the  hole  diameter  is  increased. 
In  other  words,  it  is  the  change  in  the  fracture  length  that  affects 
the  fracture  completion  time.  This  is  demonstrated  by  the  fracture  zone 
plots  shown  in  Figures  26,  29,  31,  37  -  39,  and  the  fracture  completion 
time  versus  the  hole  spacing  plots  shown  in  Figure  40. 

It  is  possible,  therefore,  to  study  the  effect  of  the  change  in 
hole  diameter  on  the  fracture  time  by  keeping  the  hole  diameter  un¬ 
changed  and  changing  the  hole  spacing  such  that  the  fracture  length 
will  remain  the  same.  Models  llA,  13B,  and  13C  as  described  in  Table  9 
were  designed  in  support  of  this  hypothesis.  These  models  were  ob¬ 
tained  by  increasing  the  diameters  of  models  lA,  3B,  and  3C  to  2.0  in. 
and  increasing  the  hole  spacing  correspondingly  so  as  to  keep  the 
fracture  lengths  unchanged. 

Plots  of  the  fracture  zones  for  models  llA,  13B,  and  13C  are 
shown  in  Figures  41-43.  Comparison  of  fracture  times  and  the  fracture 
plots  of  these  models  with  those  of  models  lA,  3B,  and  3C  shown  in 
Figures  44,  38,  and  39  confirms  the  above  mentioned  hypothesis. 
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Figure  37.  Fracture  Zones,  Model  3A 
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Figure  38.  Fracture  Zones,  Model  3B 
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Figure  39.  Fracture  Zones,  Model  3C 
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Figure  Al.  Fracture  Zones.  Model  13A 
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Figure  42.  Fracture  Zones,  Model  13B 
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rtoplh  *f«  lUtod  in  TohU  10.  Voluos  frow  thin  tobio  wer«  utod  to 
computf  the  frocturo  lenqth  rotioi  *nd  the  corrojpondlnQ  fracture  com¬ 
pletion  time  ratios  as  uiwen  in  Table  ll.  These  values  were  plotted 
on  a  loflarlthmlc  scale  as  shown  In  riqure  as  from  which  the  following 
relation  was  obtained; 

‘fl“’  (7.1) 

where  is  the  fracture  completion  time  ratio,  and  L*  Is  the  frac¬ 
ture  lenyth  ratio. 

The  usefulness  of  Equation  (7.1)11t7s  In  the  fact  that  the  frac¬ 
ture  completion  time  for  any  given  fracture  length  can  be  predicted 
fron\  the  for  a  single  test. 
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TABir  10 


frdctur©  Tit«*s  'Tor  Slot  Moduli  with  rt.'-L/? 


Mo<i«i 

Hol«  0i4in 
d,1n, 

Spacing 

c,in. 

Fracture  Length 
1.  ■  c-d 

Conv.  Depth 
in. 

Fracture  Time 
tf,  sec. 

lA 

1.0 

4.0 

3.0 

1.5 

0.97 

IB 

fl.O 

7.0 

3.0 

5.82 

li; 

12.0 

11.0 

5.0 

40.6 

2A 

2.0 

4.0 

2.0 

1.0 

0.43 

2R 

8.0 

6.0 

3.0 

4.4 

2C 

12.0 

10.0 

5.0 

25.4 

3A 

1.6 

4.0 

2.5 

1.5 

0.B5 

IR 

8.0 

6.5 

3.0 

5.0 

3C 

12.0 

10.5 

5.0 

34.0 

IIA 

r.o 

5.0 

3.0 

1.5 

0.98 

13B 

8.5 

b.5 

3.0 

5.0 

13C 

12.5 

10.5 

5.0 

34.0 
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TABLF  11 

Fracture  Length  and  Fracture  Time  Ratios  for 
Slot  Models  with  Convection  Depths  Approximately 
Equal  to  Half  the  Fr'^cture  Length 


Length 

L, 

Length  Ratio 
L*  -  ^ 

Time  Ratio 
>  -ta.. 

V  t, 

Length 

Li 

Length  Ratio 

Time  Ratio 

tt  “  ^ 

V  t, 

2.0 

5.5 

96.0 

6.0 

1.83 

9.23 

5.25 

80.2 

1.7S 

7.23 

5.0 

60.0 

1.67 

5.78 

3.5 

(3.73 

1.16 

1.322 

3.25 

11.8 

1.08 

1.136 

3.0 

10.38 

1.5 

2.28 

6.5 

1.69 

8.12 

1.25 

2.02 

1.61 

6.8 

1.53 

5.09 

2.5 

4.4 

47.4 

1.07 

1.165 

4.2 

39.7 

4.0 

29.7 

7.0 

1.57 

6.97 

2.8 

6.78 

1.5 

5.84 

2.6 

5.84 

1.42 

4.38 

2.4 

5.14 

1.2 

1.13 

10.0 

1.1 

1.59 

1.05 

1.34 

3.0 

3.67 

42.0 

3.5 

35.2 

10.5 

1.04 

1.19 

3.33 

26.3 

2.33 

6.02 

2.17 

5.17 

2.0 

4.55 

FRACTURE  LENGTH  RATIO,  I^^L,. 


Figure  45.  Dimensionless  Fracture  Time-Length  Rel 
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Chapter  VIII 
HOLE  MODEL  ANALYSIS 

The  hole  model  was  designed  to  study  the  influence  of  the  hole 
diameter  and  the  hole  spacing  parameters  on  the  fracture  which  occurs 
along  the  line  of  a  series  of  holes.  Because  of  its  orientation  with 
respect  to  the  working  face,  this  fracture  in  the  subsequent  dis¬ 
cussion  is  referred  to  as  the  perpendicular  fracture. 

The  model  is  obtained  by  passing  a  cutting  plane,  parallel  to 
the  working  face,  through  the  center  of  the  thermal  inclusions.  The 
geometry  of  this  model  is  therefore  that  of  a  thin  uniform  plate  with 
a  series  of  holes.  For  the  plane,  two-dimensional  thermal  stress 
analysis,  the  variations  in  the  temperature  and  the  stresses  across 
the  thickness  are  neglected.  For  the  actual  problem  however,  this 
is  only  true  for  the  melt  depth  section  as  the  hole  is  only  partially 
heated,  and  as  shown  in  the  previous  chapter,  the  fracture  initiates 
in  the  vicinity  of  the  point  of  transition  from  the  melt  condition  to 
the  convection  condition  and  propagates  to  the  surface.  It  was  also 
shown  that  the  fracture  times  for  the  parallel  fractures  using  small 
or  zero  convection  depths  were  much  higher  than  those  obtained  by 
using  convection  depths  approximately  equal  to  or  greater  than  half 
the  fracture  lengths.  Thus,  the  results  of  the  hole  model  analysis 
given  in  this  chapter  are  representative  of  the  upper  bound  solution 
of  the  three-dimensional  problem  with  small  hole  depths  heated  along 
their  entire  length. 

The  hole  model  studies  were  conducted  using  three  values  each  of 
the  hole  diameter  and  the  hole  spacing;  the  combinations  giving  nine 
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different  fracture  lengths  as  described  in  Table  12  .  The  hole 
diameter  and  the  hole  spacing  values  for  these  models  are  the  same 
as  those  for  the  parallel  fracture  studies  using  the  slot  models. 

A  typical  finite  element  grid  used  for  the  hole  model  analysis  is 
shown  in  Figure  46. 


TABLE  12 

Parametric  Description  of  Hole  Models 


Model 

Hole  Dia. 

Hole  Spacing 

Fracture  Length 

d,  in. 

c,  in. 

4.0 

L  =  c  -  d 

21A 

1.0 

3.0 

21B 

8.0 

7.0 

21C 

12.0 

11.0 

22A 

2.0 

4.0 

2.0 

22B 

8.0 

6.0 

22C 

12.0 

10.0 

23A 

1.5 

4.0 

2.5 

23B 

8.0 

6.5 

23C 

12.0 

10.5 

Temperature  Analysis: 


As  in  the  case  of  the  slot  model  studies,  the  fracture  time  values 


for  the  hole  models  were  observed  to  be  small  compared  to  those 
associated  with  the  steady-state  temperature  distribution.  As  a 
result,  the  temperature  field  at  the  time  of  the  fracture  is  highly 
localized  in  teh  vicinity  of  the  melt  condition  and  is  independent  of 
the  changes  in  the  boundary  conditions  at  the  far  end.  Thus,  the  hole 
models  are  characterized  by  almost  radial,  one-dimensional  temperature 
fields  as  typically  illustrated  In  Figures  47-49. 


32 


Figure  47.  »amperature  Distribution,  Model  22A 
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Figure  49.  Temperature  Distribution,  Model  22C 
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B.  Stress  and  Fracture  Analysis: 

Unlike  the  slot  models,  where  the  ratio  of  the  cold  zone  volume 
to  the  hot  zone  volume  varies  with  the  convection  depth,  the  volume 
ratio  for  a  given  hole  model  is  constant.  Thus,  the  fracture  plots 
for  the  hole  models  are  characteristic  of  those  for  the  slot  models  with 
deep  holes  heated  along  their  entire  depth.  For  slot  models  with 
very  small  convection  depths,  it  was  shown  in  Chapter  VII  that  the  frac¬ 
ture  zone  consists  of  a  McClintock-Walsh  zone  in  the  high  temperature 
region,  followed  by  a  small  Griffith  tension  zone  and  a  large  sec¬ 
ondary  McClintock-Walsh  zone.  Also,  the  fracture  completion  times 
for  small  convection  depth  models  were  found  to  be  much  higher  than  those 
given  in  Table  10  where  convection  depths  used  are  approximately 
equal  to  the  critical  values.  These  two  observations  are  also  re¬ 
flected  through  the  hole  model  analysis  as  seen  from  the  typical 
fracture  plots  shown  in  Figures  50-52. 

From  these  results,  it  can  be  Inferred  that  the  fractures  orig¬ 
inate  in  the  melt  depth  section  across  the  holes  and  propagate  to 
the  surface.  It  can  also  be  observed  that,  like  the  slot  models, 
the  fracture  completion  time  in  the  case  of  the  hole  models  is 
dependent  only  on  the  fracture  length,  and  that  the  individual  var¬ 
iations  in  the  hole  diameter  and  the  hole  spacing  parameters  will  have 
no  Influence  on  the  fracture  time  provided  the  fracture  length  is 
kept  constant.  This  is  obvious  from  the  approximately  parallel 
curves  in  Figure  53  which  show  the  fracture  time  versus  the  hole 
spacing  plots  for  the  hole  models  with  different  diameters. 

As  seen  from  the  temperature  plots  for  the  holemodels  shown  In 
Figures  47-49  ,  a  major  portion  of  the  fractured  voliane  does  not 
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Figure  51.  Fracture  Zones,  Model  22B 
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Figure  52.  Fracture  Zones,  Model  22C 
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Figure  S3.  Fracture  Time-Spacing  Curves  for  Hole  Model 
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•)(p«r1«nc«  iny  twnperiturt  ch«ngt,  ami  hfnca»  tho  In  this 

rtglon  art  out  anti  rely  to  the  load  vector  resulting  from  the  thermally 
Induced  constraint.  Thus,  C*  explanation  for  the  sharp  Increase 
In  fracture  completion  time  with  Increased  hole  spacing  In  the  cate 
of  the  slot  models  also  applies  to  the  hole  models. 

Comparison  of  Figures  40  and  S3  shows  a  remarkable  similarity 
between  the  parallel  and  the  perpendicular  fracture  characteristics » 
although  for  a  given  problem  geometry^  the  perpendicular  fracture 
time  using  the  hole  model  approximation  Is  much  higher  than  the  parallel 
fracture  time.  However*  as  mentioned  earlier,  the  hole  model  analysis 
Involves  small  hole  depths  without  any  convection  depths  and  thus, 
assumes  a  plane  stress  condition.  In  practice,  however,  the  holes  arc 
rather  deep  and  the  problem  appraoches  a  plane  strain  condition. 

Thus,  the  actual  perpendicular  fracture  times  as  compared  to  the 
actual  parallel  fracture  times  will  not  be  as  high  as  theoretically 
predicted.  It  should  also  be  realized  that  while  the  convection  depth 
will  greatly  reduce  the  perpendicular  fracture  time  as  1n  the  case 
of  the  parallel  fracture,  any  Increase  In  the  convection  depth  will 
mean  the  heat  source  will  have  a  greater  burden  against  which  to  open 
cracks  between  the  holes. 


Chapter  IX 

HEAT  TRANSFER  IN  ROCK  MATERIALS  OF  SIMPLE  GEOMETRIC  SHAPE 


The  objective  of  this  part  of  the  research  was  to  determine  the 
temperature  distribution  within  a  rock  |iater1al  resulting  from  an  in¬ 
ternal  cylindrical  heater.  The  Information  obtained  should  provide 
a  means  of  evaluating  heater  performance  as  well  as  Insight  Into  the 
appropriate  boundary  conditions  at  the  Internal  rock  surface  due  to 
primarily  radiant  heat  sources.  The  heaters  were  constructed  from 
Kanthal  wire  colls  wound  on  alumina  cylinders. 

Discussion  of  Analytical  Problem 

For  thermal  rock  fracturing  processes  there  Is  a  need  for  a 
method  for  solving  transient  heat  transfer  problems  with  variable 
thermal  properties  and/or  radiation  boundary  conditions.  The  develop 
ment  and  operation  of  a  three-dimensional  code  would  be  expensive  for 
complete  nonlinear  conditions.  The  equations  of  heat  transfer  for 
this  type  of  problem  are: 


v(k7T)  “  pcf^ 

(9.1) 

"'T-U  '  K  11 

(9.2) 

e„(r-T,')  =  k|L 

(9.3) 

The  partial  differential  equation  Is  nonlinear  due  to  the  temperature- 
dependence  of  the  material  properties.  The  radiation  boundary  condi¬ 
tion  Is  also  nonlinear. 

There  are  two-dimensionless  codes  available  which  will  solve  the 
linear  transient  heat  transfer  problem.  A  nonlinear,  two-dimensional 
steady  state  code  is  also  available.  The  available  two-dimensional, 


transient  code  is  being  modified  to  solve  problems  with  temperature- 
dependent  material  properties. 

Heat  Transfer  Study 

The  method  of  heat  transfer  from  the  Kanthal  heaters  to  the  rock 
surfaces  Is  a  complicated  combination  of  convection,  radiation  and  con¬ 
duction.  An  exact  analysis  or  even  an  experimental  measurement  of  the 
temperature  distribution  Inside  a  large  rock  mass  would  be  difficult. 
However,  the  most  Important  function  which  a  boundary  condition  performs 
Is  to  determine  how  much  heat  actually  enters  the  rock  body.  The  mode 
of  heat  transfer  which  is  employed  at  the  boundary  Is  of  secondary  im¬ 
portance. 

Discussion  of  Experimental  Problem 

An  Ideal  experiment  gives  direct  readings  of  physical  events  without 
analytical  or  experimental  manipulation.  Since  it  is  not  possible  to 
measure  the  temperature  directly  inside  a  rock  specimen,  the  experiment 
should  be  designed  as  close  to  the  ideal  experiment  as  possible.  To 
minimize  the  analytical  manipulation,  a  one-dimensional  model  was  chosen. 

A  basic  heat  transfer  study  was  devised  to  determine  the  heat  flux 
entering  the  rock  as  a  function  of  time.  This  study  consisted  of  de¬ 
signing  an  experimental  test  which  would  have  essentially  one-dimensional 
heat  flow  and  then  measuring  the  temperature  as  a  function  of  time. 

The  temperature  history  will  be  used  in  future  work  to  estimate  the 
surface  temperature  as  a  function  of  time.  Thus,  a  time-dependent  tem¬ 
perature  boundary  condition  will  be  available  for  further  analytical 
studies. 

The  shape  of  the  Kanthal  heaters  is  cylindrical  and  the  natural 
shape  for  the  rock  specimen  would  also  be  cylindrical.  Figure  54 
shows  the  heater,  rock  specimen,  four  thermocouples  and  two  end  plugs. 

The  end  plugs  are  used  to  prevent  excessive  convective  currents. 


KANTHAL  WIRE 
WOUND  HEATER 


END  PLATE 


BASALT  CYLINDER 


4  ALUMEL-CHROMEL 
THERMOCOUPLES 


END  PLATE 


Figure  54.  Experimental  Configuration  for  Cylinder  Tests 
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Since  the  basic  shape  was  determined,  the  actual  dimensions  had 
to  be  determined.  A  special  heater  was  built  whose  outside  diameter 
was  somewhat  less  than  one  inch.  Also,  the  largest  core  drill  presently 
available  at  the  RMERC  has  a  four-inch  inside  diameter.  The  outside 
diameter  of  the  rock  specimen  then  was,  from  logistical  considerations, 
four  inches  and  the  inside  hole  diameter  was  one  inch.  The  length  of 
the  rock  specimen  had  to  be  determined  so  that  nearly  one-dimensional 
heat  flow  is  experienced  at  the  mid-plane  of  the  cylinder.  This  length 
was  to  be  determined  by  considering  two  important  factors: 

a)  A  steady  state  analysis  will  show  the  two-dimensional  effects 
more  strongly  than  a  transient  analysis,  and 

b)  Through  numerous  tests  with  Kanthal  wire  heaters,  the  rock 
surface  has  never  been  melted. 

Thus,  if  a  two-dimensional  model  was  designed  based  on  a  steady  state 
analysis  and  an  absolute  maximum  temperature  condition,  (i.e.,  the  sur¬ 
face  at  the  material  melt  temperature),  then  the  design  should  be  valid 
for  transient  conditions.  The  method  of  study  was  to  select  grids 
representing  one-eighth  (from  symmetry  and  thermocouple  placement)  of 
the  rock  cylinder  with  boundary  conditions  of  a  melting  inner  surface 
and  convection  on  the  outer  surface  and  find  the  steady  state  tempera¬ 
ture  distribution.  A  typical  grid  shown  in  Figure  55  was  drawn  for  five 
length-to-depth  ratios  (L/D«2-10)  and  the  temperatures  were  calculated 
with  constant  conductivity  and  then  temperature-dependent  conductivity. 
Dresser  basalt  was  chosen  as  the  rock  material  for  two  reasons: 

a)  Material  properties  are  g^ven  in  Reference  64. 

b)  Basalt  is  easier  to  machine  than  most  other  hard  rocks. 

The  value  of  h  was  given  in  Reference  68  as  h  »  1.57  Btu/hr-ft^-^F. 

For  the  constant  conductivity  case  the  conductivity  corresponding  to 
zero  degrees  was  used  as  a  starting  point  of  iteration. 
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Figure  55.  Typical  Finite  Element  Grid  for  Cylinder  Model  (L/l>10) 
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The  isother  produced  by  the  zero  degree  value  of  conductivity 
and  the  final  value  of  conductivity  are  typically  shown  in  Figures  56 
and  57.  The  radius  of  the  center  hole  is  1/2  in.  and  the  data  indicates 
that  within  4  mm  of  the  surface,  the  heat  flow  is  essentially  one¬ 
dimensional  for  all  of  the  ratios  of  L/D. 

A  bound  on  the  error  due  to  placing  a  thermocouple  in  the  rock 
cylinder  was  needed.  A  new  grid  (Fig.  58)  was  made  which  would  repre¬ 
sent  1/8  of  the  circular  section  of  the  cylinder.  Then  the  material 
was  designated  as  solid  basalt  and  as  basaU  with  a  thermocouple,  in¬ 
sulator  and  a  small  air  space.  The  tip  of  the  thermocouple  was  allowed 
to  be  at  1  mm,  2  mm,  3  mm,  and  finally  at  4  mm.  Boundary  conditions 
were  Newton's  law  of  cooling  on  the  outer  surface  and  the  basalt  melt¬ 
ing  temperature  (1250"C.)  on  the  inner  surface.  The  values  of  conduc¬ 
tivity  for  air,  chromel-alumel  and  omegatite  were  assumed  constant 
and  are  given  in  Table  13. 


TABLE  13 


M=»terial 

Air 

Chromel-Alumel 

Omegatite 


Conductivity 
Btu/hr  ft  “F 

.0303 

24.08 

14.5 


Figures  59-  62  show  typical  isotherms  for  the  constant  conductivity  case 
and  for  the  temperature-dependent  case. 

Table  14  gives  values  of  the  error  produced  by  the  thermocouple 
in  the  rock  specimen  by  depth.  The  maximum  value  of  6.12  percent  error 
is  an  acceptable  amount  considering  that  a  steady  state  analysis  was 
used  and  as  an  assumption  in  solving  the  grid  of  Figure  58,  all 


Figure  58.  Typical  Finite  Element  Grid  for  Cylinder  Model  with  Thermocouple 


Teeperature  Distribution  for  Yariable  Themal  Conductivity  vithcvt  Tbc! 


•lemtnti  wtr#  t«ktn  i«  hftvInQ  unit  th1ckn»i«,  Thli  assumption  of  unit 
thlckntii  would  mionlfy  tho  ictunl  oi^ror  productd  by  t.h»  thtmocouple 
thlckntts. 

TAfilE  14 


Depth 

Tempereture 
Solid  Wedge 

Temperature 
Wedge  with 
Thennocouple 

Temperature 

Difference 

%  Error 

1  mm 

?lb2.783 

2048.450 

114.333 

5.27 

?  mm 

2118.149 

2020.116 

98.033 

4.67 

3  mm 

2077.643 

1950.51? 

127.131 

6.12 

4  mm 

2040.612 

1950.820 

09.792 

4.31 

The  next  step  wes  tn  construct  the  rock  specimens  end  measure  the  tem¬ 
peratures  produced  by  the  Kanthal  heaters. 

From  the  temperature  nteasureinents  a  ono>diinens1ona1  analysis  is 
to  be  made  to  determine  the  rock  surface  temperature  as  a  function  of 
time.  This  tentperature  function  will  then  be  used  in  the  two-dimensional 
programs  to  calculate  the  temperature  distribution  within  rock  materials. 
Experimental  Results 

Figure  63  shows  the  temperatures  of  the  rock  at  1  iiin»  depth  for  three 
selected  power  settings.  Figure64  shows  the  outside  surface  tempera¬ 
tures  for  the  same  power  settings  of  Figure  63.  Figure 65  shows  the 
four  thermocouple  readings  for  a  maximum  power  setting. 


Fig.  64.  Outside  Surface  Temperature  Variation  with  Time  for 
Dresser  Basalt  Cylinder. 
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Chapter  X 

LABORATORY  AND  FIELD  FRACTURE  TESTS 

Introduction 

Research  on  thermal  fragmentation  of  hard  rock  was  begun  at  UMR 
in  1969  and  has  been  carried  out  in  three  phases: 

1.  Surface  heat  weakening  and  chipping. 

2.  Mechanical  slot  cutting  and  thermal  ker*  removal. 

3.  Investigation  of  a  thermal  round  sim!'  .'.  in  design  to  an 
explosive  round. 

In  phase  one  the  surface  of  blocks  of  select^'>d  igneous  and  meta- 
morphic  rock  was  subjected  to  heat  from  a  small  flame  jet  torch  caus¬ 
ing  shallow  fractures.  The  rock  was  then  chipped  off  with  a  traversing 
jackhammer.  It  was  found  that  while  the  surface  of  most  of  the  rocks 
tested  was  weakened  by  heat,  rock  removal  by  impact  was  not  enhanced 
by  the  heating  process. 

In  phase  two  slots  1-3/4  inches  wide  by  six  inches  deep  were  cut 
in  the  surface  of  three-foot  cubes  of  granite,  quartzite  and  basalt, 
with  six-inch  ridges  in  between  the  slots.  When  heat  was  applied  to 
the  bottom  of  the  slots  the  ridges  between  were  readily  fractured  off 
from  the  rock  face. 

Phase  three  is  being  devoted  to  the  investigation  of  the  feasibility 
of  developing  a  thermal  round  which  w’!!  function  in  principle  in  a 
manner  similar  to  an  explosive  round.  The  application  of  heat  at  the 
bottom  of  deep  slots  within  the  rock  in  phase  two  demonstrated  clearly 
that  the  rock  is  much  more  effectively  fractured  if  the  heat  energy 
is  deposited  within  the  rock. 
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Thermal  Stress  Fracture  of  Rock 

It  Is  a  well  established  principle  of  bias  mg  that  a  cylindrical 
column  of  explosive  must  have  one  or  more  free  faces  to  which  to  break. 
With  stress  waves  generated  by  explosives  the  free  face  provides  a 
meens  for  spallation  to  occur  and  subsequent  fracture  to  the  explosive 
will  take  place  in  a  properly  designed  round.  Each  hole  is  placed  so 
that  it  will  provide  the  necessary  free  faces  to  which  other  holes  may 
break. 

The  process  of  fracture  from  thermal  inclusions  is  different  in 
several  ways  from  fracturing  by  explosives.  When  a  heater  is  placed 
in  a  hole  in  the  rock  a  thermal  stress  is  generated  around  the  hole. 

The  immediate  periphery  of  the  hole  is  in  biaxial  compression  with  a 
tensile  stress  beyond.  The  magnitude  of  the  tensile  stress  depends 
upon  the  properties  of  the  rock,  the  temperature  and  the  radius  of 
the  hole. 

The  nonstationary  state  of  stress  due  to  a  constant  temperature 
suddenly  applied  to  a  spherical  cavity  in  an  infinite  medium  was 
analyzed  theoretically  by  Sternberg  (83).  The  radial  displacement 
of  the  cavity  surface  is  zero  for  all  time.  The  temperature  drops 
rapidly  with  distance  away  from  the  cavity  surface  (Fig. 66)  and  does 
not  penetrate  beyond  3  cavity  radii  for  t  =  3  (T  =  -v.t  =  time, 

K  =  diffusivity  and  a  =  cavity  radius).  Radial  displacement  increases 
with  time  as  does  the  radial  stress,  which  is  always  compressive. 

The  tangential  stress  at  a  given  r/a  greater  than  one  is  compressive, 
becomes  tensile  and  then  compressive  again  with  increase  in  time. 

The  magnitude  of  the  tensile  stress  decreases  with  r/a. 
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For  constant  diffusivity,  the  magnitude  of  the  tangential  tensile 
stresses  is  largest  for  small  values  of  r/a,  or  for  large  values  of  the 
cavity  radius,  and  increases  linearly  with  the  cavity  temperature. 

The  curves  in  Figure  66 may  serve  as  approximate  guides  to  the  ther¬ 
mal  stress  history  around  the  thermal  inclusion  created  by  carbon  elec¬ 
trodes  in  a  hole  in  a  solid  granite  face.  However,  the  cavity  in  granite 
increases  in  size  with  time  and  some  of  the  rock  properties  are  tempera¬ 
ture  dependent.  Also,  the  stress  pattern  is  altered  by  the  presence 
of  the  drill  hole. 

The  analysis  by  Sternberg  (83)  agrees  qualitatively  with  that  by 
Nixon  (4),  the  latter  being  the  (cylindrical)  finite  element  thermal 
stress  analysis  for  electron  beam  penetration.  That  is,  the  magnitudes 
of  compressive  stresses  are  ten  times  those  of  the  tensile  stresses,  and 
for  longer  heatin^  times  the  zone  of  tension  extends  considerably  beyond 
the  heater  zone. 

The  initiation  and  propagation  of  fractures  in  rock  in  the  configura¬ 
tions  used  herein  due  to  internal  thermal  stresses  are  controlled  by  two 
important  factors,  (1)  the  magnitude  and  direction  of  the  tensile  stresses, 
and  (2)  strain  relief  due  to  compressibility  of  the  rock  or  nearby  trac- 
tiorifree  surfaces. 

Displacement  Relief  and  Fracture 

When  a  high  temperature  heat  source  is  placed  in  a  cylindrical  or 
spherical  cavity  in  an  infinite  (rock)  medium  thermal  stresses  are  in¬ 
duced  in  patterns  as  indicated  above  (Fig. 66).  This  may  be  shown  graphi¬ 
cally  for  a  given  time  after  a  temperature  T^  has  been  applied  to  the 
surface  of  the  cavity.  When  a„  in  the  tensile  zone  exceeds  the  tensile 
strength  of  the  material  a  crack  normal  to  the  stress  direction  will  form. 
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However,  the  distance  it  will  propagate  and  the  amount  it  will  widen 
depends  upon  strain  displacement  relief  in  body  which  may  permit  dis¬ 
placements.  In  nonporous  hard  granite  it  has  been  found  that  the  rock 
cannot  readily  move  normal  to  the  crack  unless  there  are  three  planes 
or  faces  of  strain  relief  within  critical  distances.  (In  this  respect 
thermal  stresses  function  differently  than  explosively  induced  stresses 
which  will  create  craters  with  only  one  free  face).  With  less  than  mini¬ 
mum  displacement  i^elief  only  hairline  fractures  are  formed,  extending 
into  the  cavity  and  a  short  distance  into  the  rock.  However,  where 
free  faces  are  present  and  multiple  heater  holes  furnish  stress  fields 
which  reinforce  each  other,  cracks  propagate  for  longer  distances  and 
widen  to  1/4  inch  or  more  due  to  the  expansion  of  the  rock  around  the 
cavities.  The  blocks  of  rock  formed  by  the  fracture  are  easily  removed. 
That  is,  thermal  forces  have  demonstrated  a  significant  amount  of  follow- 
through  in  dislodging  the  fractured  blocks  in  experimental  tests. 
Experimental  Results  -  Laboratory 

Thirumalai  (  8)  had  shown  earlier  that  thermal  inclusions  generated 
by  dielectric  current  would  readily  fracture  small  blocks  of  basalt  but 
not  quartzite.  Thus,  a  positive  method  of  heating  was  desired  whi-cl'  is 
independent  of  the  electrical  properties  of  the  rock.  In  the  first  ex¬ 
periments  coiled  nichrome  wire  was  wrapped  in  a  spiral  on  alundum  cores 
and  placed  in  pneumatically  drilled  holes  in  three-foot  cubes  of  quart¬ 
zite  and  granite.  The  heaters  successfully  fractured  the  rock  parallel 
and  perpendicular  to  the  heater  holes,  but  the  nichrome  wire  burned  out 
at  Sustained  temperatures  of  600®  to  700®C.  Kanthal  wire  resistance 
heaters  were  employed  at  temperatures  of  900®  to  1200°C.  These  fractured 
the  rock  more  readily,  but  also  burned  out  in  heater  holes  in  granite 
at  sustained  temperatures  of  1000®C.  They  were,  however,  employed 
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successfully  to  fracture  granite  in  place  with  ten-inch  burden  on 
the  holes  and  with  three  directions  of  strain  relief. 

Coiled  Kanthal  wire  electric  resistance  heaters  one  foot  long 
and  1-3/4  inches  in  diameter  will  generate  from  3.5  to  4.5  kw  at  tem¬ 
peratures  from  gOO^C  to  1200°C.  They  will  successfully  fracture  a 
strong  granite  in  place  with  10-inch  burden  and  spacing  to  a  depth  of 
15  inches  where  three  free  faces  are  present.  However,  there  is  a 
little  follow-through  by  radial  compression  and  fractured  blocks  are 
sometimes  difficult  to  remove.  Also,  in  spall  able  rocks  such  as  granite, 
flakes  of  minerals  spall  off  from  the  cavity  walls  and  clog  up  the  heater 
coils,  causing  them  to  burn  out  after  15  to  20  minutes  of  use. 

The  results  obtained  with  coiled  wire  heaters  indicated  that  higher 
temperatures  might  be  desirable,  and  that  the  heater  should  be  of  an 
economical  type  that  could  furnish  high  temperatures  in  a  hole  in  rock 
without  breaking  down.  Refractory  metal  heaters  of  a  type  similar  to 
those  used  for  the  LASL  melting  drill  (11)  were  considered,  but  the  cost 
was  too  high. 

The  initial  internal  heater  tests  were  made  with  30-inch  cubes  of 
quartzite  and  granite.  A  slot  3/4  inches  wide,  12  inches  long  and  8 
inches  deep  was  cut  in  one  face  of  a  cube  of  quartzite  with  six  1-1/4 
inch  diameter  heater  holes  drilled  on  either  side  of  the  slot  (Fig.  67). 
The  first  tests  were  made  with  nichrome  coiled  wire  heaters  at  a  tempera¬ 
ture  of  650®C  and  3  kw  power  level  in  one  set  of  holes.  The  heat  flux 
was  about  350  watts  per  inch  of  hole.  Cracks  connecting  the  holes  ap¬ 
peared  in  about  1-1/2  minutes.  Other  cracks  appeared  and  grew  until 
one  of  the  heaters  burned  out  in  2-1/2  minutes.  Smaller  heaters 
(150  watts/inch)  were  placed  in  the  second  row  of  holes  (Fig.  67). 
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Figure  67,  Slot  and  Hole  Pattern 


Figure  68.  Vertical  Fracture  Pattern  Parallel  to  Face 
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Initial  fractures  appr'^ed  In  3  mlnuUs  and  they  continued  to  widen 
to  15  minutes.  This  application  of  h»at  resulted  In  the  block  being 
fractured  Into  several  large  pieces.  In  the  slot  area  fractures  paral¬ 
lel  to  the  face  occurred  about  4  Inches  from  the  bottoms  of  the  drill 
holes  (fig.  68, 69). This  type  of  fracture  occurred  consistently  In  eight 
tests  of  all  of  three  kinds  of  rock,  i.e.,  quartilte.  granite,  and 
bisalt. 

Electric  arcs  with  carbon  electrodes  were  suggested  as  a  source  of 
heat  and  these  were  tested.  Depending  upon  the  slie  of  electrodes  and 
other  factors,  from  10  kw  to  20  kw  of  power  are  easily  obtainable  In  a 
2-1nch  drill  hole.  Temperature  of  the  arc  Is  in  the  neighborhood  of 
7,000®C.  While  carbon  electrodes  are  brittle,  with  reasonable  care  they 
can  be  used  under  dlfficul  vorKing  conditions. 

In  an  Initial  test  two  l/2-1nch  electrodes  were  placed  In  a  2-1/4  Inch 
diameter  hole  In  a  granite  block  with  about  5  to  6  kw  power  from  a  250 
amp  arc  welder.  The  block  broke  with  a  loud  noise  1n  6  minutes,  using 
about  0.6  kwhr,  and  producing  a  fracture  completely  through  the  block 
parallel  to  the  heater  hole,  as  well  as  breaking  the  top  half  of  the 
block  into  three  large  segments  (Fig.  70).  The  heat  developed  a  hollow, 
bulb-shaped  thermal  inclusion  about  eight  Inches  long  and  a  four  inch 
maximum  diameter.  The  mass  of  the  inclusion  was  made  up  of  partially 
fused  granite  and  chemically  altered  mineral  flakes,  of  380  cc  volume, 
representing  0.2  percent  of  the  volume  of  the  rock.  The  time  and  energy 
requirements  for  tests  at  UMR  and  those  of  Thlrumalai  (  8)  Indicate  that 
the  cube  root  law  of  energy  scaling  may  apply.  Similar  tests  were  con¬ 
ducted  on  30-inch  blocks  of  quartzite  and  basalt  with  the  following 
results: 
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The  Mlsioui'1  red  granite  la  U&a  brittle  and  stronger  than  the 
quartilte,  which  broke  easily,  probably  due  to  the  high  coefficient 
of  thermal  expansion  of  quartz  and  Its  change  from  u  to  h  quartz  as 
Its  temperature  Is  Increased.  Basalt  fractured  readily  even  though 
It  contains  little  or  no  free  quartz.  Its  lava  being  very  liquid. 

Field  Tests 

All  equipment  used  In  the  field  was  made  portable  and  included  as 
of  the  final  tests  reported  below,  a  300  cfm  cotnpressor,  a  100  kw  diesel- 
driven  generator,  four  600  ampere  dc  welders,  pneumatic  drills,  three 
trucks,  a  drill  post  mounted  on  a  truck  bumper,  high  frequency  units, 
carbon  arc  heaters,  and  accessory  equipment.  Because  of  the  Improvised 
nature  of  the  drilling  equipment,  progress  has  been  slower  than  that 
which  could  he  achieved  with  an  autanated  jumbo,  a  new  compressor,  etc. 

The  results  of  the  Internal  heating  tests  on  3-ft  cubes  of  quart¬ 
zite  and  granite  showed  conclusively  that  In  situ  conditions  of  confine¬ 
ment  and  other  mass  effects  could  not  be  easily  simulated  In  the  labora¬ 
tory.  It  was  therefore  planned  to  conduct  full-scale  field  tests  in 
an  abandoned  quarry  In  Missouri  red  granite  located  near  Granltevllle, 
Missouri.  The  granite  In  the  old  working  faces  of  the  quarry  has  a 
wide  variety  of  structural  features  suitable  for  In  place  thermal  frag¬ 
mentation  tests,  varying  from  solid  masses  to  rock  which  are  heavily 


16a 

JolhteU  «ml  »««st  A|>p»rtr1ng  f«c«»  worw  choi^n 

for  InUUl  tosts.  Theso  i-on^Ut^d  of  the  followlmj  configuration!: 

Type  1.  A  tingle  tint  l^VA  x  16  In.  with  three  heater  holei 
on  either  tide  of  the  tlot  (Fig, 71),  with  cylindrlcel 
resUtance  wire  heeters. 

Type  i'.  Two  mutually  perpendicular  slot!,  2  In,  x  6  ft  In  the 
rock  face  with  heater  holes  parallel  to  the  slots 
(Fig.  72),  with  cylindrical  resistance  wire  heaters, 

Type  3.  Three  heater  holes  in  line  with  three  free  faces 
(Fig.  73),  with  carbon  arc  heaters. 

Type  4.  Four  heater  holes  In  a  square  pattern  (Fig.  74),  with 
carbon  arc  beaters. 

Type  1.  Resistance  wire  heaters  at  900  to  1000“C  and  4.5  kw  energy 
level  were  placed  In  the  holes  on  each  side  of  the  slot  (Fig. 71)  for 
a  period  of  12  to  18  minutes.  Several  fractures  formed  between  the 
heater  holes  and  radially  from  some  of  them,  but  only  a  few  surface 
fragments  could  be  removed  with  a  pry  bar.  The  rock  between  the  holes 
and  the  slot  were  broken  out  with  sotne  difficulty  with  a  Darda  rock 
sputter.  Indicating  that  more  displacement  relief  may  be  required  for 
easy  removal.  Also  some  of  the  rock  partings  between  the  holes  In  the 
slot  had  not  been  completely  removed  and  these  offered  some  resistance 
to  fracturing  and  rock  removal. 

(It  was  found  later  In  Type  4  tests  that  the  fracturing  and  rock 
removal  was  markedly  enhanced  when  internal  pressure  was  applied  with 
a  rock  splitter  during  the  heating  of  the  rock.  This  type  of  experi¬ 
ment  is  planned  in  future  work  for  the  cut  portion  of  a  thermal  round 
to  create  additional  free  faces  for  displacement  relief.) 
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Figure  71.  Incanplete  Slot  with  Htdter  Holes  Used  in  First  Field  Test 


Figure  72.  Two  Slot  Cut  with  Example  Heater  Holes  for  Tunnel  Round 


171 


Type  2.  Two  mutually  perpendicular  slots  (Fig. 72)  1-3/4  In.  wide 
by  6  ft  long  by  30  In.  deep  were  cut  In  a  granite  face  and  heater 
holes  2  In.  diameter  were  drilled  parallel  to  the  slots  with  a  bur¬ 
den  and  spacing  of  8  In.  to  10  In.  With  small  burden  and  spacing 
much  of  the  rock  fractured  was  easily  removed  with  a  pry  bar,  In 
some  cases  the  fracture  process  being  aided  by  joints  parallel  to 
the  rock  face.  However,  the  cylindrical  resistance  heaters  failed 
after  about  12  minutes  at  a  load  at  4.5  kw. 

With  the  drilling  equipment  on  hand  the  slots  were  time  consum¬ 
ing  to  cut  and  a  natural  approximately  cubic  corner  of  a  rock  section 
In  the  quarry  was  selected  for  subsequent  tests  utilizing  carbon  arc 
heaters . 

Type  3.  In  laboratory  tests  with  arc  heaters  a  250  ampere  (7.5  kw) 
power  source  (arc  welder)  was  employed.  For  field  work  a  higher  power 
level  was  desired,  and  an  ac,  three  phase,  100  kw  generator  could  be 
used  for  three  12  kw  heaters.  However,  difficulties  were  encountered 
in  maintaining  arc  stability.  Two  tests  were  successfully  <'onducted 
with  three  holes  in  line  using  ac  current,  with  12  in.  and  14  in.  bur¬ 
dens  (Fig.  75).  Visible  fractures  appeared  in  about  4  to  6  minutes, 
and  were  widened  to  about  1/8  in.  after  the  heaters  were  left  on  for 
about  20  minutes.  The  blocks  were  easily  removed  with  a  pry  bar. 

The  thermal  stress  fields  created  around  the  heated  cavities 
mutually  augmented  the  affects  of  each  other,  creating  fractures  be¬ 
tween  the  holes  and  a  fracture  perpendicular  to  the  holes,  extending 
to  the  second  and  third  faces  and  into  the  rock  normal  to  the  axes  ^^f 
the  holes,  The  latter  fracture  aids  the  formation  of  new  fractures 
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behind  additional  blocks  of  rock  which  are  broken  off.  This  fracture 
does  not  form  normal  to  the  second  free  face  but  curves  away  from  the 
first  face.  Also,  a  set  of  fractures  appear  to  be  formed  tangent  to 
the  thermal  inclusion,  furnishing  the  primary  separation  surface  near 
the  heated  cavity. 

One  test  with  two  heater  holes  was  conducted  with  20  kw  per  hole. 

No  improvement  in  fracture  time  was  observed  although  two  holes  may  be 
less  than  minimal  for  effective  fracture. 

Type  4.  Theoretical  analyses  of  the  thermal  stress  distribution  caused 
by  multiple  cylindrical  and  spherical  thermal  Inclusions  shows  that  if 
stress  concentrations  can  be  induced  in  the  proper  orientation  and  loca¬ 
tion,  direction  and  location  of  the  fracture  surface  created  can  be  con¬ 
trolled.  Thus,  the  stress  concentrations  around  the  closed  ends  of  cy¬ 
lindrical  holes  will  result  in  a  plane  of  fracture  if  the  holes  all 
bottom  in  the  same  plane.  Likewise,  if  spherical  inclusions  lie  in  the 
same  plane,  the  stress  concentrations  will  mutually  act  to  create  a 
plane  surface  of  fracture  extending  between  and  beyond  the  heated  cavi¬ 
ties.  The  cylindrical  holes  also  serve  to  concentrate  stresses  in  such 
a  manner  as  to  cause  fractures  to  form  parallel  to  the  axes  of  the  holes. 

While  high  temperature  gradients  in  the  rock  will  cause  correspond¬ 
ingly  large  stress  gradients,  the  strain  and  the  displacement  induced 
by  thermal  expansion  accompanying  the  critical  fracture  stress  levels 
is  relatively  small,  and  fractures  may  be  formed  in  the  rock  before 
they  are  visible  to  the  naked  eye.  Hence,  in  the  later  tests  one  or 
two  Darda  rock  splitters  were  placed  in  holes  in  the  rock  between  the 
heater  holes  (1)  to  aid  in  fracture  initiation,  (2)  to  exploit  the  frac¬ 
tures  fonried  by  the  thermal  stresses,  and  (3)  to  move  fr-!ctured  blocks 
of  rock  sufficiently  so  that  they  could  be  removed. 
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Four  500  amp  (continuous  duty)  welders  were  obtained  as  power 
sources,  and  holes  drilled  on  a  rectangular  pattern  (Fig.  74).  The 
variables  investigated  were  the  burden,  hole  spacing,  depth  of  hole, 
and  time  to  fracture.  Each  heater  was  operated  at  an  average  of  12  kw 
power  level,  or  the  maximum  continuous  level  accommodated  by  single 
welders.  One  experiment  was  conducted  with  two  heaters  in  two  holes 
at  800  amps  each  to  ascertain  some  of  the  effects  of  higher  power  level 
as  noted  above. 

Test  No.  5  (Table  15)  was  conducted  with  heaters  only.  The  initial 
observable  fractures  appeared  in  about  4  minutes,  and  they  were  widened 
considerably  in  10  minutes  total  heating  time.  The  block  was  later 
removed  with  a  splitter.  The  next  two  tests  with  18-inch  spacing  were 
made  with  a  splitter  in  one  hole  only  (one  splitter  was  available  at  the 
time).  Initial  fractures  appeared  in  about  3  minutes  and  subsequent 
fracturing  was  materially  aided  by  the  single  splitter,  although  in 
one  test  the  horizontal  fracture  occurred  at  the  level  of  the  splitter. 
Vertical  fractures  occurred  parallel  to  the  working  face  and  through 
the  heater  holes. 

Two  subsequent  tests  with  24-inch  spacing  were  made  with  two 
rock  splitters.  In  the  first  test  the  arc  in  No.  4  hole  operated  only 
a  short  time  and  the  block  fractured  diagonally.  The  second  test  in 
this  series  functioned  in  a  normal  manner,  the  fracture  times  being 
decreased  considerably  by  the  mechanical  stress  induced  by  the  splitters. 
The  results  of  the  test  indicate  that  24-inch  spacing  may  be  approach¬ 
ing  the  practical  limit  for  a  12  kw  power  level  in  each  heater.  A 
total  of  seven  experimefits  was  conducted  at  the  second  site  (Figs.  76  &  77). 


Figure  76.  Comnosite  of  Experiments  Four  Through  Nine  Showing 
Successful  Controlled  Block  Fracture 
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Fig.  77.  Four  Ciectric 
arc  heaters  in  place  with 
one  Darda  rock  splitter 
(wrapped  in  asbestos  cloth). 


Fig.  78.  Drill  holes  and  frac¬ 
ture  patterns  for  heaters  and 
splitters.  A''l  blocks  were 
easily  fractured  and  removed. 
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TABLE  15 

Summary  of  Field  Tests 

Hole 


Test  No. 

Geometry 

Depth 

spacing 

Burden 

Power/Hole 

Time 

kw 

1 

Line 

3  holes 

16 

12 

12 

8-10 

15 

2 

Line 

3  holes 

22 

4 

4 

8-10 

15 

3 

Line 

3  holes 

30 

24 

24-30 

8-10 

25 

4 

Line 

(2  holes) 

30 

18 

16 

20-25 

25 

5 

Square 

4  holes 

28 

18 

10 

12 

10 

3-f 

6 

Square 

4  holes* 

28 

18 

10 

12 

10+ 

4-f 

7 

Square 

4  holes* 

28 

18 

10 

12 

11  + 

8 

Square 

4  holes** 

32 

24 

12 

12 

15++ 

9 

Square 

4  holes** 

32 

24 

12 

12 

ITH- 

★ 

1  hole 

5  min. 

only 

ifr* 

1  hole 

8  min. 

only 

+  with  one  splitter 
++  with  two  splitters 
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Round  Design  and  Displacement  Relief 

If  one  considers  a  section  of  a  heat''d  cylindrical  hole  near  the 
corner  of  a  quadrant  of  a  large  elastic  body,  after  the  heat  has  been 
applied  for  a  certain  time  there  is  a  biaxial  compressive  stress  im¬ 
mediately  around  the  hole  which  becomes  tensile  (tangential)  with  in¬ 
crease  in  distance  from  the  hole.  With  increase  in  time  the  width  of 
the  compressive  zone  increases  as  does  that  of  the  tensile  zone,  the 
magnitude  of  the  latter  decreasing  with  distance  of  constant  hole  diameter. 
If  the  medium  were  infinite,  radial  cracks  could  form,  their  orientation 
depending  upon  properties  and  imperfections  of  the  medium.  However,  if 
there  is  one  free  face  within  a  distance  such  that  the  magnitude  of  the 
tensile  stress  is  affected  by  the  free  surface,  a  crack  will  propagate 
from  the  hole  normal  to  the  free  surface,  but  it  will  not  widen  because 
of  the  constraining  effect  of  the  mass.  If  two  surfaces  exist  within  a 
critical  distance,  then  two  cracks  will  propagate  from  the  hole  normal 
to  the  two  free  surfaces  and  will  continue  to  open  up  if  there  is  suffi¬ 
cient  displacement  relief.  One  of  the  objectives  of  fur Iner  experimenta¬ 
tion  is  to  determine  more  exactly  the  value  of  optimum  ourden  for  given 
rock  properties,  hole  size  and  energy  level. 

For  three  dimensions,  an  octant  of  an  infinite  mass  s  considered 
with  multiple  holes  drilled  parallel  to  one  face.  The  problem  now  be¬ 
comes  one  of  determining  the  optimum  distance  from  a  free  surface 
(burden),  the  spacing,  and  the  depth  of  the  holes.  The  holes  act  as 
loci  for  stress  concentration  and  fracture  guides.  However,  for  this 
mode  three  planes  within  critical  distanc's  of  the  thermal  inclusion 
are  needed  for  controlled  fr'?rture  and  most  effi.,ient  breakage.  For 
breakage  of  the  cut  portion  of  a  round  in  a  normal  tunnel  face  where 
no  slots  are  cut,  the  same  principle  applies  but  only  one  free  face 
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Is  presont,  and  cut  d»s1()ns  to  those  usttd  In  bUsting  Are  r«< 

quirod.  Such  cut  design  Is  to  be  investigated  in  future  research. 
Thermal  Round  Operational  Concents 

There  appear  to  be  several  concepts  of  thermal  rounds  which  will 
function  to  break  the  rock,  The  primary  problem  is  providing  initial 
displacedient  relief.  The  possible  "cut"  portion  of  a  round  may  be 
patterned  after  explosive  tunnel  rounds,  which  use  draw,  V,  pyramid, 
burn  and  other  types  of  cut*. 

Slot  Round 

In  the  first  experimental  field  test  at  a  granite  quarry  two  per¬ 
pendicular  slots  2  inches  wide  by  30  inches  deep  were  cut  In  the  granite 
in  the  center  of  the  tunnel  face  with  a  pneumatic  (drifter)  drill.  Such 
a  slot  requires  about  five  times  as  much  drilling  as  an  explosive  round. 
At  an  estimated  cost  of  30d/foot  of  drill  hole  the  total  cost  of  drill¬ 
ing  for  a  15  by  15  foot  tunnel  216  feet  of  drill  hole  is  required  or 
$64.80/t'oot  of  tunnel.  Assumimji  for  purposes  of  illustration  that  half 
of  the  work  time  is  used  for  drilling  at  a  r'ate  of  24  Inches/min,  15 
drills  would  be  required  for  100  foot/day  tunnel  advance.  The  total 
heat  energy  for  breakage  is  estimated  at  a  maximum  of  2  kwhr  per  foot 
length  of  heater  hole  (36  holes)  to  be  72  kwhr.  At  3<t/kwhr  the  energy 
cost  for  breakage  would  be  $2. 16/foot  of  tunnel. 

The  initial  heater  hole  size  required  is  determined  largely  by  the 
space  required  for  insertion  of  the  heater.  The  most  effective  tempera¬ 
ture  for  inducing  stresses  in  the  rock  is  theoretically  the  highest 
that  can  be  applied,  because  higher  temperatures  result  in  greater 
heat  flow  and  sharper  stress  increases  at  large  distances.  Carbon 
electrodes  first  create  an  arc  (6,000“C),  causing  the  rock  to  spall, 
melt,  and  change  chemical  and  physical  form.  After  a  melt  is  formed. 
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however,  most  of  the  he«t  qertcreted  which  comes  from  the  ere  is  trans¬ 
ferred  by  several  complex  processes.  The  maximum  eff^'Ctlve  temperature 
In  the  cavity  Is  the  melting  temperature  of  the  rock,  which  Is  usually 
lower  than  1700**C  for  silicate  rocks.  With  any  type  of  heater  that 
causes  the  rock  to  melt  the  cavity  Is  being  continually  enlarged,  which 
Is  advantageous.  That  Is,  the  larger  the  cavity  the  greater  the  stresses 
are  at  larger  distances  from  the  cavity. 

A  hole  spacing  of  2  feet,  which  may  be  about  optimum  for  12  kw  car¬ 
bon  electrode  heaters.  Is  similar  to  hole  spacing  required  tor  explo¬ 
sives.  Explosive  energy  equivalent  varies  from  0.4  to  0.6  kwhr/lb  for 
40  percent  to  75  percent  dynamites,  with  about  1  1b  per  ft  of  hole  for 
strong  rocks. 

While  electrical  breakage  energy  requirements,  based  upon  tenta¬ 
tive  Information,  are  about  three  times  greater  for  thermal  fragmenta¬ 
tion  than  for  explosives,  electricity  costs  about  2  to  3</kwhr,  dyna¬ 
mites  at  50«t/lb  range  fro(n  $0.83  to  $1.25/kwhr.  Most  Inexpensive  ammonium 
nitrate-fuel  explosives  are  not  suitable  for  breaking  hard  rocks. 

The  time  required  for  automated  insertion  of  heaters  such  as  carbon 
electrodes  1n  a  hole  Is  a  very  small  fraction  of  that  required  for  load¬ 
ing  explosives.  More  Important  is  that  the  operational  characteristics 
of  the  proposed  system  will  permit  simultaneous  drilling,  fragmentation 
and  mucking.  That  1s,  the  operation  as  planned  will  be  semi -continuous 
or  continuous. 

Although  for  a  single  slot  of  2-foot  length  it  was  not  possible 
to  obtain  effective  fracture  of  rock  on  the  sides  of  the  slot,  It 
may  be  that  a  single  long  slot  in  the  full  height  of  the  face  will 
be  sufficient  to  provide  displacement  relief  for  breakage. 
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At  the  present  stage  of  devoloptnent  the  themiel  round  concept 
appears  to  have  few  of  the  disadvantages  of  the  other  novel  concepts, 
and  has  promise  of  overcotning  most  of  the  disadvantages  of  the  drill 
and  blast  method.  While  there  are  several  basic  and  operational  prob 
lenis  to  be  solved,  the  state-of-the-art  Is  sufficiently  well  advanced 
In  the  areas  Involved  to  permit  an  optimistic  prediction  of  early  ap- 
pl Icatlon. 

No  radically  new  types  of  equipment  will  be  required.  Current 
types  of  drill  assemblies  can  be  adapted  for  slot  and  heater  hole 
drilling.  Heater  elements  (carbon  electrodes)  are  currently  avail¬ 
able  at  very  low  cost,  and  will  require  only  minor  Improvements  and 
the  development  of  mobile  holders.  Capital  equipment  and  operation 
costs  appear  to  be  In  the  same  range  as  drill  and  blast,  but  higher 
rates  of  advance  are  Indicated. 

Important  characteristics  may  be  summarized  as  follows: 
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Characteristics  of  Thermal  Fragmentation  Method  of  ExcavaUon  "  Slot  Cut 
Drilling  required:  1  ft  of  holc/cu  ft  of  rock  removed 

Drilling  energy:  200  to  390  j/cc 

Perct nt  of  face  drilled:  5% 

Breakage  energy:  0.2  kwhr/cu  ft 

Rock  melted:  maximum  of  0.3?'> 

Labor  at  face:  4  men  In  15  x  15  tunnel 

Character  of  muck:  Small  fragments  to  3  ft  cubes 

Continuity:  Semicycllc  to  continuous 

Down  time:  Minimal 

Heat  developed:  Negligible 

Noise  and  vibration:  From  pneumatic  drilling  only 

Gases  created:  Small  amounts  -  easily  removed 

Environmental  problems:  Minimal 

Equipment  Costs:  Nominal,  drills,  jumbos,  etc.,  costs  for 

modification  and  redesign 
Services:  Electrical  power 

Compressed  air,  assists  ventilation 

Tunnel  wall  stability:  Some  slight  fracturing  of  surrounding  rock 

Equipment  required:  Standard  drilling  equipment  -  modified 

Mountings  for  electrodes 
Improved  electrodes 
Estimated  time  to  prototype:  3  years 
Estimated  time  to  full  operation:  4-6  years 
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Cost  Analysis 

The  direct  costs  for  excavation  are  usually  considered  under  the 
categories  of  fixed  costs,  or  those  which  are  constant  per  foot  of 
advance,  and  variable  costa,  which  are  constant  per  day  charges,  and 
consequently  vary  with  the  rate  of  advance.  Examples  of  fixed  costs  are 
power  and  bit  wear,  and  variable  costs  are  equipment  and  labor 
costs.  Overhead  costs  also  fall  Into  the  last  category. 

It  Is  not  possible  to  make  a  simple  comparison  by  category  of 
the  direct  costs  of  rock  breaking  or  i^emoval  by  such  methods  as  tunnel 
boring  machines,  drill  and  blast,  and  excavation  utilizing  melting 
processes,  because  of  differences  In  the  operational  characteristics 
of  each.  Also,  each  method  has  different  types  of  trade-off  costs. 

For  example,  the  smooth  walls  created  by  tunnel  boring  machines  may 
reduce  the  cost  of  tunnel  support,  and  the  glass  liners  created  by 
melting  may  have  the  same  value.  One  method  may  require  less  venti¬ 
lation  than  another,  or  incur  higher  muck  removal  costs. 

Hence,  a  complete  comparison  of  methods  would  need  to  Include  cost 
factors  In  all  categories  which  are  affected  by  the  overall  operation. 

Costs  for  Drilling  and  Rock  Breaking.  In  a  study  of  excavation  costs  In 
the  NE  Corridor  by  Harza  Engineering  Company  In  1968  (84)  rock  breaking 
costs  were  given  for  four  general  types  of  rock: 

X.  Rocks  assumed  to  be  mineable  with  boring  machines  with  re¬ 
lative  drill  ability  factors  between  5  and  7,  average  unconfined  com¬ 
pressive  strength  4000  to  32,000  psi. 

Y.  Same  as  X  with  higher  average  compressive  strengths,  18,000  to 
20,000  psi,  65  percent  of  rock  type  is  mineable  with  boring  machines. 
Relative  drill  ability  factor  of  1  to  4. 
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Z.  Rocks  of  somewhat  the  same  strength  as  Y,  but  requiring  con¬ 
ventional  drill  and  blast  methods. 

D.G.  Difficult  ground  associated  with  faults,  gouge,  susceptible 
to  squeezing. 

The  representative  direct  rock  breaking  costs  were  categorized  as 
follows  by  percentage  distribution 


Unit  X 

Y 

Z 

DG 

Labor 

33 

39 

50 

30 

Equipment 

62 

54 

38 

65 

Materials 

5 

7 

12 

5 

It  is  notable  that  labor  costs  are  the  highest  for  drill  and  blast  (Z) 
partly  due  to  the  cyclic  nature  of  this  method.  It  is  notable  also 
that  costs  per  foot  increase  almost  linearly  with  the  excavated  diameter 
(Fig.  79). 

Because  of  the  high  labor  costs  of  drill  and  blast,  novel  methods 
of  breaking  hard  rock  have  a  considerable  margin  of  capital  which  can 
be  expended  for  additional  energy  and  equipment  if  the  method  is  automated 
and  has  other  requisite  cost  advantages  similar  to  boring  machines. 

Thus,  for  a  tunnel  driven  by  a  thermal  fragmentation  method,  possible 
trade-off  for  increased  drilling  costs  will  be  labor  reduction,  reduction 
in  cost  of  breakage  energy,  higher  rate  of  advance,  and  other  factors 
yet  to  be  defined. 

In  analyses  of  pneumatic  drilling  costs  it  has  been  found  that  the 
factors  of  greatest  cost  are  for  labor,  equipment,  and  maintenance  (85). 
Amortization  or  rental  charges  per  foot  for  equipment  are  less  for  two 
or  three  shift  operation. 
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For  a  single  2"  pneumatic  drill  the  costs  are  approximately  as 


follows : 

Per  hr 

Compressor 

$  3.75 

Bit  Cost 

$25.00 

Drill 

4.25 

Life 

2,000 

Labor  (2  men) 

7.00 

Reconditioning 

$15.00 

Fuel  (4  gal) 

0.60 

Rod  cost/ ft 

$  0.03 

2"  carbide  bit  & 

rod 

3.00 

$18.60 


Rate  of  drilling  in  hard  granite  (overall)  60  ft/hr 

Cost/ft  =  $0.35 

A  method  of  excavation  employing  thermal  fragmentation  may  be 
semicyclic  or  continuous.  One  procedure  would  permit  simultaneous 
drilling  of  heater  holes  and  slots,  or  another  the  cutting  of  a  slot 
first  and  simultaneous  drilling  of  heater  holes  and  fragmentation. 

The  latter  would  probably  require  dry  drilling.  Also,  after  the  rock  is 
fractured,  rock  removal  is  accomplished  by  rock  splitter  or  impact 
tool.  Direct  costs  may  be  estimated  as  follows: 

Example:  10  x  10  tunnel 
Hole  spacing:  2  ft  centers 
Depth  of  round:  2  ft 

Number  of  heater  holes:  36  -  2-1/4"  -  72  ft  total 

Slots:  2  -  2"  X  2'  X  10'  240  ft  of  2"  diam  hole 

16  ft/shift  =  576+ 

1920 

ft  of  hole 


Advance: 
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Equipment  Costs: 

Cost/hr 

Total/hr 

3  Shifts 

10  small  drifter  drills  150  cfm 

3.50  ea. 

35.00 

11.67 

4  large  drifter  drills  400  cfm 

4.75  ea. 

19.00 

6.33 

1  compressor  3000  cfm 

4.50/drill 

63.00 

21.00 

Jumbo 

5.00 

5.00 

1.66 

Fuel  (compressor  5  gal/drill 

10.50 

10.50 

10.50 

Rods 

0.40 

5.60 

5.60 

Bits  .40/drill 

0.40 

5.60 

5.60 

Maintenance  &  Misc. 

0.50 

2.00 

2.00 

Labor  (3  men) 

15.00 

15.00 

160.70 

15.00 

79.36 

Cost  per  ft  of  hole 

0.50 

0.25 

Drilling  costs/ft  of  tunnel 

$40.00 

Drilling  costs  would  be  reduced  for  other  types  of  cuts  under  investigation. 
Electrical  energy  costs  are  given  under  Thermal  Round  Concepts. 

Other  factors  which  are  critical  in  the  design  of  a  thermal  round 
are  the  hole  placement  with  respect  to  geologic  structure,  effect  of 
water  in  small  and  large  amounts.,  tip  most  effective  type  of  strain 
relief  that  can  be  obtained  ■ 

From  fracture  phenomerio  '  '  vp.i  itte  the  investigators  feel 
that  heater  holes  can  be  placed  fliiv.-ini .^(jeous  ly  in  relation  to  existing 
cracks.  Fractures  parallel  to  a  heater  hole  relieves  the  stresses 
locally  but  are  exploited  by  the  expansion  of  the  rock  around  the  cavity. 

In  granites  the  joint  systems  are  seldom  orthogonal  and  may  present 
problems  in  obtaining  the  most  effective  use  of  stress  fields. 
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Thus,  the  effects  of  hole  placement  and  direction  with  respect  to  natural 
joint  patterns  in  the  rock  will  be  investigated  with  the  objective  in 
mind  of  utilizing  the  joints  in  aiding  the  thermal  fracture  process. 

The  overall  effects  of  moisture  on  heaters,  heat  coupling,  and 
relat**d  factors  are  not  known.  However,  for  small  amounts  of  water, 
few  problems  are  expected.  Where  large  amounts  of  water  are  present, 
loss  of  heat  for  water  vaporization  will  be  minor  unless  water  is  flowing 
into  a  heater  hole.  On  the  other  hand,  1t  may  be  necessary  to  protect 
high  temperature  portions  of  the  heaters,  refractory  metals  or  carbon 
electrodes,  from  sudden  cooling  and  consequent  fracture  due  to  thermal 
shock. 

To  date,  crossed  slots  in  the  working  face  appear  to  be  one 
effective  means  for  providing  displacement  relief  for  rock  fracture.  Ar¬ 
rangements  of  heater  and  splitter  holes  similar  to  those  in  a  V-cut, 
a  pyramid  cut,  or  a  burn-Cit  utilized  in  explosive  round  offer  real 
promise.  These  configuration^,  are  being  analyzed  and  investigated 
experimentally. 


I 


190 


Rate  of  Advance 

Observations  of  pneumatic  drilling  operations  as  an  Integral  part  of 
a  fragmentation  system  Indicate  that  for  an  optimum  number  of  drills  in 
a  tunnel  heading,  the  rate  of  advance  will  be  determined  by  the  absolute 
and  overall  rate  of  drilling.  For  a  larger  number  of  drills,  the  rate 
of  advance  will  be  determined  by  a  balance  of  the  number  of  drills  and 
heaters  which  can  be  efficiently  and  economically  used  In  a  heading  of 
a  given  size  and  shape. 

The  overall  rate  of  drilling  Is  fixed  by  the  absolute  rate  of  drill¬ 
ing  adjusted  for  the  time  required  to  position  the  jumbo,  the  drills, 
and  for  delays,  i.e.,  for  effective  drilling  time.  A  further  adjustment 
must  be  made  for  the  time  required  for  drilling  and  breaking  the  cut 
portion  of  a  round.  Thus,  a  drilling  time  factor,  the  percent  of  the 
total  working  time  at  the  face  during  which  the  drill  is  actually  opera¬ 
ting,  will  govern  the  rate  of  advance. 

In  the  current  experiments  in  Missouri  red  granite  a  2-3/4  Inch  dia¬ 
meter  hole  Is  drilled  for  heater  holes  and  1-3/4  inch  diameter  hole 
for  splitters.  The  absolute  rate  drilling  Is 

2-3/4  inch:  10  In. /min 
1-3/4  Inch:  20  in. /min 

The  following  example  calculations  can  be  made  for  seven  and  nine  drills 


per  heading  in  a  20  ft  diameter  tunnel: 

Tunnel  size:  20  ft  diam  20  ft  diam 

Heater  hole  spacing:  2.5  5.0 

No.  of  heater  holes:  81  65 

No.  of  splitter  holes:  40  34 

No.  of  drills 

large:  9  7 

small:  3  2 


Max. 

rate  of  adv.  ft/day: 

150 

150 

75^.  d.t.f.: 

113 

113 

50'/.  d.t.f.: 

75 

75 

No. 

of  heaters: 

8 

8 

with  a  heater  hole  spacing  of  ?.5  ft  Ideally  3  x  3  x  2  ft  block 
would  be  broken.  For  maximum  advance  at  2  ft/round,  this  would  require 
1.8  blocks  per  minute.  With  four  sets  of  four  heaters  each,  and  a 
fracture  time  of  three  minutes  per  block,  and  simultaneous  heating  of 
two  or  more  blocks,  the  required  advance  rate  would  be  achieved.  With 
three  or  four  heat  sets  operating  simultaneously,  the  power  level  would 
be  approximately  150  kw  at  12.5  kw  per  heater.  Rates  at  75^.  or  50K  drill 
time  factor  appear  to  be  readily  attainable. 

With  automated  jumbos  and  manual  splitter  and  heater  placem\ent,  and 
efficient  muck  removal,  a  high  drilling  time  factor  can  be  obtained.  It 
can  be  further  improved  If  holes  normal  to  the  face  can  be  drilled  for 
more  than  one  round  from  one  drill  setup. 

It  is  estimated  that  drill  time  factors  will  vary  from  about  50 
to  75  percent,  which  would  give  a  rate  of  advance  varying  from  75  to 
113  ft/day  in  20  ft  diameter  tunnels  for  2.5  and  3  ft  spacing.  Th*. 

Is  1.5  to  2.8  times  the  rate  of  advance  In  similar  hard  rocks  which 
require  drill  and  blast  methods,  as  reported  by  Harza  (Fig.  80)> 
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A  close  correUtlon  between  the  theorotlcAlly  previicted  enti 
the  iKne»'1mentrtHy  obse>'ved  fracture  patterns  and  the  fracture  lenqth- 
fracture  t<ine  relationships  Indicate  that  the  conclusloni  drawn  re- 
gord1?>g  the  Influence  of  various  process  parameters  using  the  two  two- 
dimensional  plane  models  can  be  applied  fairly  accurately  to  a  three- 
dimensional  configuration.  Also,  for  rock  types  with  characterU-  > 
similar  to  Dresser  basalt,  the  approximate  fracture  completion  time 
corresponding  to  any  given  fracture  length  can  be  approximated  frw« 
the  data  for  a  single  test  using  the  power  relationship 
between  the  dimensionless  fracture  length. 

For  a  given  geometry,  the  predicted  fracture  time  Is  only  a  small 
fraction  of  the  time  necessary  to  reach  a  steady-state  temperature 
distribution.  The  temperature  field  at  the  time  of  the  fracture  Is 
thus  highly  localized,  and  the  major  portion  of  the  fractured  volume 
experiences  only  a  small  Increase  In  temperature.  The  fracture- 
inducing  stresses  In  this  region  are,  therefore,  entirely  due  to  the 
load  vector  resulting  from  the  thennal  constraint.  Since  the  severity 
of  thermal  constraint  depends  on  the  relative  volumes  of  the  cold  and 
the  hot  zones,  the  fracture  time  is  significantly  influenced  by  both 
the  hole  spacing  and  the  melt-free  depth. 

For  very  small  melt-free  depths,  the  fracture  1s  characterized  by 
mainly  a  compressive  failure  according  to  the  McCl intock-Walsh  modi¬ 
fied  SrlffUh  criterion.  This  fracture  mode  Is  naturally  undesirable 
as  the  compressive  strength  of  hard  rocks  is  many  times  greater  than 
their  tensile  strength. 


\n 

WhiW  the  McCt iMtotk-WeUh  I va  felUire  /ona  liNtwinii  r^tn* 

ttvely  uM(-haM\}e(t  tn  the  tag i ah  of  the  high  teaiiterAtgra  u>'e(MenUi  the 
aompmitve  coaiponenu  ti>  ihe  lot^e  dact'eaue  ttie«d1ty  with 

An  Incpea^a  In  the  melt-free  depth.  With  further  Intreetie  In  the  welt- 
free  depth,  A  conditinn  U  reeched  when  the  major  principal  «trfi«  cow- 
ponenttk  become  teniille,  Thuti,  with  a  continuous  Increase  In  the  welt- 
free  depth,  the  HcCl 1ntock-Wal»h  fracture  >one  In  the  cold  tniglon  trens- 
foreu  first  Into  a  pertlal  Griffith  contpresslon  jone  and  finally  Into 
the  Griffith  tension  lone.  The  effectiveness  of  th*s  transformation 
Is  evident  fr'om  the  drastic  reduction  In  the  fracture  time  as  the  frac¬ 
ture  now  resu’*s  fr««  the  tensile  failure. 

The  melt'fi'ee  depth  associated  with  the  cojnpletlon  of  the  trenstor- 
watlun  of  the  compressive  fracture  Into  the  tensile  fracture  mode  Is 
approximately  equal  to  half  the  fracture  length,  that  Is,  half  the  dif¬ 
ference  between  the  hole  specing  and  the  hole  diameter.  Any  further 
Increase  in  this  paranwtec  has  only  a  slight  effect  on  the  fracture  time 

Another  significant  effect  of  the  melt-free  d'pth  Involves  the  loca¬ 
tion  of  the  parallel  fracture  (subsurface  fracture  parallel  to  the  work¬ 
ing  face).  The  optimum  location  of  this  fracture  requires  the  crack 
Initiation  to  occur  very  close  to  the  hole  base  and  the  crack  propaga¬ 
tion  In  a  plane  approximately  parallel  to  the  hole  base,  from  the  slot 
model  analysis,  to  obtain  such  a  location,  It  is  necessary  that  the 
thermal  Inclusions  should  be  concentrated  at  the  very  base  of  the  holes 
(requiring  only  a  small  melt  dopth)  while  the  melt-free  depth  should, 
at  least,  be  equal  to  half  the  fracture  length. 

The  hole  model  analysis  I*  performed  using  a  plane  stress  assump¬ 
tion.  The  actual  three-dimensional  fragmeni.it ion  configuration,  however 
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Itwolvt^s  A  t’oodlilon  tho  pUn^  (vhA))ow 

hoUik  h»At»(t  Alonp  tht»li'  unth'o  Unqihi)  Anil  th^  plAn«  i^trAlo  ((it»A|> 
hQl«»  hi>At9(1  Al«>nu  iheir  lengths)  forwulAllons ,  AUu.  thU 

AMAlysU  (10A«  ml  involve  the  m«lt>f«'«e  depth  pArAineter  which 
thown  fiHW  the  Hlot  model  titiidle&  cao  be  u^ed  AdVAntAqeoiisl v  to  trtn^" 
fonn  the  mode  of  feilure  from  cceopresiive  to  tensile  And  thereby  to 
yrtAtly  reduce  the  frecture  time,  Thus»  it  is  naturAl  to  expeot  a 
precioAiinAotly  compressive  feilure  in  the  cose  of  the  hole  models, 
the  results  of  which,  in  view  of  the  Above  consideretions ,  ere  ob¬ 
viously  representAtive  of  the  upper  bound  of  the  theoreticdl  solu¬ 
tion. 

The  results  of  the  pUne  Analysis  indicate  that  since  the  fracture 
initiates  !•'  the  vicinity  of  the  point  of  transition  from  the  melt  con¬ 
dition  to  the  convection  condition  (alonq  the  hole  depth),  only  small 
melt  depths  should  be  used.  Also,  in  order  to  optimiie  the  location 
of  the  parallel  fracture  and  the  fracture  time,  the  melt-free  depth 
should,  at  least,  be  equal  to  half  the  difference  between  the  hole 
spacing  and  the  hole  diameter.  Although  any  further  Increase  in  the 
melt-free  depth  shows  an  insignificant  influence  on  the  location  and 
time  of  the  parallel  fracture,  it  will  mean  that  the  heat  source  will 
have  a  greater  burden  against  which  to  open  cracks  between  the  holes. 
The  optimum  depth  therefore  seems  to  be  one  associated  with  a  melt-free 
depth  equal  to  approximately  half  the  difference  between  the  hole  spac¬ 
ing  and  the  hole  diameter. 

Both  the  slot  and  the  hole  model  analyses  Indicate  the  hole  spac¬ 
ing  to  be  the  most  influential  parameter  governing  the  fracture  time. 
The  effects  of  the  changes  in  the  hole  diameter,  on  the  other  hand,  are 
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only  in  nml  cun  hut  tAken  into  Account  by  Adjust- 

inq  tht>  hoU»  HOAiiny  such  thAt  th«  fcACture  length  r'emrtins  unchenged. 
Thus,  the  effects  of  the  hole  speclng  And  the  hole  dlAineter  ten  br 
enpressed  In  tei^is  of  a  single  vAriAble,  the  fracture  length. 

In  practice,  the  hole  diameter  will  be  determined  mainly  front 
heater  siye  considerations,  and  since  the  optimum  hole  depth  is  ex¬ 
pressed  In  terms  of  the  fracture  length,  the  thennal  fragmentation 
configuration  can  bo  optimiied  by  a  proper  choice  of  the  single  para¬ 
meter,  the  hole  spacing. 

While  the  posver  relationship,  t!^  L*'*’  between  the  dimension¬ 

less  fracture  time  and  the  dimensionless  fracture  length  (for  Dresser 
basalt)  implies  faster  friigmentation  for  smaller  spacings,  it  is  im¬ 
portant  to  note  that  smaller  spacing  also  means  higher  drilling  costs 
as  well  as  smaller  rock  volume  that  can  be  removed  through  each  cycle 
of  operation.  The  choice  of  an  optimum  spacing,  therefore,  will  also 
involve  factors  such  as  the  drilling  rate,  the  cost  of  heating  and 
drilling,  the  mucn  removal  capacity  and  the  overall  efficiency  of  the 
actual  excavation  machine. 

Thus,  while  the  analysis  presented  here  is  not  sufficient  to  en¬ 
able  one  to  optimally  choose  the  process  parameters,  It  serves  two 
very  useful  purposes.  First,  it  reduces  the  optimal  choice  of  various 
process  parameters  to  the  choice  of  a  single  parameter,  the  hole  spac~ 
Ing,  and  also  provides  a  relationship  between  the  fracture  time  and 
the  hole  spacing.  Secondly,  since  this  investigation  considers  frac¬ 
tures  due  entirely  to  the  effects  of  thermal  inclusions,  the  results 
presented  here  can  be  used  with  advantage  to  approximately  determine 
the  effectiveness  of  the  mechanical  loading  to  be  incorporated  in  de¬ 
signing  an  excavation  prototype. 
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(;xp«r1m®nt<ition  and  Analysis 

An  Investigation  of  the  technical  and  nconpailc  feasibility  of 
thermal -mechanical  fragmentation  of  hard  rock  is  being  carried  out 
In  the  field  to  advance  the  research  toward  an  early  proof  of  concept. 
Research  has  advanced  through  phases  of  surface  heating  and  chipping, 
deep  slot  cutting  and  kerf  removal,  to  hole  drilling  and  Insertion 
of  heaters  Inside  the  rock.  Colled  wire  heaters  successfully  frac¬ 
tured  granite  blocks,  but  carbon  arcs  were  more  effective. 

Carbon  arc  heaters  are  Inexpensive,  conslsti  ig  of  two  carbon  rods 
secured  to  asbestos  board.  These  furnish  efficient  heat  sources  up 
to  a  12  kw  level  for  5/8  in.  diameter  carbon  rods.  Problems  encountered 
with  the  arc  heaters  If’clude  arc  stability,  electrode  erosion,  restart¬ 
ing  arcs  In  a  solidified  melt,  and  brittleness  of  the  carbon.  All  of 
these  difficulties  appear  to  be  amenable  to  correction  with  proper 
heater  design  and  application  of  arc  engineering  techniques. 

Blocks  of  rock  can  be  readily  fractured  along  orthogonal  planes 
If  three  directions  of  displacement  relief  are  provided.  That  1s, 
four  thermal  inclusions  generate  a  fracture  plane  at  depth  and  parallel 
to  the  working  face  as  well  as  longitudinal  fractures  between  the  drill 
holes,  forming  rectangular  blocks  of  rock.  The  temperature  field  creates 
Intense  stresses,  but  little  displacement  without  prolonged  heating. 
Hence,  Darda  rock  splitters  are  employed  to  enhance  and  exploit  the 
early  thermal  fractures,  and  to  remove  the  blocks,  which  are  often 
keyed  in  place.  This  type  of  fracturing  lends  itself  to  the  design 
of  a  thermal  round,  similar  to  a  blasting  round,  with  heat  employed 
instead  of  explosives.  As  with  a  blasting  round  a  critical  factor  is 
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the  tieslijn  «nd  proper  functioning  of  the  cut  portion  of  the  round  In 
such  a  manner  that  adequate  din placement  relief  Is  provided  for  sub¬ 
sequent  fracture  and  removal  of  blocks  of  rock  from  the  remainder  of 
the  working  face. 

Several  types  of  cuts  will  be  Investigated  In  future  research  such 
as  pyramid,  V,  draw,  simulated  burn  cuts,  as  well  as  a  breakage  process 
employing  a  continuous  spiral  drill  and  fracture  process. 

Preliminary  economic  studies  Indicate  that  costs  of  excavation  should 
be  approximately  the  same  as  for  conventional  drill  and  blast  methods, 
and  that  projected  rates  of  advance  will  be  two  to  three  times  as  gr'iat. 
Much  of  the  drilling  equipment  required  can  be  obtained  by  minor  modi¬ 
fications  of  available  Items.  About  the  same  amount  of  drilling,  or 
slightly  more,  may  be  required  as  for  drill  and  blast  procedures.  Con¬ 
siderable  heater  development  may  be  required  for  a  type  of  element 
needed  for  continuous  use  underground.  Breakage  energy  costs  are  about 
three  times  as  great  as  for  explosives,  but  electrical  energy  costs  are 
only  one-third  those  of  explosives.  Mobile,  multiple  purpose  jumbos 
will  be  required  for  mounting  drills,  heaters,  and  rock  splitters. 

A  different  type  of  mucking  system  may  be  required  for  handling  and 
removal  of  large  blocks  of  rock.  Shields  will  be  needed  at  the  face 
so  that  drilling  and  breaking  can  be  carried  on  simultaneously. 

To  date  experimentation  and  analysis  indicate  that  this  method  of 
fragmentation  of  both  hard  and  soft  rock  for  tunneling  and  other  types 
of  excavation  is  both  technicallv  i^nd  economically  feasible. 
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